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ABSTRACT 


The  external  nucleate  pool-boiling  heat-transfer 
coefficient  of  both  smooth  and  enhanced  horizontal  tubes 
in  R-114-oil  mixtures  (0  to  10  mass  percent  oil)  were 
measured  for  heat  fluxes  from  1  to  91  kW/m2  at  a  saturation 
temperature  of  2.2  C.  The  enhanced  tubes  tested  were  a 
Wieland  GEWA-T  finned  tube  containing  1.02  fins/mm,  a 
Hitachi  Thermoexce 1 -E  and  a  Hitachi  Thermoexce 1 -HE  tube. 
The  Thermoexce 1 -E  and  -HE  tubes  with  their  re-entrant 
cavity  designs  were  found  to  improve  the  heat-transfer 
coefficient  over  the  smooth  tube  value  at  a  constant  heat 
flux  by  a  factor  of  approximately  7  in  oil-free  R-114. 
while  the  GEUA-T  tube  improved  the  coefficient  bv  a  factor 
of  about  4.  Whi le  al 1  of  the  tubes  showed  a  general ly 
decreasing  performance  with  the  presence  of  oil.  the  GEWA-T 
tube  resulted  in  the  minimum  reduction.  This  tube  showed 
up  to  20%  and  35%  reduction  (compared  to  the  oil-free  case) 
at  3  percent  and  10  percent  oil,  respectively.  The 
Thermoexce 1 -E  and  Thermoexce 1 -HE  tubes  showed  performance 
reductions  of  up  to  40%  and  80%,  respectively,  with  3 
percent  oil.  At  a  practical  heat  flux  of  about  30  kW/m2 
with  3  percent  oil,  the  GEWA-T,  Thermoexce 1 -E  and 
Thermoexce 1 -HE  tubes  outperformed  the  smooth  tube  by 
factors  of  4.8,  4.6  and  4.0,  respectively. 
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L_  INTRODUCTION 

A.    BACKGROUND 

In  the  design  of  refrigeration  systems  for  use  aboard 
Naval  vessels,  the  desire  is  to  make  the  svstem  as 
efficient  as  possible.  To  achieve  this  higher  efficiency, 
one  must  consider  the  use  of  a  high  pressure  refrigerant. 
since  this  gives  a  higher  energy  transfer  per  unit  volume 
of  vapor  (i.e..  Ev  =  h*  3 /v,  )  CRef.  13,  and  results  in  less 
of  a  pressure  drop  along  the  vapor  path  for  a  given  heat 
load.  R-12,  which  is  a  high  pressure  refrigerant,  has  an 
energy  transfer  per  unit  volume  of  vapor  of  8570  kJ/'m3 
CRef.  23  and,  since  it  is  higher  pressure,  it  requires 
heavier  gauge  components.  R-ll,  on  the  other  hand,  is  a 
low  pressure  refrigerant  with  a  low  energy  transfer  per 
unit  volume  of  vapor  of  2086  kJ/m3  and  requires  lighter 
gauge  mater  i  a  1 . 

Table  1  provides  a  short  summary  of  some  of  the 
available  refrigerants  and  their  respective  energy  transfer 
per  unit  volume  of  vapor,  their  relative  toxicity  and  their 
required  system  pressure  for  optimum  operation.  As  can  be 
seen  from  Table  1,  R-12  and  R-22  possess  the  highest  values 
of  energy  transfer  per  unit  volume  of  vapor,  and  therefore 
would  provide  systems  of  higher  efficiency  than  the  other 
refrigerants  listed  since  pressure  losses  in  the  system 
would    be    less.    However,   since  both   R-12  and  R-22  are 

11 


hi gh- pressure   refrigerants,  their   respective  systems  would 

be   significantly    heavier  larger   than  those   of  the  other 

refrigerants.   However,   R-114   being   a   moderate-pressure 

refrigerant   with  an   energy  transfer   per  unit   voulume   of 

vapor   of  3800   kJ/m3,  it   would    allow   for   a   relatively 

lighter   component   system   than   those  of   R-12   and  R-22. 

TABLE  1  CRef.  31 
SUMMARY  OF  REFRIGERANTS 

Operat  i  ona 1 
System 
Ref  ri  gerant _ kJ/m3 Toxicity _ Pressure 

R-ll  2086  moderate  (Group  5)          low 

R-12  8570        low  (Group  6)  high 

R-22  12480  moderate  (Group  5)  high 

R-113  1128       high  (Group  4)          low 

R-114  3800  moderate  (Group  6)  moderate 

Additionally,   R-114  belongs  to  the  refrigerant  group  (Group 

6)    having   the  lowest   toxicity  CReF.   2],  and  it  has  also 

been    proven   fairly   stable   with   temperature   [Ref.   21. 

Therefore,   R-114  is   a  viable   alternative   to   the   higher 

pressure  refrigerants. 

B.     REVIEW  OF  REFRIGERANT-OIL  MIXTURE  BOILING  BEHAVIOR 

The  presence  of  oil  in  the  evaporators  of 
refrigeration  systems  is  a  common  occurrence.  This  is  due. 
in  part,  to  the  use  of  lubricants  in  the  hermetically 
sealed  refrigeration  compressors.  Over  a  period  of  time,  a 
considerable  amount  of  oil  is  introduced  into  the 
refrigerant,  subsequently  changing  the  physical 
properties   of  the  refrigerant.   Jensen  and  Jackman  [Ref.  3] 


12 


reported  that  the  density  and  the  specific  heat  behave 
linearly  in  the  refrigerant-oil  mixtures,  however,  the 
viscosity  and  the  surface  tension  do  not.  Henrici  and 
Hesse  [Ref.  4]  experimentally  determined  in  1971  that  the 
surface  tension  of  mixtures  first  decreased  for  R-114-oil 
mixtures  of  up  to  2.5  percent  oil,  and  at  higher  oil 
concentrations,  the  surface  tension  increased.  This  non- 
linear behavior  of  the  properties  of  the  mixtures  makes  it 
difficult  to  explain  the  changes  in  the  heat- transfer 
coefficient  for  these  mixtures. 

The  most  evident  example  of  the  addition  of  oil  to  the 
refrigerant  is  the  presence  of  foam  in  the  mixture  when 
boiling  takes  place.  With  the  oil  concentration  above  1 
percent,  there  is  a  substantial  amount  of  foam  generated 
from  the  nucleate  boiling  surface.  This  foam  is  generated 
because  the  refrigerant  in  the  ref r i gerant -o i  1  mixture  is 
more  volatile  than  the  oil  and  therefore  can  vaporize, 
creating  a  gas  bubble  surrounded  by  an  oil-rich  layer 
[Ref.  53.  Figure  1.1  shows  the  ideal  model  of  the  growth  of 
a  bubble  in  a  refrigerant-oil  mixture.  These  generated 
bubbles  are  at  lower  density  than  that  of  the  surrounding 
refrigerant-oil  mixture.  The  bubbles  ascend  to  the  liquid- 
vapor  interface  and  collect  at  the  interface  and  produce  a 
foam  layer.  As  reported  by  Henrici  and  Hesse  [Ref.  4],  the 
foaming    action   is    most   pronounced   in   refrigerant-oil 
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Figure  1.1      Oil  Concentration  Gradient  in  a  Bubble  in 

Refrigerant-Oil  Mixtures  CRef.  53. 
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mixtures  of  1  to  10  percent  oil,  and  provide  a  the 
presence  of  oil  results  in'  a  subsequent  decrease  in  the 
heat- transfer  performance   of  the  boiling  tube. 

Numerous  researchers  have  tried  to  explain  the  reasons 
for  the  decrease  in  the  heat- transfer  coefficient  in 
refrigerant-oil  mixtures  of  increasing  oil  content.  Thome 
CRef.  63,  after  an  extensive  literature  search,  found  that, 
in  1956,  Van  Wijk  et  al.  explained  the  effect  as  being  the 
result  of  the  evaporation  of  the  more- vo 1  at i 1 e  component 
leaving  an  oil-rich  layer  which  has  a  higher  resultant 
local  boiling  point.  This  increased  local  boiling  point 
requires  an  increase  in  the  amount  of  superheat  necessary 
to  continue  the  vaporization  and  bubble  growth.  The  result 
is  a  decrease  in  the  heat- transfer  coefficient. 

Stephan  and  Preusser  [Ref.  73  determined  that  the  work 
required  to  form  bubbles  in  the  refrigerant-oil  mixtures 
was  greater  than  in  an  equivalent  amount  of  pure 
refrigerant.  Therefore,  the  heat- transfer  coefficient  was 
lower  for  the  refrigerant-oil  mixture  than  for  the  pure 
refrigerant.  Chongrungreong  and  Sauer  [Ref.  83  reported 
that  the  rate  of  heat  diffusion,  which  is  governed  by  the 
thermal  properties  of  the  oil-rich  layer,  limits  the  bubble 
growth  and  that  the  effects  of  surface  tension  are 
negligible.  Thome  [Ref.  63  agreed  with  the  above  factors 
and   concluded  that   the  viscosity   variations  are  important 
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in      explaining    the   decrease    in    the    heat- trans f er 
coefficient   for   some   refrigerant-oil   mixtures. 

Although  researchers  have  not  decided  on  any  one  set 
of  factors  that  determine  the  reasons  for  the  decrease  in 
the  heat- transfer  coefficient  as  oil  content  increases  in 
refrigerant-oil  mixtures,  all  agreed  that  the  physical  and 
thermal  properties  are  important  factors  in  explaining  the 
heat- transfer  behavior  of  refrigerant-oil  mixtures. 

C.     REVIEW  OF  BOILING  PERFORMANCE  FROM  SMOOTH  TUBES 

The  U.  S.  Navy  attempts  to  reduce  the  size  of  its 
shipboard  refrigeration  systems  by  the  use  of  alternate 
refrigerants  (such  as  R-114).  At  the  same  time,  the  Navy 
desires  to  increase  the  performance  of  these  systems  by  the 
use  of  enhanced  surfaces.  However,  as  discussed  by 
Wann i arachchi  et  al.  CRef.  93,  only  a  limited  amount  of 
data  is  available  regarding  the  boiling  of  refrigerants  in 
the  presence  of  oil. 

Stephan  CRef.  103  conducted  experiments  using  a 
smooth  horizontal  plate  as  a  boiling  surface  in  R-12 
and  R-22  refrigerants  and  reported  up  to  50  percent 
reductions  in  the  boiling  heat- transfer  coefficient 
with  a  refrigerant-oil  mixture  of  9  percent  oil  (by 
mass),  while  50  percent  oil  resulted  in  a  90 
percent  reduction  in  the  boiling  heat- transfer  coefficient. 
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In  1971,  Henrici  and  Hesse  [Ref.  4]  conducted 
experiments  involving  the  boiling  of  R-114  in  the  presence 
of  oil  from  a  smooth  copper  tube  and  reported  a  maximum 
decrease  of  20  percent  decrease  in  the  boiling  heat- 
transfer  coefficient  with  1  percent  oil  by  mass.  In  1981, 
Stephan  and  Mitrovic  [Ref.  11]  reported  that  for  the  GEWA-T 
tube  (manufactured  by  the  Wieland  Company)  in  R-114-oil 
mixtures,  the  boiling  heat- transfer  coefficient  of  the 
enhanced  surface  was  significantly  altered  by  the  presence 
of  the  oil.  In  1985,  Reilly  [Ref.  5]  reported  a  drop  in 
the  heat- transfer  coefficient  of  0  to  35  percent  (depending 
on  the  heat  flux)  for  a  R-114-oil  mixture  of  10  percent  oil 
for  a  smooth  copper  tube. 

D.     REVIEW  OF  BOILING  PERFORMANCE  FROM  ENHANCED  SURFACES 

Arai  et  al.  [Ref.  12]  tested  a  200-ton  R-12 
centrifugal  water  chiller.  The  chiller  which  used  an 
enhanced-surf ace  tube,  commercially  known  as  Thermoexce 1 -Em 
made  by  Hitachi,  was  28  percent  shorter  in  length  and  had 
an  approximate  improvement  in  the  overall  heat- transfer 
coefficient  of  50  percent.  Additional  research  by  Yilmaz 
and  Westwater  [Ref.  13],  Marto  and  Lepere  [Ref.  14],  and 
Carnavos  [Ref.  15]  on  commercially  available  enhanced 
surface  tubes,  in  various  refrigerants  other  than  R-114, 
indicated  that  a  porous -coated  surface  exhibited  the  best 
boiling   heat- transfer  coefficient   in  a   pure   refrigerant. 


17 


Comparisons  by  Reilly  CRef.  5]  indicated  that  in  R-114, 
Union  Carbide's  "High  Flux"  tube  with  a  porous-coated 
surface  resulted  in  an  improvement  of  at  least  a  factor  of 
seven  over  a  smooth  copper  tube  at  the  same  operational 
parameters.  Reilly's  findings  complemented  the  research  of 
Yilmaz  and  Westwater  CRef.  13].  Additionally,  Reilly 
reported  that  for  a  porous-coated  surface  in  an  R-114-oil 
mixtures  of  1  percent  oil,  the  oil  was  seen  to  degrade  the 
nucleate  pool-boiling  heat- trans fer  performance  by  about 
20  percent. 

1.    Surfaces  with  Reentrant  Cavities 

Reentrant  cavities  have  been  investigated  for  a 
considerable  period  of  time.  Griffith  and  Wa 1 1 i s  [Ref.  15] 
determined  that  the  reentrant  cavity  geometry  is  important 
in  the  generation  of  nucleate  boiling  sites  in  two  ways; 
(1)  the  diameter  of  the  reentrant  cavity  mouth  determines 
the  amount  of  superheat  required  to  initiate  boiling, 
and  (2)  the  cavity  shape  determines  the  stability  once 
the  boiling  has  begun.  Webb  CRef.  16]  stated  that  the 
key  to  the  high  performance  of  the  reentrant  grooved 
structures  can  be  attributed  to  three  factors:  (1) 
reentrant  cavity  within  a  critical  size  range,  (2) 
interconnected  cavities,  and  (3)  nucleation  sites  of  a  re- 
entrant shape.  If  the  cavities  are  interconnected, 
adjacent   cavities  can   activate  each  other.   Webb  CRef.  16] 
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also   stated  that  reentrant  cavities  provide  a  stable   vapor 
trap,   which   can  remain  active  at  low  values  of  superheat. 

The  GEWA-T  (1.02  fins/mm)  enhanced-sur f ace  tube, 
manufactured  by  Wi e 1 and-Werke  AG,  is  a  copper  tube  with 
spiral  fins  of  a  T-shaped  profile.  It  has  a  reentrant- type 
cavity  formed  by  cold  working  an  integral  finned  tube.  The 
two  adjacent  T-fins  form  a  spiral  cavity  with  an  opening 
of  0.25-mm  (see  Figure  1.2). 

The  Thermoexce 1 -E  and  Thermoexce 1 -HE  enhanced-sur face 
tubes  are  made  by  the  Hitachi  Company.  Like  the  GEWA-T 
tube,  the  Thermoexce 1 -E  (Figure  1.3)  and  Thermoexce 1 -HE 
(Figure  1.4)  tubes  are  also  formed  by  cold  working  integral 
fins;  however,  these  tubes  have  low  fins  which  have  a 
small  space  cut  out  at  the  fin  tips,  giving  the  appearance 
of  a  sawtooth.  These  sawtooth  fins  are  bent  parallel  to 
the  tube  axis,  to  a  horizontal  position,  forming  tunnels 
with  evenly  spaced  surface  pores.  This  cold  working 
procedure  results  in  a  high  area  density  of  reentrant 
nucleation  sites.  The  Thermoexce 1 -HE  tube  is  a  variation 
of  the  Thermoexce 1 -E  tube,  designed  for  improved 
performance  in  1 ow-heat- f 1 ux  regions. 

Experimental  data  showing  the  effect  of  oil 
concentrations  on  the  heat  -  transfer  coefficient  of  the 
GEWA-T,  the  Thermoexce 1 -E  and  Thermoexce 1 -HE  tubes  in  R-114 
are   lacking.     In  an   effort  to   gain  more  understanding  of 
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(a)    Schematic  of  GEWA-T  Surface  Cross  Section 


UMIII. 


Cb)    Cutaway  o-p  the  n^WA-T  Surface 


(c)    Micrograph  of  the  GE^A-T  Surface 
Figure  1.2      Schematic,   Cut-away  Cross  Section  and 

Micrograph  of  the  GEWA-T  Surface. 
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Figure  1.3 


Surface  Photograph  and  Enlarged  Outer 
Surface  of  the  Thermoexce  1 -E  Surface. 
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Fi  gure  1 . 4 


Surface  Photograph  and  Enlarged  Outer 
Surface  of  the  Thermoexce 1 -HE  Surface 
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the   pool-boiling  performance   of  these   tubes,  the  David  W. 

Taylor   Naval  Ship  Research  and  Development  Center  sponsored 

this   investigation.  The   smooth  copper   tube  and  the  GEWA-T 

tube   were  provided  by  the  Wieland  Company.  The  Thermoexcel- 

E   and  Thermoexce 1 -HE   tubes  were   provided  by   the   Hitachi 

Company. 

E.    THESIS  OBJECTIVE 

The  objectives  of  this  thesis  are: 

1.  take  data  on  a  smooth  tube  with  the  use  of 
auxiliary  heaters  together  with  a  convection  shield 
to  enable  the  collection  of  more  data  with  a  minimum 
of  operator  effort, 

2.  take  boiling  data  on  a  smooth  copper  tube  in  R-114 
to  provide  a  baseline  for  follow-on  comparison 
for  enhanced  tubes  at  oil  concentrations  of  0,  1, 
2,  3,  6  and  10  percent  (by  mass)  at  a  boiling 
temperature  of  2.2  C  and  a  heat-flux  range  from  500 
W/m2  to  91  kW/m2 . 

3.  take  data  on  the  following  tubes  for  the  above- 
-mentioned  conditions: 

a.  GEWA-T  (1.02  fins/mm), 

b.  Thermoexce 1 -E,  and 

c.  Thermoexce 1 -HE. 
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II.  DESCRIPTION  OF  EXPERIMENTAL  APPARATUS 

A.    OVERALL  APPARATUS 

The  equipment  used  for  this  study  consisted  of  nine 
basic  components:  (a)  a  Pyrex-glass  tee  for  the  pool 
boiling  of  the  R-114  liquid;  (b)  a  Pyrex-glass  tee  for  the 
condensing  of  the  R-114  vapor;  (c)  an  R-114  liquid 
reservoir;  (d)  a  water -ethy 1 ene  glycol  mixture  sump;  (e)  an 
R-12  refrigeration  system  with  an  external  proportioning 
valve;  (f)  oil  reservoir;  (g)  a  vacuum  pump; 
(h)  instrumentation  and  data-acquisition  system;  and  (i)  a 
convection  shield.  Figure  2.1  shows  schematically  the 
arrangement  of  the  various  components.  Complete  details  of 
the  design,  construction  and  operation  of  the  apparatus 
are  provided  by  Karasabun  [Ref.  17],  Reilly  [Ref.  5] 
and  Wanni archch i  et  al.  [Ref.  ID.  The  original  design  of 
Karasabun  was  modified  by  the  inclusion  of  the  convection 
shield  in  the  pool-boiling  section  as  shown  in  Figures  2.2 
and  2.3. 

The  operation  of  the  apparatus  consisted  of  the 
boiling  of  R-114  liquid  in  the  lower  glass  tee  (1), 
and  condensing  of  the  R-114  vapor  in  the  upper 
glass  tee  (2).  The  resultant  condensate  was  fed  to 
the  distribution  tube  within  the  boiling  section  by 
gravity.    A   water -ethy 1 ene  glycol  mixture  from  a  30-gallon 
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Schematic  of  Experimental  Apparatus 
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Pyrex  Boiling  Tee 


Convection 
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Boiling  Tube 
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Figure    2.2 


Schematic  of  Convection  Shield  in 
Pyrex  Boiling  Tee. 


26 


(a)    Bottom  View 


(b)    Side  View 


Figure  2.3      Photograph  of  Convection  Shield 
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sump  was  used  to  cool  the  copper  condensing  coil,  thereby 
condensing  the  R-114  vapor.  This  water -ethy 1 ene  glycol 
mixture  sump  was  maintained  between  -18  and  -14  C  by  a 
one-half  ton  R-12  air-conditioning  system.  The  mixture  was 
pumped  through  the  condensing  coil  and  controlled  manually 
with  valve  VC. 

Valve  VC  controls  the  amount  of  the  wa ter-ethy 1 ene 
glycol  mixture  to  the  condensing  coil,  thereby  controlling 
the  pressure  within  the  test  loop.  For  example,  the  opening 
of  valve  VC  allows  more  of  the  water -ethy 1 ene  glycol 
mixture  to  pass  through  the  condensing  coil,  resulting  in 
additional  condensing  of  the  R-114,  and  subsequently 
lowering  the  temperature  and  pressure  within  the  entire 
system. 

As  the  subcooled  R-114  condensate  enters  the  boiling 
section  of  the  apparatus  through  the  distribution  tube,  it 
strikes  the  convection  shield.  This  shield  provides 
increased  mixing  of  the  subcooled  R-114  condensate  with  the 
remaining  R-114  (or  R-114/oil  mixture)  in  the  lower  portion 
of  the  boiling  section.  The  shield  deflects  the  incoming 
condensate  away  from  the  boiling  surface  and  allows  for 
further  mixing  of  the  condensate  and  the  existing  R-114 
liquid  by  means  of  the  auxiliary  cartridge  heaters.  This 
deflection  of  the  incoming  condensate  allowed  for  a  more 
even  temperature  distribution  of  the  liquid  R-114  and 
therefore    decreased   the  previously   reported   temperature 
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gradients  within  the  liquid  region.  This  additional  mixing 
of  the  liquid  R-114,  coupled  with  the  ability  to  use 
auxiliary  cartridge  heaters  in  the  liquid  and  allowed  for 
the  timely  achievement  of  steady-state  conditions  without 
the  previously  observed  thermal  gradients  in  the  liquid 
region  of  the  tests  done  without  the  aid  of  the  convection 
shield.  It  was  also  observed  throughout  the  operation  of 
the  apparatus  that  the  data  obtained  using  the  convection 
shield  had  a  higher  degree  of  repeatability  than  the  data 
obtained  without  the  shield. 

B.     BOILING-TUBE  CONSTRUCTION 

The  tubes  tested  during  the  course  of  this  study  were 
the  smooth  copper  tube,  the  GEWA-T  (1.02  fins/mm),  the 
Thermoexce 1 -E  and  the  Thermoexce 1 -HE  tubes.  All  of  these 
tubes  had  an  active  test  section  (i.e.,  heated)  of  length 
of  203.2-mm  (8  in),  with  an  outside  diameter  of  15.9  mm 
(5/8  in)  and  an  inside  diameter  of  12.7  mm  (1/2  in).  The 
remaining  114.3-mm  section  on  either  side  of  the  heated 
test  section  was  smooth  and  unheated  and  did  not  nucleate 
under  any  heat  flux  or  any  oil  concentration  condition. 
Figure  1.2  shows  the  schematic  of  the  GEWA-T  finned 
surface.  Karasabun  [Ref.  17]  described  the  analysis  used 
in  the  data-reduction  program,  allowing  for  the  two  end 
surfaces  to  be  treated  as  extended  fins  from  the  center 
section  and  subsequently  accounting  for  the  heat  losses. 


The  center  section  of  each  tube  tested  was  heated  by  a 
1000-W  240-V  stainless  steel  cartridge  heater  inserted  in 
the  center  section  of  each  tube.  The  cartridge  heater  was 
inserted  in  the  boiling  tube  and  surrounded  by  a  copper 
sleeve  with  eight  1.3  mm  by  1.3  mm  thermocouple  channels. 
Figure  2.4  shows  the  layout  of  the  installed  tube  in  the 
boiling  section.  Figure  2.5  shows  the  details  of  the 
thermocouple  layout.  The  thermocouple  hot  junctions  were 
inserted  into  the  copper  sleeve  and  the  copper  sleeve  was 
then  tinned;  and  heat  was  maintained  to  keep  the  solder  in 
a  molten  state  as  the  sleeve  was  inserted  and  positioned  in 
the  center  section  of  the  boiling  tube.  Upon  completion  of 
the  thermocouple  installation,  calibration  of  the 
thermocouples  was  "achieved.  Reilly  [Ref.  5]  determined  in 
the  course  of  his  investigation  that  since  the  data- 
reduction  program  utilized  the  differences  between  the 
thermocouples  in  all  computations,  such  as  wall 
temperatures  minus  saturation  temperatures,  the  corrections 
provided  by  the  thermocouple  calibration  were  basically 
unnecessary  since  the  calibration  procedures  were  only 
necessary  for  items  dependent  on  absolute  temperature. 
Appendix  A  summarizes  the  procedure  used  by  Karasabun 
[Ref.  17]  and  Reilly  [Ref.  5]  to  calibrate  the 
thermocoupl es. 
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(a)  View  of  the  boiling  tube  thermocouple  locations  as  seen 
from  the  front  of  the  experimental  apparatus. 
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(b)  Thermocouple  sleeve  unwrapped  (at  section  A-A)  to  show 
the  relative  locations  of  the. thermocouple  channels 
(all  dimensions  in   millimeters). 
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(c)  End  views  of  the  boiling  tube. 

Figure  2.5      Boiling  Tube  Thermocouples  Channels 
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C.    CONVECTION  SHIELD 

Initial  tests  on  the  smooth  copper  tube  revealed  an 
axial  thermal  gradient  in  the  liquid  over  the  length  of  the 
boiling  tube.  This  thermal  gradient  increased  with 
increasing  oil  concentration  and  with  increasing  heat  flux, 
with  the  thermal  gradient  ranging  from  0.5  K  to  1.5  K  over 
the  length  of  the  boiling  section.  This  thermal  gradient 
resulted  in  questionable  data  since  steady-state  conditions 
could  not  be  maintained  over  the  course  of  the  individual 
data  runs. 

During  the  early  stages  of  this  investigation,  it  was 
found  that  the  steady-state  conditions  were  fairly 
difficult  to  achieve.  This  was  especially  true  after  making 
a  step  change  in  the  heat  flux.  For  example,  if  the  heat 
flux  was  increased,  the  cooling  flow  rate  to  the  condenser 
must  be  increased  to  maintain  the  system  pressure.  Since 
the  system  showed  about  a  30-second  time  lag  between  the 
change  of  the  valve  position  and  the  subsequent  system 
response,  the  pressure  was  often  seen  to  undergo  cyclic 
variations.  With  considerable  experience  of  the  operator, 
it  was  possible  to  achieve  steady-state  conditions  within  5 
to  10  minutes,  as  demonstrated  by  Reilly  CRef.  53. 

In  order  to  minimize  the  efforts  necessary  by  the 
operator,  two  600-W  auxiliary  cartridge  heaters  were 
installed  within  the  boiling  tee.  With  this  arrangement, 
steady-state   conditions  could   be  maintained  very  easily  by 
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maintaining  the  total  heat  load  in  the  boiling  tee  at  a 
constant  value.  For  example,  if  the  heat  duty  for  the 
boiling  tube  was  increased,  the  heat  duty  on  the  auxiliary 
heaters  must  be  decreased  by  the  same  magnitude.  In  this 
manner,  steady-state  conditions  could  be  maintained  with 
minimum  attention.  However,  this  arrangement  required  that 
the  heat  duty  on  the  auxiliary  heaters  be  a  maximum  when 
the  heat  duty  on  the  boiling  tube  is  at  a  minimum.  As 
determined  by  the  early  data  runs  (see  Figure  2.6),  the 
convective  effects  created  by  the  bubbles  of  the  auxiliary 
heaters  created  an  artificially  enhanced  heat  -  transfer 
coefficient  at  low  heat-flux  settings. 

In  order  to  alleviate  the  problem  discussed 
above,  while  maintaining  the  advantages  offered  by 
the  auxiliary  heaters,  the  present  investigator 
considered  the  use  of  a  convection  shield.  As 
shown  in  Figure  2.2,  this  shield  was  positioned 
around  the  test  tube,  thus  isolating  it  from  the 
convective  bubble  effects  created  by  the  auxiliary  heaters. 

The  convection  shield  also  provided  increased 
mixing  of  the  refrigerant  and  condensate  without 
interfering  with  the  boiling  phenomenon  of  the  boiling 
tube  or  the  circulation  of  the  oi 1 -refrigerant 
mixture  in  the  boiling  section  by  both  the  boiling 
tube  and  the  auxiliary  heaters. 
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D.     DATA  ACQUISITION  AND  REDUCTION 

The  data-acquisition  system  used  for  this  research  was 
a  Hewlett-Packard  9826A  computer  and  a  Hewlett-Packard 
3497A  Automatic  Data  Acquisition  Unit.  The  HP-9826A 
computer  was  used  to  control  the  HP-3497A,  to  provide  hard- 
copy  storage  capacity  and  to  analyze  the  data  obtained 
during  this  investigation. 

The  data-acquisition  system  provided  the  means  by 
which  the  researcher  was  able  to  review  system  performance 
and  stability  prior  to  the  taking  of  any  data.  System 
parameters  were  monitored  on  the  video  screen  following  the 
scanning  of  all  of  the  channels  by  the  HP-3497A.  Following 
the  automatic  scanning  of  all  of  the  channels,  the  computer 
would  display  the  desired  saturation  temperature  (2.2  °C), 
the  observed  saturation  temperature,  the  liquid  temperature 
at  both  ends  of  the  boiling  section  and  their  average,  the 
vapor  temperature,  and  the  temperature  of  the  water- 
ethylene  glycol  mixture  sump  (-15  to  -18  °C).  Once  the 
steady-state  conditions  were  met  for  the  required  satura- 
tion temperature  within  the  prescribed  variance  (tO.l  K) 
of  the  desired  saturation  temperature,  the  data-acquisition 
system  was  initiated  for  the  individual  data  run. 
Since  there  was  a  time  lag  in  the  request  for  the 
taking  of  data  and  the  actual  printing  of  the  data  due  to 
the  sampling  of  each  channel  twenty  times  before  recording 
the   data,  the  data-acquisition  system  would  print  the  data, 

36 


al  lowing  the  researcher  to  review  the  results  prior  to 
storing  the  data.  If,  during  the  sampling  of  the  data,  a 
transient  trend  was  observed  due  to  the  cycling  of  the 
R-12  refrigerant  system  or  the  temperature  difference 
between  the  thermocouples  measuring  the  liquid  R-114  was 
not  within  the  specified  temperature  range,  the  researcher 
had  the  capability  to  reject  the  affected  data  set.  This 
capability  allowed  the  researcher  to  apply  a  cursory 
analysis  of  each  data  point  and  ascertain  if  the 
temperatures  measured  during  that  individual  data  set  were 
within  the  prescribed  limits  of  variance. 

The  information  scanned  on  each  channel  of  the  HP- 
3497A  were  thermal  emf's  provided  by  copper-cons tantan 
(type-T)  thermocouples  made  of  0. 25-mm-d iameter  (30  gauge) 
wire.  Two  additional  measurements  were  provided  to  the  HP- 
3497A  by  the  power -sens i ng  device,  as  described  by 
Karasabun  [Ref.  173,  which  converted  the  AC  voltage  dialed 
into  the  variac  of  the  main  boiling  tube  heater  into  a  DC 
signal,  allowing  the  HP-3497A  to  accurately  scan  the 
converted  value  and  then  allow  the  HP-9826A  to 
automatically  apply  this  value  to  the  analysis  calculation. 

Upon  the  initiation  of  the  command  to  take  data,  the 
data-acquisition  system  would  sample  each  channel  and  then 
compute,  according  to  the  step-wise  procedure  outlined  by 
Karasabun  [Ref.  17]  and  print  the  results.  Appendix  B 
contains   a  summary   of  the   data-acquisition  program   and  a 
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listing  of  the  data-reduction  program  utilized  during  the 
course  of  this  research.  Appendix  C  contains  a 
representative  sample  of  the  data  printout  received  for 
each  data  run. 
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III.  EXPERIMENTAL  PROCEDURES 

A.     INSTALLATION  OF  TUBE  IN  BOILING  APPARATUS 

The  external  surfaces  of  each  of  the  boiling  tubes 
were  cleaned  with  Nitol  (2  percent  nitric  acid  and  98 
percent  ethyl  alcohol)  solution  and  then  rinsed  with 
acetone  to  ensure  tube  cleanliness.  This  treatment  was  also 
effective  in  the  removal  of  oil  and  grease  residues  from 
the  tube  surfaces. 

After  air  drying,  the  individual  tube  to  be  tested  was 
installed  in  the  boiling  section  of  the  apparatus.  The 
system  was  sealed  and  then  evacuated  to  approximately  29 
inHg  by  means  of  a  portable  mechanical  vacuum  pump  (6), 
shown  in  Figure  2.1.  The  apparatus  remained  under  the 
evacuated  condition  for  a  period  of  about  one  hour  to  check 
if  there  were  any  major  leaks.  System  pressure  was  measured 
by  means  of  a  Marsh  pressure  gauge  (30  inHg  to  150  psi 
range  with  ±2.5  inHg  and  ±0.5  psi  accuracy, 
respectively).  If  there  was  not  a  noticeable  drop  in  the 
vacuum  reading  as  observed  on  the  pressure  gauge,  the 
system  was  charged  with  R-114  vapor,  which  was  condensed  by 
the  water-ethy 1 ene  glycol  mixture  and  then  the  absolute 
pressure  was  raised  to  approximately  0.186  MPa  (27  psi). 
Upon   reaching  this   pressure,   an   Automatic   Halogen   Leak 
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Detector,    TIF   5000,    with   a    sensitivity   of    3    ppm 
concentration,  was  used  to  check  for  system  leaks. 

Upon  successful  completion  of  the  leak-detection  tests 
and  taking  the  necessary  steps  to  fix  any  leaks  found,  the 
R-114  liquid  level  was  adjusted  to  be  20  mm  above  the  top 
of  the  boiling  tube  surface.  Liquid-level  adjustments  were 
made  by  either  boiling  off  the  excess  refrigerant  to  the 
R-114  reservoir  (3),  or  drawing  the  required  amount  of 
refrigerant  liquid  into  the  system  from  the  R-114 
reservoir.  The  required  liquid  level  for  boiling  was 
maintained  constant  for  all  tubes  tested  by  means  of  a 
permanent  scribe  mark  on  the  inside  of  the  oil  intake 
endbe 1 1  of  the  boiling  section.  At  this  initial  refrigerant 
level,  the  system  was  now  ready  for  boiling  tests  in  pure 
re  f  r  i  gerant . 

B.    GENERAL  SYSTEM  OPERATION 

Table  2  provides  a  summary  of  al 1  of  the  83  data  runs 
made  during  this  investigation.  Each  data  run  is  annotated 
to  indicate  the  conditions  involved  in  the  respective  run. 
The  data  runs  were  numbered  sequential ly  with  an  indicator 
code  for  each  tube  preceding  the  number  of  the  data  run. 
Since  this  research  is  a  follow  up  of  previous  research, 
the  initial  data  run  was  indicated  as  number  168.  Al  1  data 
runs   were  undertaken   at  a  saturation  temperature  of  2.2°C. 
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A   listing  of   the  boiling   tube  indicator  codes  for  Table  2 
is  as  foil ows : 

WH  --  Wieland  Hard  Copper  Smooth  Tube 

GTB  --  Wieland  Gewa-T  (1.02  fins  per  mm) 

TXE  --  Hitachi  Thermoexce 1 -E 

THE  --  Hitachi  Thermoexce 1 -HE. 

Initially,   all  data   runs  were   taken  over   a  range  of 

nine   different  heat-flux   settings,  ranging  from  2  kW/m2  to 

91   kW/m2 .    Additional  heat   fluxes  of   0.5,  0.8,  1.2  kW/m2 

were   later  included  for  all  tubes  with  enhanced  surfaces  to 

study   the  onset   of  nucleate  boiling  in  the  oil-refrigerant 

mixtures   and  to   study  the   effects  of   hysteresis  on  these 

tubes.    The   initiation  of  all  data  runs  was  carried  out  in 

the   same  manner  to  ensure  consistency  in  the  obtained  data. 

Each   run  was   started  by  an  initial  system  "warm-up"  period 

allowing   the  R-12   refrigeration  system  to  obtain  a  minimum 

coolant   starting  temperature   of  about   -15°C.     This   was 

accomplished    by   lighting    off   of   the   R-12   system   in 

conjunction   with  the   effluent   reci rcu 1  at i on/coo  1  i ng   pump 

(8)    and    slight    opening   of  the  control  valve  VC.   This 

procedure   allowed  the   system  to   gradually  decrease  system 

temperature   and   pressure,   while   minimizing   the   adverse 

effects   of   shocking   the   boiling   tube.   The   system   was 

operated   in  this   condition  until   a  stable   condition   was 

established   with  the  saturation  temperature  of  2.2°C.   This 

operation   was   most   critical   when   the   testing   sequence 

involved      increasing     heat-flux      conditions     since 
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TABLE  2 
Summary  of  Data  Runs 


Run 

% 

Heat 

No. of  Data 

Aux  i 1 iary 

Convect  i  on 

No. 

Oi  1 

F  1  ux 

Po  ints 

Heater 

Shield 

WH168 

0 

DEC 

20 

off 

wi  thout 

WH169 

0 

INC 

22 

off 

w  i  thout 

WH170 

1 

DEC 

20 

off 

w  i  thout 

WH171 

1 

INC 

22 

off 

w  i  thout 

WH172 

2 

DEC 

20 

off 

wi  thout 

WH173 

2 

INC 

20 

off 

w  i  thout 

WH174 

3 

DEC 

20 

off 

w  i  thout 

WH175 

3 

INC 

22 

off 

wi  thout 

WH176 

3 

DEC 

20 

off 

wi  thout 

WH177 

6 

DEC 

20 

off 

wi  thout 

WH178 

6 

INC 

21 

off 

w  i  thout 

WH179 

10 

DEC 

18 

off 

w  i  thout 

WH180 

10 

INC 

13 

off 

wi  thout 

WH181 

10 

INC 

12 

off 

wi  thout 

GTB182 

0 

DEC 

24 

on 

with 

GTB183 

0 

INC 

23 

on 

with 

GTB184 

0 

DEC 

22 

on 

with 

GTB185 

1 

DEC 

20 

on 

with 

GTB186 

1 

INC 

22 

on 

with 

GTB187 

2 

DEC 

34 

on 

with 

GTB188 

2 

INC 

18 

on 

-with 

GTB189 

3 

INC 

22 

on 

with 

GTB190 

3 

DEC 

20 

on 

with 

GTB191 

6 

DEC 

27 

on 

with 

GTB192 

6 

INC 

20 

on 

with 

GTB193 

10 

INC 

20 

on 

with 

GTB194 

10 

DEC 

21 

on 

with 

WH195 

0 

INC 

22 

on 

with 

WH196 

0 

DEC 

20 

on 

with 

WH197 

0 

DEC 

17 

on 

with 

WH198 

0 

INC 

22 

on 

with 

WH199 

1 

INC 

20 

on 

with 

WH200 

1 

DEC 

20 

on 

with 

WH201 

2 

INC 

20 

on 

with 

WH202 

2 

DEC 

20 

on 

with 

WH203 

3 

INC 

20 

on 

with 

WH204 

3 

INC 

20 

on 

with 

WH205 

3 

DEC 

20 

on 

with 

WH206 

3 

INC 

21 

on 

with 

WH207 

6 

DEC 

20 

on 

with 

UH208 

6 

INC 

20 

on 

with 

WH209 

10 

INC 

20 

on 

with 

42 


TABLE  2 
Summary  of  Data  Runs  (cont'd) 


Run 

% 

Heat 

No.  of 

No. 

Oi  1 

F  1  ux 

Poin 

WH210 

10 

DEC 

21 

TXE210 

0 

DEC 

18 

TXE211 

0 

DEC 

12 

TXE212 

0 

DEC 

12 

TXE213 

0 

DEC 

14 

TXE214 

0 

DEC 

12 

TXE215 

0 

INC 

27 

TXE216 

0 

DEC 

14 

TXE217 

1 

INC 

27 

TXE218 

1 

DEC 

14 

TXE219 

2 

DEC 

16 

TXE220 

2 

DEC 

16 

TXE221 

2 

INC 

25 

TXE222 

3 

DEC 

15 

TXE223 

3 

INC 

22 

TXE224 

3 

INC 

24 

TXE225 

6 

INC 

30 

TXE226 

6 

DEC 

16 

TXE227 

6 

INC 

25 

TXE228 

10 

INC 

24 

TXE229 

10 

DEC 

25 

TXE230 

10 

INC 

20 

THE231 

0 

INC 

20 

THE232 

0 

DEC 

12 

THE233 

1 

DEC 

14 

THE234 

1 

INC 

24 

THE235 

1 

INC 

27 

THE236 

2 

DEC 

14 

THE237 

2 

INC 

12 

THE238 

2 

INC 

8 

THE240 

3 

DEC 

14 

THE241 

3 

INC 

18 

THE242 

3 

INC 

8 

THE243 

3 

INC 

5 

THE244 

3 

INC 

5 

THE245 

6 

DEC 

14 

THE246 

6 

INC 

24 

THE247 

6 

INC 

18 

THE248 

10 

DEC 

14 

THE249 

10 

INC 

24 

:  i  1  i  ary 

Convect  i  on 

[eater 

Shield 

on 

with 

on 

with 

on 

with 

on 

with 

on 

with 

on 

with 

on 

with 

on 

with 

on 

with 

on 

with 

on 

with 

on 

with 

on 

with 

on 

with 

on 

w  i  th 

on 

with 

on 

with 

on 

with 

on 

w  i  th 

on 

with 

on 

with 

on 

with 

on 

with 

on 

with 

on 

with 

on 

with 

on 

with 

on 

with 

on 

with 

on 

with 

on 

with 

on 

with 

on 

with 

on 

with 

on 

with 

on 

with 

on 

with 

on 

with 

on 

with 

on 

with 
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failure  to  strictly  adhere  to  this  would  result  in 
premature  nucleate  boiling  and  the  subsequent  purging  of 
the  data.  This  initial  warm-up  period  required  a  minimum 
of  one  hour  to  ensure  adequate  system  performance. 

Upon  reaching  stable  conditions,  the  initial  heat  flux 
for  the  boiling  tube  was  set  by  means  of  the  resistor 
variac.  Also,  the  heat  flux  for  the  auxiliary  cartridge 
heaters  was  set  by  a  second  variac.  The  heat  flux  of  the 
auxiliary  cartridge  heaters  was  limited  to  settings  that 
allowed  the  heating  of  the  liquid  R-114  without  creating 
excess  turbulence  on  the  boiling  tube.  As  discussed 
earlier,  the  boiling  tube  was  shielded  from  the  effects  of 
the  auxiliary  heater  by  means  of  the  convection  shield.  As 
the  oil  concentration  of  the  r ef r i gerant-o i i  mixture 
increased,  the  heat-flux  setting  for  the  auxiliary  heater 
was  additionally  limited  to  lower  heat  fluxes.  This 
precluded  the  generated  oil  foam  from  carrying  over  the 
sides  of  the  convection  shield  and  interfering  with  the 
operation  of  the  boiling  tube. 

After  setting  of  the  desired  heat  flux  for  the  data 
run,  the  system  was  monitored  for  the  vapor  temperature, 
sump  temperature  and  the  liquid  refrigerant  temperature 
(two  readings;  left  side  and  right  side  of  the  boiling 
tube).  Control  valve  VC  was  adjusted  to  allow  the  system  to 
reach  steady-state  conditions.  Additionally,  the  liquid 
refrigerant   temperatures  were   compared  to   ensure   minimum 
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variance.  After  fulfilling  the  requirements  of  the  steady- 
state  conditions,  the  data-acquisition  system  was  initiated 
to  take  data. 

The  data-acquisition  system  would  scan  each  assigned 
channel,  compute  the  heat- transfer  coefficient  and  print 
the  results  for  operator  viewing.  If  the  data  were  within 
the  prescribed  limits,  they  were  stored  on  a  floppy  disk 
for  permanent  records.  A  representative  sample  of  a  data 
printout  is  contained  in  Appendix  B. 

After  each  run  was  completed,  oil  was  added  to 
the  boiling  section  by  means  of  a  graduated  cylinder 
through  valve  VI.  Table  3  contains  the  listing  of 
the  required  volume  of  oil  added  to  obtain  the 
proper  mass  percent  of  oil  for  each  series  of  data 
runs.  Mixing  of  the  oil  and  refrigerant  liquid  was 
accomplished  by  means  of  the  auxiliary  heater,  to 
ensure  a  homogeneous  mixture.  Foaming  of  the  oil- 
refrigerant  mixtures  was  observed  over  the  entire  heat- 
flux  range  of  each  data  run,  although  carryover  of 
the  foam  over  the  convection  shield  was  limited  to  the 
highest  heat  flux.  As  the  oil  concentration  increased, 
carryover  increased,  but  was  limited  to  the  vertical 
section  of  the  apparatus.  Carryover  oil  was  not 
observed  in  the  condenser  section  of  the  apparatus. 
Figures    3.1    to    3.3   demonstrate  the   foaming  action  of 
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the    oil-refrigerant    mixtures    at    a    heat   flux  of  37 

kW/m2   for  oil  concentrations  of  0,  6  and  10   percent. 

Data   runs  were   made  for  each  boiling  tube  tested  with 

pure    R-114   and   followed   by     runs   of   o i 1  - r ef r i gerant 

mixtures    of   increasing    oil   content.     Oi 1 -ref r i gerant 

mixtures    of   1,  2,   3,  6   and  10   percent  oil  by  mass  were 

used  for  all  tubes. 

TABLE  3 
OIL  CONCENTRATIONS  VALUES 

Total  Volume  (cc)   Step  Change  (cc) 
of  Oil  in  Volume  of  Oil 


26.8  26.8 

54.2  27.4 

82.2  28.0 

169.6  87.4 

295.2  125.6 


Percen 

t 

Oi  1 

0 

1 

2 

3 

6 

10 
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37  kW/m' 


91  kW/m 


Figure  3.1    Foaming  Action  for  Pure  R-114 
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Figure  3.2 


Foaming  Action  for  R-114  with 
6  Percent  Oil 
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37  kW/m' 


91  kW/m' 


Figure  3.3 


Foaming  Action  for  R-114  with 
10  Percent  Oil 
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IY.-_    RESULTS  AND  DISCUSSION 

A.    EFFECT  OF  CONVECTION  SHIELD 

Two  sets  of  data  runs  were  taken  on  the  smooth  copper 
tube;  one  set  without  the  convection  shield  to  validate  the 
performance  in  comparison  with  the  data  obtained  by  Reilly 
(at  a  saturation  temperature  of  -2.2°C)  [Ref.  5];  and  a 
second  data  run  using  the  convection  shield  and  the 
auxiliary  cartridge  heaters.  These  data  runs  allowed  for 
an  accurate  comparison  of  the  data  taken  over  a  period  of 
time  and  provided  justification  for  the  use  of  the 
convection  shield.  As  can  be  seen  from  Figures  4.1,  4.2 
and  4.3,  for  0,  3  and  10  percent  oil  concentration, 
respectively,  the  present  data  with  and  without  the  shield 
in  place  are  essentially  the  same.  However,  as  discussed 
in  section  II. C,  the  use  of  the  shield  provided  consistent 
steady-state  conditions  over  the  entire  range  of  heat 
fluxes  and  oil  concentrations,  thereby  minimizing  the 
efforts  required  by  the  operator. 

Figure  4.1,  which  contains  the  data  run  for 
pure  refrigerant  for  the  smooth  copper  tube  at  a 
saturation  temperature  of  2.2°C,  with  and  without  the 
shield,  demonstrates  the  repeatability  of  the  data. 
Also  plotted  on  this  figure  are  Reilly's  smooth-tube 
data    for    pure  refrigerant  at  a  saturation  temperature  of 
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-2.2°C.  This  data  set  represents  a  series  of  data 
runs,  starting  with  a  heat  flux  of  100  kW/m2  and 
decreasing  the  heat  flux  to  a  minimum  of  1  kW/m2 .  The 
three  sets  of  curves  are  nearly  identical  with  the 
deviation  between  the  Reil ly  data  and  the  most  recent  data 
being  attributed  to  the  pressure  effect  and  the  possible 
aging  of  the  tinned  interface  between  the  thermocouple 
sleeve  and  the  inner  surface  of  the  smooth  copper  tube. 
The  maximum  discrepancy  of  the  computed  wall  superheat 
between  Reilly's  data  and  the  present  data  is  about  twenty 
percent  at  the  maximum  heat-flux  setting.  This  figure  also 
shows  the  typical  uncertainty  of  about  10  percent  computed 
for  the  wall  superheat  values  during  this  investigation. 

B.     BOILING  PERFORMANCE  OF  THE  SMOOTH  TUBE 

Figure  4.4  shows  the  performance  of  the  smooth  copper 
tube  in  pure  R-114  and  refrigerant-oil  mixtures  of  3  and 
10  mass  percent  oil.  The  behavior  of  the  smooth  copper 
tube  in  pure  R-114  liquid  demonstrates  typical 
nucleate  pool-boiling  performance.  As  can  be  seen  in 
Figure  4.4,  the  line  between  points  A  and  B  represents 
a  region  of  constant  slope  (on  a  log-log  scale)  for 
the  variation  of  heat  flux  with  the  tube  wa 1 1  superheat 
(Two-Tsat).  This  section  of  the  curve  represents  natural 
convection   (i.e.,  no  bubbles  were  generated). 


54 


Upon  reaching  point  B  (see  Figure  4.4).  the  boiling 
tube  wall  superheat  begins  to  decrease  even  though  the 
heat  flux  is  continually  increased.  At  this  point,  the 
tube  begins  to  demonstrate  the  characteristics  of 
incipient  nucleate  boiling.  This  region  is  known  as  the 
mixed  boiling  region,  where  the  transition  from  natural 
convection  to  nucleate  boiling  occurs.  This  event  is 
characterized  by  the  increasing  number  of  active 
nucleation  sites  with  increasing  heat  flux.  It  should  be 
noted  that  during  al 1  data  runs,  the  unheated  ends  of  the 
tubes  did  not  show  any  nucleation  sites.  Reilly  CRef.  5] 
also  stated  that  the  unenhanced  ends  of  the  boiling  tube, 
at  all  observed  heat  fluxes,  underwent  only  natural 
convection  despite  a  small  amount  of  heat  axial ly  conducted 
along  the  tube  walls.  It  was  observed  that  this  transition 
from  natural  convection  to  nucleate  boiling  occurs  very 
quickly,  usually  within  a  few  seconds  after  the  activation 
of  the  first  nucleation  site. 

At  the  heat  flux  represented  by  point  C  in  Figure 
4.4,  all  the  available  activation  sites  of  the  boiling  tube 
appear  to  be  active  and  the  wa 1 1  superheat  begins  to 
increase  with  increasing  heat  flux.  From  point  C  to 
point  D,  where  the  heat  flux  is  increased  to  a  maximum 
of  91  kW/m2 ,  new  activation  sites  are  generated  (maximum 
heat   flux  was   limited   by   cartridge  heater  rated  output). 
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As    the   oil   concentration   was   increased,   there   was   a 
significant  increase  in  the  amount   of   foam    generated. 

After  obtaining  a  maximum  heat  flux  at  point  D,  the 
heat  flux  was  decreased  and  is  represented  as  the  region 
between  points  D  and  E  on  Figure  4.4.  This  region  of  the 
curve  represents  decreasing  wall  superheat;  however,  more 
activation  sites  that  were  created  at  the  higher  heat-flux 
settings  remain  active  to  lower  heat  fluxes.  This 
increased  range  over  which  the  boiling  tube  continues  to 
nucleate  results  in  an  increase  in  the  heat- transfer 
performance  since  the  wall  superheat  temperature  for  the 
nucleate  boiling  at  any  particular  heat  flux  is  less  than 
the  tube  wall  superheat  for  the  natural  convection  heating 
at  the  same  heat  flux. 

Stephan  CRef.  9]  reported  that  the  effects  of  adding 
oil  introduces  a  mass-diffusion  resistance  to  the  liquid 
and  therefore  lowers  the  heat- transfer  coefficient.  Reilly 
CRef.  5],  in  his  discussion  of  his  smooth-tube  data  at  a 
saturation  temperature  of  -2.2°C  (Figure  4.5),  stated  that 
the  3  mass  percent  oil  and  10  mass  percent  oil 
concentrations  demonstrated  higher  wall  superheat 
temperature  than  that  observed  in  the  pure  R-114  liquid. 
However,  Reilly  observed  that  the  wall  superheat  for  the  3 
percent  oil  concentration  was  higher  than  that  for  the 
10    percent   oil   concentration.   He  attributed  this  to  the 
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non-linear  physical  property  characteristics  of  the 
refrigerant-oil  mixtures.  As  can  be  seen  in  Figure  4.4, 
the  general  trend  of  the  data  is  the  same  with  the 
exception  of  the  3  percent  oil  concentration  wall 
superheat,  which  nearly  overlays  the  data  for  pure  R-114 
liquid.  This  deviation  from  the  Reilly  data  may  be  due  in 
part  to  the  difference  in  the  saturation  temperatures.  The 
addition  of  oil  to  the  pure  R-114  liquid,  up  to  a  maximum 
of  10  percent,  should  delay  the  transition  from  natural 
convection  to  nucleate  pool  boiling  on  the  tube.  With  a 
refrigerant-oil  mixture  of  10  mass  percent  oil,  it  was 
observed  that  the  surface  nucleation  sites  did  appear  to 
spread  more  slowly  as  the  heat  flux  was  increased.  This 
condition  continued  until  the  tube  was  fully  engulfed  in 
the  nucleate  boiling  mode.  This  observation  was  also  noted 
by  Reilly  [Ref.  5]  and  it  agrees  with  the  contention  of 
Chongrungreong  and  Sauer  [Ref.  8]  and  Thome  [Ref.  6]  that, 
with  the  exception  of  the  surface  tension,  the  non-linear 
variations  in  the  physical  properties  of  the  refrigerant- 
oil  mixtures  explain  the  variation  in  the  boiling  heat 
transfer  of  the  tube. 

Figure  4.6  shows  the  heat- transfer  coefficient  of  the 
smooth  tube  in  refrigerant-oil  mixtures,  plotted  as  a 
function  of  heat  flux.  Figure  4.7  shows  the  same  plot  for 
Reilly's    data   for   a   saturation   temperature   of   -2.2°C 
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[Ref.  5].  As  can  be  seen,  the  two  graphs  are  nearly 
identical,  representing  good  repeatability  of  the  data. 
Also,  the  effect  of  the  addition  of  up  to  6  percent  oil 
(Figure  4.7)  consistently  represents  an  approximate 
decrease  of  10  percent  in  the  heat- transfer  coefficient. 
As  the  oil  concentration  was  increased  from  6  to  10 
percent,  the  decrease  in  the  heat- transfer  coefficient  was 
a  maximum  of  about  35  percent.  These  decreases  in  the 
heat-transfer  coefficients  were  dependent  on  the  comparison 
of  the  pure  R-114  and  the  10  percent  oil  concentration 
at  equal   heat-flux  settings. 

C.     BOILING  PERFORMANCE  OF  THE  GEUA-T  TUBE 

Figure  4.8  shows  the  nucleate  pool-boiling  performance 
of  the  GEWA-T  (1.02  fins/mm)  tube  in  R-114-oil  mixtures. 
The  magnitude  of  the  tube-wall  superheat  was  considerably 
less  than  the  values  obtained  for  the  smooth  tube  discussed 
in  Section  IV. B.  This  lower  wall  superheat  can  be 
attributed  to  the  special ly  formed  fins  of  the  GEWA-T  tube. 
The  performance  of  the  GEWA-T  tube  is  similar  to  that  of 
the  smooth  tube  with  the  curves  of  the  increasing  heat  flux 
to  the  point  of  incipient  nucleate  boiling  (point  B), 
paralleling  the  performance  of  the  smooth  copper  tube 
(Figure  4.4).  The  same  condition  exists  throughout  the 
entire  range  of  decreasing  heat  fluxes. 
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The  region  bounded  by  points  A  and  B  shows  the 
characteristic  natural-convection  heating  behavior.  The 
onset  of  incipient  nucleate  boiling  at  point  B  is  clearly 
defined  for  the  pure  R-114  liquid.  However,  as  the  oil 
concentration  was  increased,  the  point  of  incipient 
nucleate  boiling  becomes  less  clearly  defined  and  the 
region  of  mixed  boiling  (point  B  to  C)  increased  with 
increased  heat  flux.  It  was  also  observed  that  in  the 
region  from  A  to  B,  the  3  percent  oil  mixture  showed  the 
smallest  wall  superheat  and  the  pure  refrigerant  showed  the 
highest  wall  superheat.  The  reason  for  this  behavior  is 
not  fully  known  at  present,  but  may  be  due  to  convection 
effects  within  the  channels  of  the  GEUA-T  tube  which  are 
enhanced  by  the  foaming  action  in  the  presence  of  oil. 

The  region  bounded  by  the  points  C  and  D  for  the 
increasing  heat-flux  condition  is  indicative  of  fully 
developed  nucleate  pool  boiling.  The  region  between  points 
D  and  E  represents  nucleate  pool  boiling  for  decreasing 
heat  flux.  Figure  4.8  shows  that,  as  the  oil 
concentration  was  increased,  the  wall  superheat  increased 
at  a  given  heat  flux.  This  is  contrary  to  the  results  for 
the  condition  of  increasing  heat  flux  in  the  region  of 
natural-convection  heating. 

When  comparing  the  performance  of  the  smooth  tube 
(Figure   4.4)  to   that  of   the  GEWA-Tube   (Figure  4.8),   the 
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maximum  adverse  effect  of  the  increased  oil  content  for  the 
smooth  tube,  occurs  at  the  higher  heat-flux  settings.  For 
the  GEUA-T  tube,  the  maximum  decrease  in  the  tube 
performance  for  increased  oil  content,  occurs  in  the  lower 
and  intermediate  heat-flux  range. 

Figure  4.9  shows  the  heat- transfer  coefficient  as  a 
function  of  heat  flux  over  the  range  of  R-114-oil 
mixtures  (0,  1,  2,  3,  6  and  10  percent  oil).  Data  observed 
over  the  range  of  R-114-oil  mixtures  of  1,  2,  3  and  6 
percent  oil  are  nearly  coincident.  The  10  percent  oil 
concentration  data  represent  an  approximate  30  percent 
reduction  of  the  heat- transfer  coefficient  when  compared  to 
the  heat- transfer  coefficient  obtained  during  the  data  for 
pure  R-114.  Figure  4.10  shows  the  direct  comparison  of  the 
GEWA-T  tube  with  the  smooth  tube  for  R-114-oil  mixtures  of 
0,  3  and  10  percent  oil.  The  GEWA-T  tube  represents  an 
overall  enhancement  of  4.2  over  the  baseline  smooth  copper 
tube  at  a  heat  flux  of  40  kW/m2  with  pure  R-114  and  an 
overall  enhancement  of  3.8  with  10  percent  oil.  These 
results  compare  favorably  with  the  results  of  Yilmaz, 
Hwalek  and  Westwater  [Ref.  18]  for  a  12.3-mm  OD  GEWA-T  tube 
in  p-xylene,  in  which  they  obtained  an  enhancement  of  5.3 
in  the  pool-boiling  heat- transfer  coefficient  over  their 
baseline  smooth  copper  tube.  Additionally,  Yilmaz  and 
Westwater     [Ref.  12]  tested  a  similar  surface  in  isopropyl 
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alcohol  in  which  an  enhancement  of  2.0  was  observed.  Marto 
and  Lepere  CRef.  13]  tested  a  17.9-mm  OD  GEWA-T  surface  in 
both  R-113  and  FC-72  and  reported  enhancements  of  2.8  and 
2.5,  respectively. 

Figure  4.11  shows  more  clearly  the  degradation 
that  oil  causes  in  the  boiling  heat-transfer 
performance  for  the  GEWA-T  tube.  This  figure  plots 
the  heat- transfer  coefficient  of  the  GEWA-T  tube  in 
the  oil  mixture  relative  to  the  heat- transfer 
coefficient  in  pure  R-114  as  a  function  of  oil 
concentration.  From  the  shape  of  the  curves,  the 
oil-caused  degradation  does  not  show  a  clear  trend 
with  the  oil  concentration  at  all  practical  heat 
fluxes.  Over  the  range  of  oil  concentrations,  except 
at  a  heat  flux  setting  of  37  kW/m2 ,  there  is  only 
a  smal 1  degradation  (up  to  a  maximum  of  35%)  in  the 
boiling  performance.  This  is  probably  due  to  the 
large  channels  and  reentrant  cavities  in  the  tube, 
allowing  the  "pumping"  action  of  the  bubbles  to  continue 
to  remove  the  oil-rich  liquid  throughout  the  entire 
heat-flux  range  without  a  significant  reduction  in 
performance.  The  only  significant  decrease  in 
performance  occurred  at  heat  fluxes  greater  than  5 
kW/m2   with   10   percent   oil. 
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D.     BOILING  PERFORMANCE  OF  THE  THERMOEXCEL-E  TUBE 

Figure  4.12  shows  the  nucleate  poo  1 -bo i 1 i ng 
performance  of  the  Thermoexce 1 -E  tube  in  R-114-oil 
mixtures.  As  can  be  seen  in  the  figure,  the  region  of 
natural-convection  heating  for  the  condition  of  increasing 
heat  flux,  points  A  to  B,  for  the  pure  R-114  is  very  small, 
with  the  point  of  incipient  nucleate  boiling  occurring  at  a 
reasonably  low  heat  flux  (1.5  kW/m2 ) .  At  this  point  of 
incipient-  nucleate  boiling,  the  tube  performs  in  the  mixed 
boiling  region  as  a  limited  number  of  nucleation  sites 
become  active.  From  point  C  to  C,  the  tube  appears  to 
operate  in  more  of  a  natural-convection  heating  mode  as  no 
new  nucleation  sites  are  activated,  although  previously 
activated  sites  continue  to  nucleate.  This  phenomenon  is 
probably  due  to  the  inability  of  the  cartridge  heaters  to 
produce  uniform  heat  flux  at  these  lower  settings.  At 
point  C  with  a  heat-flux  setting  of  5  kW/m2 .  more 
nucleation  sites  are  generated  and  as  the  heat  flux  is 
increased  to  60  kU/m2 ,  the  tube  becomes  fully  nucleated. 

For  the  3  percent  oil  concentration,  the  point  of 
incipient  nucleate  boiling,  point  B',  decreases  as  the  non- 
linear physical  properties  of  the  R-114-oil  mixture  allow 
the  onset  of  the  mixed  boiling  region  to  occur.  Also,  the 
region    of   the    mixed   boiling,    point   B*    to   D',   is 
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far  larger  than  that  for  pure  R-114,  as  the  generation  of 
new  nucleation  sites  is  retarded  by  the  increased  oil 
content . 

For  the  10  percent  oil  concentration,  on  the  other 
hand,  the  point  of  incipient  nucleate  boiling,  point  B'*, 
is  delayed  considerably.  At  approximately  14  kW/m2 ,  the 
transition  from  natura 1 -convect i ve  heating  to  nucleate 
boiling  occurs.  The  exact  reason  for  the  difference  between 
the  3  percent  and  10  percent  behavior  is  not  known,  however 
it  is  possibly  due  to  the  change  in  the  physical  properties 
of  the  R-114-oil  mixture. 

The  region  of  nucleate  pool  boiling  for  the  condition 
of  decreasing  heat  flux  represented  by  points  D  to  E,  D' 
to  E'  and  D,f  to  Ef '  for  0,  3  and  10  percent  oil 
concentrations,  respectively,  for  the  Thermoexce 1 -E  tube, 
represent  classic  nucleate  pool-boiling  characteristics. 
As  the  oil  concentration  increases,  the  wall  superheat  also 
increases  throughout  the  range  of  heat-flux  settings.  At 
the  lower  heat  fluxes,  the  largest  increase  in  the  wall 
superheat  occurred  at  3  percent  oil  concentration  when 
compared  to  the  pure  R-114,  while  the  10  percent  oil 
concentration  was  observed  to  undergo  continual  degradation 
in  tube  performance  over  the  entire  range  of  heat-flux 
settings,  although  a  more  noticeable  reduction  occurred  at 
the    higher   heat  fluxes.    At   a  heat-flux   setting  of   91 
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kW/m2 ,  the  R-114-oil  mixture  of  10  percent  oil  was  observed 
to  undergo  a  reduction  in  tube  performance  by  a  factor  of 
3.  For  the  practical  heat-flux  range  of  30  to  40  kW/m2 , 
there  is  an  increase  in  wa 1 1  superheat  by  factors  of  1.4 
and  2.6  for  R-114-oil  mixtures  of  3  and  10  percent, 
respect i  ve 1 y . 

Figure  4.13  shows  the  heat-transfer  coefficient  of  the 
Thermoexce 1 -E  tube  in  r ef r i ger ant -o i 1  mixtures,  plotted  as 
a  function  of  heat  flux.  As  can  be  seen  in  the  figure,  the 
heat- transfer  coefficient  remains  nearly  constant 
throughout  the  range  of  heat-flux  settings  for  pure  R-114, 
while  there  are  slight  reductions  for  the  1,  2  and  3 
percent  oil  concentrations.  Additionally,  it  can  be  seen 
that  at  these  lower  oil  concentrations,  the  heat- trans fer 
coefficients  are  nearly  coincident.  As  the  R-114-oil 
mixture  oil  concentration  is  increased  to  6  percent,  the 
tube  performance  increases  with  increasing  heat  flux  until 
reaching  a  heat-flux  setting  of  60  kW/m2 ,  at  which  time 
tube  performance  begins  to  deteriorate.  This  is  contrary 
to  the  tube  performance  for  0,  2,  and  3  percent  oil 
concentrations  in  which  the  heat- transfer  coefficient  was 
observed  to  increase  slightly  over  the  range  of  heat 
fluxes.  For  the  10  percent  oil  concentration,  the  heat- 
transfer  coefficient  decreased  as  the  heat  flux  was 
increased . 
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Figure  4.14  shows  a  direct  comparison  of  the 
Thermoexce 1 -E  tube  to  the  smooth  copper  tube  for  R-114-oil 
mixtures  of  0,  3  and  10  percent.  The  Thermoexce 1 -E  tube 
presents  enhancements  of  6.6  and  4.1  over  the  performance 
of  the  smooth  copper  tube  in  pure  R-114  and  an  R-114-oil 
mixture  of  10  percent  oil,  respectively,  at  a  heat-flux 
setting  of  40  kW/m2 .  The  results  for  the  pure  R-114  compare 
favorably  with  the  results  obtained  by  Yilmaz,  Palen  and 
Taborek  [Ref.  19]  for  the  Thermoexce 1 -E  in  p-xylene. 

Figure  4.15  shows  the  degradation  that  oil  causes  in 
the  boiling  heat- transfer  performance  for  the  Thermoexce 1 -E 
tube.  This  figure  plots  the  ratio  of  the  heat  -  trans fer 
coefficient  of  the  Thermoexce 1 -E  tube  with  oil  to  the  value 
without  oil  as  a  function  of  oil  concentration.  The  shapes 
of  the  curves  continue  to  demonstrate  the  non-linear 
variation  of  the  physical  properties  of  refrigerant-oil 
mixtures.  However,  the  significant  degradation  of  the  tube 
performance  is  probably  due  to  the  inability  of  the 
generated  bubbles  to  scavenge  the  interior  tunnels  of  the 
Thermoexce 1 -E  tube.  This  is  in  contrast  to  the  GEWA-T 
tube,  where  more  scavenging  is  possible  due  to  its  larger 
channel  size. 

E.     BOILING  PERFORMANCE  OF  THE  THERMOEXCEL-HE  TUBE 

Figure  4.16  shows  the  nucleate  pool-boiling 
performance    of   the    Thermoexce 1 -HE   tube    in   R-114-oil 
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mixtures.  In  the  region  between  points  A  and  B  /  B  '  , 
natural-convection  heating  occurs  for  increasing  heat  flux. 
Although  similar  in  performance  to  the  Thermoexce 1 -E  tube, 
the  Thermoexce 1 -HE  tube  presented  a  more  linear 
characteristic  than  that  of  the  Thermoexce 1 -E  tube, 
representing  a  more  uniform  heat  distribution.  At  point  B 
on  the  figure,  the  point  of  incipient  nucleate  boiling  for 
pure  R-114  liquid  occurs  at  approximately  2  kW/m2 .  For 
R-114-oil  mixture  of  3  and  10  percent  oil  concentrations, 
the  point  of  incipient  nucleate  boiling,  point  B',  occurs 
at  approximately  8  kW/m2 .  Also  in  this  region,  the  pure 
R-114  and  the  R-114-oil  mixture  of  3  percent  oil 
concentration  are  nearly  coincident  in  the  range  of  heat 
flux  from  500  W/m2  to  2  kW/m2. 

From  points  B  to  C,  the  Thermoexce 1 -HE  tube  in  the 
pure  R-114  liquid  achieved  full  nucleate  boiling  at 
approximately  37  kW/m2 ,  with  intermediate  data  points  in 
the  region  of  mixed  boiling.  These  intermediate  data 
points  indicate  that  the  Thermoexce 1 -HE  tube  was  not  able 
to  activate  all  nucleation  sites  simultaneously  for 
increasing  heat-flux  conditions.  The  R-114-oil  mixture  of 
10  percent  oil  was  able  to  achieve  activation  of  all 
nucleation  sites  upon  reaching  22  kW/m2 ,  without  any  data 
points  contained  in  the  mixed  boiling  region.  With  3 
percent   oil,  the   Thermoexce 1 -HE  tube   was  able   to  achieve 
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full  nucleate  boiling  at  37  kW/m2 .  The  Thermoexce 1 -HE 
tube  in  the  R-114-oil  mixtures  appeared  to  activate 
all  nucleation  sites  simultaneously. 

The  region  from  points  C  to  D  and  from  C '  to  D* 
represent  the  nucleate-  pool  boiling  for  increasing  heat 
flux.  The  performance  of  the  Thermoexce 1 -HE  tube  in  this 
region  appears  to  decrease  as  the  oil  content  was 
increased.  The  10  percent  oil  concentration  increased  the 
wall  superheat  by  a  factor  of  approximately  2  above  the 
wall  superheat  for  the  pure  R-114  liquid  for  heat-flux 
settings  of  40  kW/m2 . 

The  regions  between  the  points  D  and  E  and  D'  and  E' 
represent  nucleate  pool-boiling  for  decreasing  heat-flux 
conditions.  As  can  be  seen  from  the  figure,  the  greatest 
reduction  in  performance  occurs  in  the  lower  heat-flux 
settings  between  the  pure  R-114  and  the  3  percent  oil 
concentration.  This  is  similar  to  the  results  for  the 
Thermoexce 1 -E  tube,  where  the  tube  performance  decreased 
as  oil  content  increased.  There  is  a  significant  decrease 
in  the  Thermoexce 1 -HE  tube  performance  over  the  range  of  - 
heat-flux  settings  as  the  oil  concentration  is  increased  to 
10  percent. 

Figure  4.17  shows  the  heat- transfer  coefficient  of  the 
Thermoexce 1 -HE  tube  in  refrigerant-oil  mixtures,  plotted 
as  a  function  of  heat  flux.  As  can  be  seen,  the 
performance   is  very   similar  to   that  of   the  Ther moe xce 1 -E 
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tube  (Figure  4.13).  As  was  observed  for  the  Thermoexce 1 -E 
tube,  the  Thermoexce 1 -HE  tube  undergoes  a  slight  reduction 
in  the  heat- transfer  coefficient  as  the  oil  concentration 
was  increased  from  1  to  3  percent.  As  the  oil 
concentration  is  increased  to  6  and  10  percent  oil,  the 
performance  reduction  is  quite  noticeable.  With  the  10 
percent  oil  concentration,  tube  performance  decreases  at 
both  the  lower  and  higher  heat-flux  settings.  This  is 
contrary  to  the  findings  for  the  Thermoexce 1 -E  tube  where 
the  performance  decreased  over  the  entire  range  of  heat- 
flux  settings,  from  low  heat  flux  to  high  heat-flux 
sett i  ngs . 

Figure  4.18  shows  a  comparison  of  the  Thermoexce 1 -HE 
tube  to  the  smooth  copper  tube  at  R-114-oil  mixtures  of  0, 
3  and  10  percent  oil  concentrations.  The  Thermoexce 1 -HE 
tube  displayed  an  enhancement  of  6.1,  4.6  and  4.1  at  a  heat 
flux  setting  of  40  kW/m2  and  R-114-oil  mixtures  of  0.  3  and 
10  percent  oil,  respectively,  over  the  same  parameters  for 
the  smooth  copper  tube.  These  results  are  nearly  identical 
to  those  found  for  the  Thermoexce 1 -E  tube. 

Figure  4.19  shows  the  degradation  that  the  oil  has  on 
the  boiling  heat- transfer  performance  of  the  Thermoexce 1 -HE 
tube.  At  the  lower  heat  fluxes,  the  effect  of  the  oil 
concentration  is  slightly  more  pronounced  than  that 
observed   for  the  Thermoexce 1 -E  tube  throughout  the  range  of 
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oil  concentrations,  although  more  severe  at  the  lower  oil 
concentrations.  At  the  higher  heat  fluxes,  the  trend 
appears  to  reverse. 

F.     COMPARISON  TO  REILLY'S  DATA  FOR  THE  HIGH  FLUX  TUBE 

Figures  4.20,  4.21  and  4.22  provide  a  comparison  of 
the  smooth,  High  Flux,  GEWA-T,  Thermoexce 1 -E  and 
Thermoexce 1 -HE  tubes  for  0,  3  and  10  percent  oil 
concentrations,  respectively.  For  pure  R-114  liquid,  the 
High  F 1 ux , Thermoexce 1 -E  and  Thermoexce 1 -HE  tubes  operating 
at  a  heat  flux  of  approximately  20  kW/m2 ,  show  nearly 
identical  performance.  As  the  heat  flux  was  increased,  the 
performance  of  the  High  Flux  tube  continued  to 
increase  while  the  performance  of  Thermoexce 1 -E  and 
Thermoexce 1 -HE  tubes  remained  nearly  constant.  At  heat- 
flux  settings  lower  than  20  kW/m2 ,  the  Thermoexce 1 -E  and 
Thermoexce 1 -HE  tubes  showed  higher  heat  -  transfer 
coefficients  than  that  for  the  High  Flux  tube.  As  the  oil 
concentration  was  increased  to  3  percent,  the  performance 
of  all  five  tubes  decreased,  as  can  be  seen  in  Figure  4.21. 
Additionally,  the  GEWA-T  tube's  performance  appeared  to 
decrease  the  least  since,  as  discussed  earlier  in  section 
IV. C,  the  finned  surface  of  the  GEWA-T  tube  is  not 
adversely  effected  by  the  small  increase  in  the  oil  content 
due  to  its  large  reentrant  cavities.  From  this 
observation,   it  can   be  stated   that  the  performance  of  the 
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tubes  with  porous  coatings  or  small  reentrant  cavities  is 
adversely  affected  by  the  increase  in  oil  content  due  to 
the  generation  of  an  oil-rich  film  in  the  nucleation  sites 
impeding  the  tube's  performance.  This  is  further  evidenced 
by  Figure  4.22.  which  depicts  the  performance  of  all  five 
tubes  in  R-114-oil  mixtures  of  10  percent  oil.  The 
performance  of  all  the  tubes  is  reduced  significantly, 
while  the  performance  of  the  High  Flux  tube  is  severely 
decreased  in  the  range  of  heat-flux  settings  from  22  kW/m2 
to  the  maximum  setting  of  91  kW/m2 .  At  approximately  45 
kW/m2 ,  the  performance  of  the  High  Flux  tube,  the  GEWA-T, 
the  Thermoexce 1 -E  and  the  Thermoexce 1 -HE  tubes  are 
approximately  equal,  although  the  performance  of  the 
Thermoexce 1 -E  and  the  Thermoexce 1 -HE  tubes  remain  nearly 
constant  over  the  range  of  heat-flux  settings  of  14  kW/m2 
to  50  kW/m2 .  Although  the  performance  of  the  High  Flux 
tube  exceeds  that  of  the  GEWA-T,  the  Thermoexce 1 -E  and  the 
Thermoexce 1 -HE  tubes  over  this  range,  its  performance 
begins  to  severely  degrade  at  the  midpoint  of  the  heat-flux 
range  indicated.  In  the  heat-flux  range  above  45  kW/m2 , 
the  GEWA-T  tube  outperformed  all  the  other  tubes.  This  is 
attributed  to  the  ability  of  bubbles  generated  in  its 
larger  channels  to  continually  scavenge  the  oil-rich  layer 
from  the  tube,  and  therefore  remove  the  insulating  oil  from 
the  tube's  channels  and  exterior  surface. 


89 


V.   CONCLUSIONS  AND  RECOMMENDATIONS 


A.    CONCLUSIONS 

1.  In  pure  R-114,  at  a  heat  flux  of  30  kW/m2  the 
pool-boiling  heat- transfer  coefficient  of  the 
Thermoexce 1 -E,  Thermoexce 1 -HE  and  GEWA-T  tubes 
were  approximately  6.6,  6.1  and  4.2  times, 
respectively,    larger  than  the  baseline   smooth  tube. 

2.  In  R-114-oil  mixtures  of  3  percent  oil,  at  a 
heat  flux  of  30  kW/m2  ,  the  poo.l -bo  i  1  i  ng  heat- 
transfer  coefficient  of  the  Thermoexce 1 -E,  Thermo- 
excel-HE  and  GEWA-T  tubes  were  approximately  4.8, 
4.6  and  4.0  times,  respectively,  larger  than  the 
smooth   tube   under  the   same  conditions. 

3.  In  R-114-oil  mixtures  of  10  percent  oil,  at  a  heat 
flux  of  30  kW/m2 ,  the  performance  of  the  Thermoexcel- 
E  and  Thermoexce 1 -HE  tubes  were  nearly  identical, 
with  an  enhancement  of  4.1  over  the  pool-boiling 
heat- transfer  coefficient  of  the  smooth  tube.  The 
GEWA-T  tube,  under  the  same  conditions,  produced  an 
enhancement  of  3.8  over  the  smooth  tube. 

4.  As  oil  content  in  the  refrigerant-oil  mixture 
increased,  the  performance  of  enhanced  surfaces 
decreased;  however,  the  amount  of  this  decrease  was 
directly  related  to  the  size  of  the  subsurface 
channels  (i.e.,  the  larger  the  reentrant  cavity  or 
channel,  the  less  adverse  effect  of  the  increased  oil 
content ) . 

5.  As  the  oil  content  increases  in  the  channels,  the 
greater  the  adverse  effect  of  the  oil-rich  layer  has 
on  the  surface,  insulating  the  surface  from  the 
refrigerant  and  therefore  decreasing  the  performance 
of  the  surface. 

6.  The  GEWA-T  surface  in  refrigerant-oil  mixtures,  with 
its  larger  subsurface  channels,  was  able  to 
scavenge  the  oil-rich  layer  from  its  cavities  by 
the  bubbling  action  generated  within  the  channels. 

7.  The  performance  of  the  Thermoexce 1 -E  and  Thermoexce 1 - 
HE  tubes,  with  their  smaller  pore  size  and  smaller 
tunnels,  showed  significant  performance  decreases 
as  oi  1  content  was  increased.  The  smal  ler  pore  and 
tunnel  size  under   these   conditions,   were  unable  to 
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effectively  disperse  the  oil-rich  film  from  the  sub 
-surface  cavities  and  therefore  allowing  the  oil-rich 
film  to  insulate  the  tube's  surface. 


RECOMMENDATIONS 

1.  The  single-tube  testing  should  be  expanded  to  include 
multi-tube  bundles.  This  would  provide  more  realistic 
data  with  regard  to  the  effects  of  tubes  operating  in 
close  proximity  to  each  other. 

2.  A  uniform  heat  flux  is  necessary  within  the  boiling 
tube  and  within  the  boiling  section  itself.  There- 
fore, the  present  cartridge  heaters  need  to  be 
replaced  with  more  precise  heaters,  capable  of  pro- 
ducing a  uniform  heat  flux  in  the  axial  direction. 

3.  The  testing  of  surfaces  should  be  expanded  to  study 
the  effects  of  various  refrigerant  oils. 

4.  Conduct  experiments  to  establish  the  physical  proper- 
ties of  refrigerant-oil  mixtures  to  enable  further 
understanding  into  the  boiling  mechanisms  involved. 

5.  There  needs  to  be  additional  studies  undertaken  on 
the  effects  of  additional  boiling  temperatures  on  the 
mixture  properties  of  the  refrigerant-oil  mixtures. 
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APPENDIX  A 

THERMOCOUPLE  CALIBRATION   CRef.  53 

Karasabun  [Ref.  17]  described  the  thermocouple 
calibration  procedure  in  great  detail.  Reilly  [Ref.  51 
stated  that  the  data-reduction  program  utilized  differences 
in  thermocouple  readings  in  al 1  computations,  such  as  the 
amount  of  the  tube  wall  superheat  (Two  -  Tsat).  Reilly 
stated  that  the  calibration  of  the  thermocouples  would  be 
necessary  if  the  thermocouples  were  measuring  quantities 
such  as  absolute  temperature  to  determine  actual  fluid 
properties.  The  procedure  utilized  by  Karasabun  and  Reilly 
for  the  initial  instrumentation  for  the  experimental 
apparatus  is  as  contained  in  Appendix  A  of  Reference  5. 

To  determine  the  calibration  for  the  thermocouples, 
two  thermocouples  were  made,  one  from  the  beginning  of  the 
wire  rol 1  and  the  second  from  the  end.  The  seventh-order 
polynomial  equation  provided  by  the  wire  manufacturer  was 
corrected  by  adding  to  it  a  second-order  curve  fit  of  the 
variation  of  the  manufacturer's  computed  temperature  for  a 
given  emf  from  a  known  set  of  reference  temperatures  as 
measured  by  a  Hewlett-Packard  2804A  quartz  thermometer.  The 
thermometer  had  a  temperature  resolution  of  to. 0001  K  and 
an  accuracy  of   ±0.03  K. 
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The  manufacturer's  emf  to  temperature  conversion 
equation  is: 

T  =  aQ  +  a,  E  +  a,  E2  ♦  a,  E3  +  a*  E4  +  (A.  1) 

a5  Es  +  a6  E6  +  a7  E7 
where 

E   =   thermocouple  reading  in  volts 

T   =   temperature 

a0  =   0.100860910 

a,  =   25727.94369 

a,  =  -767345.8295 

a3  =   78025595.81 

a.  =  -9247486589 

as     =       6.97688X101  l 

a6     =     -2.66192xlOl 3 

a7     =       3. 94078x10* * 

Figure  A.l  CRef.  5]  shows  the  difference  in  the  quartz 
thermocouple  readings  and  the  actual  thermocouple  readings 
(discrepancy)  plotted  as  a  function  of  temperature.  As  can 
be  seen  from  the  figure,  the  two  thermocouples  agreed  to 
within  0.05  K  of  each  other.  The  manufacturer's 
calibration  equation  needed  an  increase  of  approximately 
0.1  K  to  more  accurately  convert  the  emf's  to  the  true 
temperature.  The  second-order  polynomial  correction  was  as 
foil ows ; 

DCP  =  b0  +  b, T  +  b2T:  (A. 2) 

where 
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DCP  =   discrepancy  (K) 

b0  =   8. 6268968x10- 2 

b,  =   3. 7619902x10- 3 

b?  =  -5.0689259x10-= 

T2  =  thermocouple  readings  from 
equa t  i  on  A .  1  in  C 
Therefore,  the  temperatures  computed  by  the  data- r educt i on 
program  (DRP5)  were  emf's  converted  directly  to 
temperatures  by  equation  A.l  and  corrected  by  equation  A. 2. 
This  calibration  was,  thus  estimated  to  result  in 
temperature  readings  within  ±0.05  K  of  the  true 
temperature . 
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APPENDIX  B 
DATA-REDUCTION  PROGRAM  CREF.  5] 

The    data-reduction   program   below   consists   of   the 
following  sections: 

Main  Program  -  Menu  of  subprogram  options 
Sub  Main      -  Takes  or  reprocesses  data 
Sub  plot      -  Plots  data  on  log-log  scale 
Sub  Poly      -  Computes  least  squares  curve  fit 
Sub  Plin      -  Plots  data  on  1 i near  -  1 i near  scale 
Sub  Stats     -  Computes  average  and  standard 

deviation  of  data 

Subprogram  Main  consists  of  the  following  steps: 

1.  Creates  data  file  for  data   and   plotting   file. 

2.  Annotates  tube  type. 

3.  Monitors  heat  flux  or  saturation  temperature  to 
establish  steady-state  conditions. 

4.  Scans  HP3497  channels  and  stores  in   data   file. 

5.  Converts   raw  emf's  to  temperatures,  current    and 
voltages. 

6.  Computes  the  heat- transfer  rate  for  the  tube's 
internal  cartridge  heater. 

7.  Computes  the  average  wall  temperature  of  the 
tube,  the  average  wall  superheat  and 
the    film      temperature. 
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8.  Computes    the    physical    properties    of   the 
R-114  using  the  given  correlation  temperature. 

9.  Computes   the   natural-convection   heat- transfer 
coefficient     of    R-114     for   the   non-boiling 
tube  ends. 

10.  Computes     the    heat   loss   from  the  non-boiling 
tube  ends. 

11.  Calculates    the    corrected   heat   fluxes   from 
the   testing  tube  to  the  R-114. 

12.  Calculates  the  boiling   heat- transfer  coefficient 
of  the  R-114  from  the  tube  being  tested. 

13.  Prints  the  data  and  stores  the  heat  flux  and  tube 
wall  superheat  values  in  the  plot  file. 

The  following  is  the  listing  of  the  complete  data-reduction 
program  (DRP5)  written  in  Basic  3.0  for  the  Hewlett-Packard 
9826  computer. 
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1000!  FILE  NAME:  DRPS 

1003!  DATE:      Octooer  19.  1934 

1006!  REVISED:   April  25.  1986 

1009! 

1012  COM  /Idp/  IdP 

1015  PRINTER  IS  1 

1018  CALL  Select 

1021  INPUT  "WANT  TO  SELECT  ANOTHER  OPTION  ( 1 -Y. 0-N)?" . Isel 

1024  IF  Isel-1  THEN  GOTO  1018 

1027  BEEP 

1030  BEEP 

1033  PRINTER  IS  1 

1036  PRINT  "DATA  COLLECTION/REPROCESSING  COMPLETED" 

1039  END 

1042  SUB  Main 

1045  COM  /Idp/  Idp 

1048  COM  /Cc/  C(7). leal 

1051  COM  /Wil/  D2.Di,Do.L.Lu.Kcu 

1054  DIM  Emt" (12)  ,T(12),D1a(  11  ).D2a(  11  )  .Dia(  11  ).Doa(11  ).La(11)  ,Lua(11 ) ,Kcua( 1 1  ) 

Et( 19) ,Tn$(4)[15] 

1057  DATA  0.10086091. 25727. 94369. -767345. 3295. 78025595. 31 

1060  DATA  -9247486589. 6. 97688E+11 , -2. S6192E+13, 3. 94078EH4 

1063  READ  CO) 

1066!  DATA  "Smooth" , "High  Flux"  ,  "Thermoexel-E" , "Thermoexel-HE" 

1069  DATA  Smooth. Hiah  Flux , Turbo-B.Hiah  Flux  Mod,Turbo-B  Mod 

1072  READ  TnS<+> 

1075  PRINTER  IS  701 

1078  BEEP 

1081  IF  Idp-4  THEN  PRINTER  IS  1 

1084  IF  Idp-4  THEN  GOTO  1774 

1087  INPUT  "ENTER  MONTH.  DATE  AND  TIME  (MM:DD:HH:MM:SS) " .DateS 

1090  OUTPUT  709:"TD":DateS 

1093  OUTPUT  709; "TD" 

1096  ENTER  709:DateS 

1099  PRINT 

1102  PRINT  "  Month,  date  and  time  :";Date$ 

1105  PRINT 

1108  PRINT  USING  "10X,""MOTE:  Program  name  :  DRP4""" 

1111  BEEP 

1114  INPUT  "ENTER  DISK  NUMBER", Dn 

1117  PRINT  USING  "16X."*Dis*  numoer   -  "",ZZ";Dn 

1120  BEEP 

1123  INPUT  "ENTER  INPUT  MODE  (0-3054A, 1 -FILE) " , Im 

1126  BEEP 

1129  INPUT  "SELECT  HEATING  MODE  (0-ELEC:  1 -WATER) ". Ihm 

1132  BEEP 

1135  INPUT  "ENTER  THERMOCOUPLE  TYPE  ( 0-NEW. 1 -OLD) ", leal 

1138  IF  Im-0  THEN 

1141  BEEP 

1144  INPUT  "GIVE  A  NAME  FOR  THE  RAW  DATA  FILE" .D2_f  lleS 

1147  PRINT  USING  "16X.""Neu  file  name:  "  .  1  4A" :D2_f l leS 

1150  Sizel-20 

1153  CREATE  BDAT  D2_f ileS.Sizel 

1156  ASSIGN  $File2  TO  D2  fileS 

1159! 

1162!  DUMMY  FILE  UNTIL  Nrun  KNOWN 

1165  D1_f ileS-"DUMMY" 

1168  CREATE  BDAT  D1  fileS.Sizel 

1171  ASSIGN  AFilel  TO  D1  fileS 

1174  OUTPUT  -Filet ;DateS 
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1177  IF  Ihm-0  THEM 

1180  BEEP 

1183  INPUT  'ENTER  NUMBER  OF  DEFECTIVE  TCS  (O-DEFAULT) " .Idtc 

1186  IF  Idtc-0  THEN 

1189  Ldtc!-0 

1192  Ldtc2-0 

1195  PRINT  USING  "16X,""No  defective  TCs  exist""" 

1198  END  IF 

1201  IF   Idtc-1    THEN 

1204  BEEP 

1207  INPUT  "ENTER  DEFECTIVE  TC  LOCATION" , Ldtcl 

1210  PRINT  USING  "16X.""TC  is  defective  at  location  "",Du;Ldtc1 

1213  Ldtc2-0 

1216  END  IF 

1219  IF  Idtc-2  THEN 

1222  BEEP 

1225  INPUT  "ENTER  DEFECTIVE  TC  LOCATIONS" .Ldtcl ,Ldtc2 

1228  PRINT  USING  "16X.""TC  are  defective  at  Locations  "  ,D.4X,D" jLdtcl ,Ldtc2 

1231  END  IF 

1234  IF   Idtc>2   THEN 

1237  BEEP 

1240  PRINTER  IS  1 

1243  BEEP 

1246  PRINT  "INVALID  ENTRY" 

1249  PRINTER  IS  701 

1252  GOTO  1180 

1255  END  IF 

1258  END  IF 

1261  OUTPUT  C*File1;Ldtct,Ldtc2 

1264!  Im-1  option 

1267  ELSE 

1270  BEEP 

1273  INPUT  "GIVE  THE  NAME  OF  THE  EXISTING 'DATA  FILE" ,D2_f l le$ 

1276  PRINT  USING  "16X,""0ld  file  name:  " " , 1 4A" ;D2_f i le$ 

1279  ASSIGN  *File2  TO  D2_fileS 

1282  ENTER  SFile2:Nrun 

1285  ENTER  »File2:Dold$ 

1288  PRINT  USING  *16X.""This  data  set  talcen  on  :  "" ,  14A";DoldS 

1291  ENTER  »File2:Ldtc1 .Ldtc2 

1294  IF  Ldtcl >0  OR  Ldtc2>0  THEN 

1297  PRINT  USING  M6X.  ""Thermocouples  were  defective  at  locations: "" ,2(3D,4X) " ; 

Ldtcl ,Ldtc2 

1300  END  IF 

1303  ENTER  »File2;Itt 

1306  END  IF 

1309!  IF  Im-0  AND  Ihm-1  THEN  1595 

1312  BEEP 

1315  INPUT  "WANT  TO  CREATE  A  PLOT  FILE?  ( 0-N, I -Y) " , Iplot 

1318  IF  Iplot-1  THEN 

1321  BEEP 

1324  INPUT  "GIVE  NAME  FOR  PLOT  FILE" ,P_f i leS 

1327  CREATE  BOAT  P_fileS.4 

1330  ASSIGN  3Plot  TO  P  fileS 

1333  END  IF 

1336  IF  Ihm-1  THEN 

1339  BEEP 

1342  INPUT  "UANT  TO  CREATE  Uo  FILE?  ( 0-N, ! -Y) " , Iuf 

1345  IF  Iuf-1  THEN 

1348  BEEP 

1351  INPUT   "ENTER   Uo   FILE   NAME".UfileS 
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1354 

!3S7 

1380 

13S3 

1366 

1369 

I372 

1375 

1378 

1381 

1384 

1387 

1390 

1393 

1396 

1399 

1402 

1405 

1408 

1411 

1414 

1417 

1420 

1423 

1426 

1429 

1432 

1435 

1438 

1441 

1444 

1447 

1450 

1453 

1456 

1459 

1462 

1465 

1468 

1471 

1474 

1477 

1480 

1483 

1486 

1489 

1492 

1495! 

1498! 

1501  ! 

1504! 

1507 

1510 

1513 

1516 

1519! 

1522! 

1525 

1528 

1531 


BDAT  UfUeS.4 
?Ufile  TO  UfileS 


'WANT  TO  CREATE  Re  FILE? 
•1  THEN 


CO-N, l-Y)".Ire 


'ENTER  Re  FILE  NAME" 
BDAT  RefileS.10 
SRefile  TO  RefileS 


,Ref ileS 


IS  1 
THEN 


'Select  tube  number "! 


■6X.""1 

"SX,""2 
■SX.""3 
"6X.""4 

"6X,'"5 
"6X."-6 
"6X."7 
"6X.""3 
"6X.""9 
*6X.'"10 


Smooth  4  inch 
Smooth  4  inch 
Soft  Solder  4 
Soft  Solder  4 
Uieland  Hard  8 


Ret*"'" 

Cu  (Press/Slide) 
inch  Cu""" 
inch  HIGH  FLUX"" 
inch""" 


CREATE 

ASSIGN 

END  IF 

BEEP 

INPUT 

IF  Ire 

BEEP 

INPUT 

CREATE 

ASSIGN 

END  IF 

END  IF 

PRINTER 

IF  Im-0 

BEEP 

PRINT  USING  "4X, 

IF  Ihm-0  THEN 

PRINT  USING  "6X,""0 

PRINT  USING 

PRINT  USING 

PRINT  USING 

PRINT  USING 

PRINT  USING 

PRINT  USING 

PRINT  USING 

PRINT  USING 

PRINT  USING 

PRINT  USING 

PRINT  USING  "5X,""11 

ELSE 

PRINT  USING 

PRINT  USING 

PRINT  USING 

PRINT  USING 

PRINT  USING 

END  IF 

INPUT  Itt 

OUTPUT  AFilel ; Itt 

END  IF 

PRINTER  IS  701 

IF  Itt<10  THEN  PRINT  USING  "16X, 

IF  Itt>9  THEN  PRINT  USING  "16X," 

IF  Ihm-1  THEN  PRINT  USING  "16X," 

BEEP 

INPUT  "ENTER  OUTPUT  VERSION  CO-LONG. 

BEEP 

INPUT  "SELECT  (0-LIQ, 1 -VAP,2-(LIQ+VAP) /2) " , Ilqv 


HIGH  FLUX 

GEUA-K  19 

GEWA-K  26 

GEWA-T  19 

GEWA-T  26  Fins/in" 

THERMOEXCEL-E""" 

THERMOEXCEL-HE""" 


8  inch" 
Fins/ in' 
Fins/in' 
Fins/in" 


"SX,*"0 
"6X.""1 
"6X.""2 
"SX.""3 
"6X,""4 


Smooth  tube""" 
High  Flux""" 
Turbo-B*"" 
High  Flux  Mod""" 
Turoo-3  Mod""" 


"Tube  Number: 
'Tube  Number: 
'Tube  Type: 


"".D":Itt 
"",DD":ltt 

"",15A":TnS(Itt) 


SHORT. 2-N0NE)",Iov 


DIMENSIONS  FO  TESTED  TUBES 

ELECTRIC  HEATED  MODE 

Dl-Diameter  at  thermocouple  positions 

DATA  . 0 1 1 1 1 2S . . 0 1 1 1 1 25 , . 0 1 1 1 1 25 . . 0 1 29540 . . 0 1 2446 . . 0 1 29540 ,. 0 1 0 0 965 

DATA  .01 00965.. 01 157.. 01 157,. 01 157.. 01 157 

READ  D1a(*) 

DI-Dla(Itt) 

D2"Diameter  of  test  section  to  the  base  of  fins 

DATA  .01 5875.. 01 5875.. 01 5375.. 01 5824.. 01 5875,. 01 5824.. 01 270 

DATA  .0127,.0138,.0138..0138..0138 

READ  02a(*> 
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1534! 

1537 1  Di'Tnside  diameter  of  unenhanced  ends 

1540*  DATA  .0127, .0127. .0127, .0132..0127, .0132. .0111125. .01 1 1 12S, .0118, .01 18. .01 

18..  01 111 

IS43  READ  Dia(+> 

1546! 

1549'  Do*0utside  diameter  of  unenhanced  ends 

1552*  DATA  .01 5875,. 01  5375.. 01 5375,. 01 5824.. 01 5875.. 01 5824.. 01 270,. 01 270,. 01 331  . 

.01331.  .01331, .01331 

1555  READ  Doa(*> 

1558! 

1561 '  L"Length  of   enhanced  surface 

1564 "  DATA  . 1016. . 1016, . 1016, . 1016, .2032. .2032. .2032, .2032. .2032, .2032. .2032, .20 

32 

1567  READ  La(*) 

1570! 

1573'  Lii"Length  of    unenhanced  surface  at  the  ends 

1576'  DATA  .0254, .0254, .0254, .0254, .0762. .0762, .0762. .0762. .0762. .0762, .0762, . 07 

62 

1579  READ  Lua(+) 

1582! 

15351  Kcu"Thermal  Conductivity  of  tube 

1588'  DATA  398,344,344.45,344,45,344,944,398.398,398,398 

1591  READ  Kcua(+) 

1594  IF  Ihm-1  THEN 

1597! 

1600 '  Data  statements  for  water  heating  mode 

1603! 

1606  DATA  0.015375,0.015875,0.0169.0.0138.0.0169.0.0,0.0.0 

1609  READ  D2a(*) 

T612  DATA  0.0127,0.0127,0.0145,0.0127,0.0145,0,0,0,0.0 

1615  R£AD  Dia(*) 

1613  DATA  0.015875,0.015875,0.0169,0.015875,0.0169,0.0,0,0,0 

1621  READ  Doa(*) 

1624  DATA  0.3048,0.3048,0.3048.0.3048,0.3048.0.0,0,0,0 

1627  READ  La(*> 

1 630  DATA  0 . 0254 , 0 . 0254 , 0 . 0254 . 0 . 0254 , 0 . 0254 ,0.0,0.0.0 

1633  READ  Lua(*) 

1636  DATA  398,45,398,45,398.0.0,0,0.0 
1639  READ  Kcua(*) 

1642  END  IF 

1645  D2-02a(Itt) 

1648  Di-Dia(Itt) 

1651  Do-Doa(Itt) 

1654  L-La(Itt) 

1657  Lu-Lua(Itt) 

1660  Kcu-Kcua(Itt) 

1663  Xn-.8 

1666  Fr-.3 

1669  IF  Itt-0  THEN  Cf-l.70E+9 

1672  IF  ltt>0  THEN  Cf-3. 7037E+1 0 
1675  A*PI*<Do*2-Di*2)/4 

1673  P-PI*0o 

1681  IF  Ihm-1  THEN 

1634  BEEP 

1637  INPUT   "TUBE    INITIATION   MODE.    (1-HOT  WATER. 2-STEAM.3-C0LD   WATER) ".Itim 
1690      If    Itim-1    THEN   PRINT   USING    "16X."Tube    Initiate:    Hot    Water"" 

1693  IF   Itira-2    THEN  PRINT  USING    "16X."Tube    Initiate:    Steam"" 

1636      IF    Itim-3    THEN   PRINT   USING    "16X."Tube    Initiate:    Cold   Water"" 

1699  INPUT    "TEMP/VEL   MODE:    ( 0-T-CONST, v-DEC: I -T-DEC, V-CONST;    2-T-INC, V-CONST) "  . 

Itv 
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1702  IF  Itv-0  THEM  PRINT  USING  " 1 6X. " "Temp/Vel  Mode:  Constant/Decreasing 

1705  IF  Itv-1  THEN  PRINT  USING  "1SX . " "Temp/Vei  Mode:  Decreasing/Constant 

!708  IF  Itv-2  THEN  PRINT  USING  " 1SX. ""Temp/Vel  Mode:  Increasing/Constant 

1711  INPUT  "WANT  TO  RUN  WILSON  PLOT?  (1 -Y.0-N) " ,Iuil 

1714  IF  Ihm-1  AND  Iuil-0  THEN 

1717  IF  Itt-0  THEN  Ci-.032 

1720  IF  Itt-I  OR  Itt-3  THEN  Ci-.059 

1723  IF  Itt-2  OR  Itt-4  THEN  Ci-.062 

172S  BEEP 

1729  INPUT  "ENTER  CI  (DEF:  WH- . 032. HF- .059, TB-.062) ' ,Ci 

1732  PRINT  USING  "16X. ""Sieder-Tate  """ 

1735  PRINT  USING  "16X,""   Constant   -  "",Z.4D":Ci 

1738  END  IF 

1741  END  IF 

1744  IF  Ihm-1  AND  Im-1  AND  lull-!  THEN 

1747  IF  Itt-0  THEN  Ci-.032 

1750  IF  Itt-1  OR  Itt-3  THEN  Ci-.059 

1753  IF  Itt-2  OR  Itt-4  THEN  Ci-.062 

1756  ASSIGN  *File2  TO  * 

1759  CALL  Hilson(Cf.Ci) 

1762  ASSIGN  SFii«2  TO  D2_fileS 

1765  ENTER  *File2:Nrun,Doid$.Ldtc1 ,Ldtc2, Itt 

1763  END  IF 
1771  Nsub-0 

1774  IF  Idp-4  THEN  Ihm-1 

1777  IF  Ihm-1  THEN  Nsub-8 

1780  Ntc-6 

1783  IF  Ihm-0  THEN  Ntc-12 

1786  J-1 

1789  Sx-0 

1792  Sy-0 

1795  Sxs-0 

1798  Sxy-0 

1801  Repeat:  ! 

1804  IF  Im-0  THEN 

1807  Dtld-2.22 

1810  Ido-2 

1813  ON  KEY  0.15  RECOVER  1801 

1816  PRINTER  IS  1 

1819  PRINT  USING  "4X. ""SELECT  OPTION""" 

1822  PRINT  USING  *SX.""0-TAKE  DATA""" 

1325  IF  Ihm-0  THEN  PRINT  USING  "6X.""1-SET  HEAT  FLUX""* 

1828  IF  Ihm-1  THEN  PRINT  USING  "6X,""1-SET  WATER  FLOW  RATE""" 

1831  PRINT  USING  "6X."*2-SET  Tsat""" 

1834  PRINT  USING  "4X.""N0TE:  KEY  0  -  ESCAPE""" 

1837  BEEP 

1840  INPUT  Ido 

1843  IF  Ido>2  THEN  Ido-2 

1846  IF  Ido-0  THEN  2296 

1 849 ! 

1852!  LOOP  TO  SET  HEAT  FLUX  OR  FLOWMETER  SETTING 

1855  IF  Ido-1  THEN 

1358  IF  Ihm-0  THEM 

1861  OUTPUT  709:"AR  AF12  AL13  VRS" 

1864  BEEP 

1867  INPUT  "ENTER  DESIRED  Qdp" .Dqdp 

1870  PRINT  USING  "4X, ""DESIRED  Qdp  ACTUAL  Qdp""" 

1873  Err-1000 

1876  FOR  1-1  TO  2 

1879  OUTPUT  709: "AS  SA" 
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1382  Sum-0 

1885  FOR  Jt=!  TO  5 

1888  ENTER  709 :E 

1891  Sum-Sum+E 

1894  NEXT  Ji 

1897  IF  1-1  THEN  Volt-Sum*5 

1900  IF  1-2  THEN  Amp-E 

1903  NEXT  I 

190G  Aqdp-Volt*Amp/(PI*D2*L) 

1909  IF  ABS(Aqdp-Dqdp)>Err  THEN 

1912  IF  Aqdp>0qdp  THEN 

1915  BEEP  4000,. 2 

1918  BEEP  4000, .2 

1921  BEEP  4000,. 2 

1924  ELSE 

1927  BEEP  250,. 2 

1930  BEEP  250,. 2 

1933  BEEP  250,. 2 

193S  END  IF 

1939  PRINT  USING  MX,MZ.3DE,2X,MZ.3DE" :Dqdp,Aqdp 

1942  WAIT  2 

1945  GOTO  187S 

1948  ELSE 

1951  BEEP 

1954  PRINT  USING  MX,MZ.3DE.2X,MZ.3DE" ;DqdP,Aqdp 

1957  Err-500 

1960  WAIT  2 

1963  GOTO  187S 

1966  END  IF 

1969  ELSE 

1972  BEEP 

1975  INPUT  "ENTER  FLOWMETER  SETTING" . Fins 

1978  GOTO  1819 

1981  END  IF 

1984  END  IF 

1987! 

1990!  LOOP  TO  SET  Tsat 

1993  IF  Ido-2  THEN 

1996  IF  Ikdt-1  THEN  2011 

1999  BEEP 

2002  INPUT  "ENTER  DESIRED  Tsaf.Dtld 

2005!  PRINT  USING  "4X,""  DTsat  ATsat   Rate    Tv     Rate""" 

2008  Ikdt-1 

2011  Old1-0 

2014  01d2-0 

2017  Nn-1 

2020  Nrs-Nn  MOD  15 

2023  Nn-Nn+1 

2026  IF  Nrs-1  THEN 

2029  PRINT  USING  MX.""   Tsat    Tldl     Tld2     Tv    Tsump   Tinlet    Tpile 

Tout""" 

2032  END  IF 

2035  IF    Ihm-0    THEN   OUTPUT   709: "AR   AF3   AL11    VRS" 

2038  IF    Ihm-1    THEN   OUTPUT   709: "AR   AFO    AL5   VRS" 

2041  FOR    1-1    TO   6 

2044  IF    Ihm-0   AND    I>4   THEN   2125 

2047  Sum-0 

2050  OUTPUT    709: "AS   SA" 

2053  FOR   Ji-1    TO   20 

2056  ENTER   709:Eliq 
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2059  Sum-Sum+Eliq 

2062  NEXT  Ji 
2065  Eliq-Sum/20 

2063  Tld-FNTvsv(Eliq) 
2071  IF  1-1  THEN  Tldl-Tld 
2074  IF  1-2  THEN  Tld2»Tld 
2077  IF  1-3  THEN  Tv-Tld 
2080  IF  1-4  THEN  Tsump-TLd 
2083  IF  1-5  THEN  Tinlet-Tld 
2086  IF  1-6  THEN  Tout-Tld 
2089  NEXT  I 

2092  IF  Ihm-1  THEN 

2095  OUTPUT  709: "AR  AF20  AL20  VR5" 

2098  OUTPUT  709: "AS  SA" 

2101  Sum-0 

2104  FOR  Klc-1  TO  20 

2107  ENTER  709:£ 

2110  Sum-Sum+E 

2113  NEXT  Kk 

2116  Erat ( 7) -ABS( Sum/20) 

2119  Tpile-Emt(7)/3.96E-4 

2122  END  IF 

2125  Atld-<Tld1*TLd2)*.S 

2128  IF  ABS(Atld-Dtld)>.2  THEN 

2131  IF  Atld>Dtld  THEN 

2134  BEEP  4000,. 2 

2137  BEEP  4000,. 2 

2140  BEEP  4000,. 2 

2143  ELSE 

2146  BEEP  250.. 2 

2149  BEEP  250.. 2 

2152  BEEP  250,. 2 

2155  END  IF 

2158  Err1-Atld-0ld1 

2161  Oldl-Atld 

2164  Err2-Tv-01d2 

2167  01d2-Tv 

2170  IF  Ihm-0  THEN  PRINT  USING  *4X.S(ttDD.DD.2X) * tDtld, Tldl , Tld2 . Tv. Tsump 

2173  IF  Ihm-1  AND  Idp-0  THEN  PRINT  USING  "4X. 7(ttDD.DD,2X) " :Dt Id. Tldl , Tld2. Tv , Ts 

ump . Tiniet . Tpi le 

2176*  IF  Ihm-1  AND  Idp-4  THEN  PRINT  USING  "4X,S(MDD.DD,2X) ,3(M3D.DD,2X) " ;Dtld, Tl 

d1 , Tld2, Tv. Tsump, Tin 

2179  UAIT  2 

2182  GOTO  2020 

2185  ELSE 

2188  IF  ABS(Atld-Dtld)>.1  THEN 

2191  IF  Atld>Dtld  THEN 

2194  BEEP  3000.. 2 

2197  BEEP  3000. .2 

2200  ELSE 

2203  BEEP  800.. 2 

2206  BEEP  800,. 2 

2209  END  IF 

2212  Errl-Atld-Oldl 
2215  Oldl-Atid 

2213  Err2-Tv-Qld2 
2221  01d2-Tv 

2224  IF    Ihm-0    THEN   PRINT   USING    "4X.5( MDD.DD.2X) * :Dtld. Tldl . Tld2. Tv, Tsump 

2227*   IF    Ihm-1    THEN   PRINT   USING    *4X ,S(MDD.DD ,2X) ,3(M3D.DD. 1X) " ;Dtld,Tld1 . Tld2, Tv 

. Tsump, Tiniet.Tpile, 
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2230  WAIT  2 

2233  GOTO  2020 

223S  ELSE 

2239  BEEP 

2242  Errl-Atld-Oldl 

2245  Oldl-Atld 

2248  Err2-Tv-Qld2 

2251  01d2-Tv 

2254  IF  Ihii-0  THEN  PRINT  USING  "4X,S(MDD.DD,2X)  "  :D  tld,  Tidl  ,Tld2.  Tv,  Tsump 

2257  IF  Ihm-1  THEN  PRINT  USING  *4X.3(MDD.DD ,2X) " ;Dtld, TLd1 , Tld2. Tv, Tsump, Tinlet 

, Tp ile, Tout 

2260  WAIT  2 

22S3  GOTO  2020 

22S6  END  IF 

2269  END  IF 

2272  END  IF 

2275!  ERROR  TRAP  FOR  Ido  OUT  OF  BOUNDS 

2278  IF  Ido>2  THEN 

2281  BEEP 

2284  GOTO  1819 

2287  END  IF 

2290! 

2293!  TAKE  DATA  IF  Im-0  LOOP 

2296  IF  Ikoi-1  THEN  2308 

2299  BEEP 

2302  INPUT   "ENTER   BULK  OIL  7.", Bop 

2305  Ikoi-1 

2308  IF  Ihra-0  THEN  OUTPUT  709: "AR  AFO  ALII  VRS" 

2311  IF  Ihm-1  THEN  OUTPUT  709: "AR  AFO  ALS  VRS" 

2314  IF  Ihm-0  THEN  Ntc-12 

2317  FOR  1-1  TO  Ntc 

2320  OUTPUT  709: "AS  SA" 

2323  Sum-0 

2326  FOR  Ji-1  TO  20 

2329  ENTER  709 :E 

2332  Sum-Sum+E 

2335  IF   I-(9-Nsub)    OR   X 1 O-Nsubr  THEN  Et( Ji-t ) -E 

2333  NEXT   Ji 
2341  Kdl-0 

2344  IF   I-<9-Nsub)    OR    I-UO-Nsub)    THEN 

2347  Eave-Sum/20 

2350  Sum-0. 

2353  FOR  Jk-0  TO  19 

2356  IF  ABS(Et(Jk)-EaveK5.0E-S  THEN 

2359  Sum-Sum+EtUk) 

2362  ELSE 

2365  Kdl-Kdi+1 

2368  END  IF 

2371  NEXT  Jk 

2374  IF  I-(9-Nsub)  OR  I-UO-Nsub)  THEN  PRINT  USING  *4X,""Kdl  -  "",DD":Kdl 

2377  IF  Kdl>10  THEN 

2380  BEEP 

2383  BEEP 

2386  PRINT  USING  "4X.""Too  much  scattering  in  data  -  repeat  data  set""" 

2389  GOTO  1816 

2392  END  IF 

2395  END  IF 

2398  Emt(I)-Sum/(20-Kdl) 

2401  NEXT  I 

2404  IF  Ihm-1  THEN 

2407  OUTPUT  709: "AR  AF20  AL20  VRS" 
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2590  NEXT  I 

2533  Tu»-Tua/(8-Idtc) 

2536  END  IF 

2599  Tld-T(9-Nsub) 

2602  Tld2-T(10-Nsub) 

2605  Tlda-(Tld+Tld2)+.S 

2608  Tv-T(H-Nsub) 

2611  IF  Itt<3  AND  Ihm-0  THEN 

2614  Tld2--99.99 

2617  Tv-<T(10)+T( 11 ) )/2 

2620  END  IF 

2623  Tsump-T(12-Nsub) 

2626  IF  Ihm-0  THEN  2635 

2629  Tinlet-T(13-Nsub) 

2532  Tout-T(14-Nsub) 

2635  IF  Ihm-0  THEN 

2638  Amp-ABS(Ir) 

2641  Volt-ABS(Vr)*2S 

2644  Q-Volt*Amp 

2647  END  IF 

2650  IF  Itt-0  AND  Ihm-0  THEN 

2653  Kcu-FNKcu(Tw) 

2656  ELSE 

2659  Kcu-Kcua(Itt) 

2562  END  IF 

2665! 

2668!  FOURIER  CONDUCTION  EQUATION  WITH  CONTACT  RESISTANCE  NEGLECTED 

2671  IF  Ihm-0  THEN  Tu»-Tw-Q*L0G(  D2/D1  )  /(2*PI*Kcu*L) 

2674  IF  Ilqv-0  THEN  Tsat-Tlda 

2677  IF  Ilqv-1  THEN  Tsat-( Tlda+Tv)*.S 

2680  IF  Ilqv-2  THEN  Tsat-Tv 

2633  IF  Ihm-1  THEN 

2686  Tavg-TinLet 

2689  Grad-37.98S3+.104388*Tava 

2692  Tdrop-ABS(Emt  (7)  )*1  .E+6/U0*Grad) 

2695  Tavgc-Tinlet-Tdrop*.5 

2698  IF  ABS(Tavg-Tavgc)>.01  THEN 

2701  Tavg-( Tavg+Tavgc J*.S 

2704  GOTO  2689 

2707  END  IF 

2710! 

2713!  COMPUTE  UATER  PROPERTIES 

2716  IF  Ihm-1  THEN 

2719  Ku-FNKw(Tavg) 

2722  Muua-FN(1uiii(  Tavg) 

2725  Cpu-FNCpu(Tavg) 

2728  Pruj-FNPru(Tavg) 

2731  Rhow-FNRhou(Tavg) 

2734  Twi-Tavg 

2737! 

2740!  Compute  MDOT 

2743  Mdot-3 . 9657E-3+Fms* ( 3 . 6 1 9S5E-3-Fms* (  8 . 32006E-6-Fins*( 1 . 23688E-7-Fms*4 . 31 897 

E-10) ) ) 

2746!  Mdot-Mdot*< 1 .0365-Tinlet*< 1 .96644E-3-Tiniet*5.2S2E-6) ) /I .0037 

2749  Kdt-0 

2752  Q-f1dot*Cpw*Tdrop 

2755  Lmtd-Tdrop/L0G( ( Tinlet-Tsat) /( Tinlet-Tdrop-Tsat) ) 

2753  Uo-Q/(PI*0o*L*Lmtd) 

276!  Ru-0o*LQG(Do/Di)/(2.*Kcu) 

2764  Tu-Tsat+Fr+Lmtd 
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2410  OUTPUT  709: "AS  SA" 

2413  Sum-0 

2416  FOR  Klc-1  TO  20 

2419  ENTER  709:E 

2422  Sum-Sum+E 

2425  NEXT  Kk 

2428  Emf (7)-ABS(Sum)/20 

2431  END  IF 

2434  IF  Ihm-0  THEN 

2437  OUTPUT  709: "AR  AF12  AL13  VRS" 

2440  FOR  1-1  TO  2 

2443  OUTPUT  709;  "AS  SA" 

2446  Sum-0 

2449  FOR  Ji-1  TO  2 

2452  ENTER  709 :£ 

2455  Sum-Sum+E 

2458  NEXT  Ji 

2461  IF   1-1    THEN  Vr«Sum/2 

2464  IF   1-2   THEN   Ir-Sum/2 

2467  NEXT   I 

2470  END   IF 

2473  ELSE 

2476  IF   Ihm-0   THEN  ENTER   »Fi  le2:Bop,ToldS.Emf  (*)  ,Vr .Ir 

2479  IF   Ihm-1    THEN   ENTER   9File2;Bop.  ToldS. Emf  (*), Fins 

2482  END   IF 

2485! 

2488!  CONVERT  emf'S    TO    TEMP, VOLT. CURRENT 

2491  Twa-0 

2494  .  FOR   1-1    TO  Ntc 

2497  IF   IdtOO   THEN 

2500  IF   I-Ldtc!    OR    I-Ldtc2   THEN 

2503  T(D— 99.99 

2506  GOTO  2536 

2509  END  IF 

2512  END  IF 

2515  IF  Itt<4  AND  Ihm-0  THEN 

2513  IF  I>4  AND  I<9  THEN 
2521  T(D— 99.99 

2524  GOTO  2536 

2527  END  IF 

2530  END  IF 

2533  T(I)-FNTvsv(Emf (I) ) 
2536  NEXT  I 

2539  IF  Itt<4  THEN   - 

2542  FOR  1-1  TO  4 

2545  IF   I-Ldtcl    OR    I-Ldtc2   THEN 

2548  Tua-Tua 

255 1  ELSE 

2554  Twa-Tuia+TU) 

2557  END  IF 

2560  NEXT  I 

2563  Tua-Tua/ (4-Idtc) 

2566  ELSE 

2569  IF  Ihm-1  THEN  2599 

2572  FOR  1-1  TO  3 

2575  IF  I-Ldtc 1  OR  I-Ldtc2  THEN 

2578  Twa-Tua 

2531  ELSE 

2534  Tuia-Tua+T(I) 
2587  END   IF 
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2767  Vu-fldot/(Rhou»*PI*0i*2/4) 

2770  Reu-Rhou*^iii*0  L/Muua 

2773  Hi-Ci*Ku/Di*Reu"  .8*Prw*  (1/3. ) +  (Muuia/FN«uu(  Tui )  )  '  .14 

2776  Twic-Tavg-Q/(PI*Oo*L*Hi) 

2779  IF  ABS(Tui-Twic)>.01  THEN 

2782  Tu»i-(Tu>i  +  Twic)*.S 

2735  GOTO  2773 

2788  END  IF 

2791  Tu»i-(Tii»i+Tuuc)*.5 

2794  Ho-1/<l/Uo-Do/(Di*Hi)-Rw) 

2797  END  IF 

2800  END  IF 

2803  IF  Ihm-1  THEN 

2806  Thetab-Q/(Ho*PI*0o*L) 

2809  Tw-Tsat+Thetab 

2812  ELSE 

2815  Thetab-Tu-Tsat 

2813  END  IF 

2821  IF  Thetab<0  THEN 

2824  BEEP 

2827  INPUT  'TUALL<TSAT  (0-CONTINUE.  t -END)  ".lev 

2830  IF  Iev-0  THEN  GOTO  1804 

2833  IF  lev- I  THEN  3130 

2836  END  IF 
2839! 

2842!  COMPUTE  VARIOUS  PROPERTIES 

2845  Tfilm-(Tw+Tsat)*.5 

2848  Rho-FNRho(Tfiim) 

2851  Mu-FNMu(Tf ilm) 

2854  K-FNK(Tfilm) 

2857  Cp-FNCp(Tf ilm) 

2860  8eta-FMBeta( Tf ilm) 

2863  Hfg-FNHfg(Tsat) 

2866  Ni-Mu/Rho 

2863  Alpha-K/(Rho*Cp) 

2872  Pr-Nt/Alpha 

2875  Psat-FNPsat(Tsat) 

2873! 

2881!  COMPUTE  NATURAL-CONvECTIvE  HEAT- TRANSFER  COEFFICIENT 

2884!  FOR  UNENHANCED  END(S) 

2837  Hbar-!90 

2890  Fe-(Hbar+P/(Kcu*A) >* .S*Lu 

2893  Tanh-FNTanh(Fe) 

2896  Theca-Thetab*Tanh/Fe 

2899  Xx-(9.81*8eta+Thetab+0o-3*Tanh/(Fe*Ni*Alpha) ) * .166667 

2902  Yy-(1+( .559/Pr)* (9/16)) "(8/27) 

2905  Hbarc-K/Do*( .S+.387*Xx/Yy) '2 

2908  IF  ABS( (Hbar-Hbarc) /Hbarc) > . 001  THEN 

2911  Hbar-(Hbar+Hbarc)*.S 

2914  GOTO  2890 

2917  END  IF 

2920! 

2923!  COMPUTE  HEAT  LOSS  RATE  THROUGH  UNENHANCED  END(S) 

2926  QI-(Hbar+P*Kcu*A) * .S+Thetab*Tanh 

2929  Qc-Q-2*Q1 

2932  As-PI*02*L 

2935! 

2938!  COMPUTE  ACTUAL  HEAT  FLUX  AND  BOILING  COEFFICIENT 

2941  Odp-Qc/As 

2944  Htube-Odp/Thetab 
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2947  Csf-(Cp*Thetab/Hfg)/(Qdp/(Mu*Hfg>+( .Q14/(9.81*Rho) * .5) " ( 1/3. )*Pr" 1 . 7) 

2350! 

2953!  RECORD  TIME  OF  DATA  TAKING 

2956  IF  Im-0  THEN 

2959  OUTPUT  709;  "TD" 

2962  ENTER  709:Told$ 

2965  END  IF 

2968! 

2971!  OUTPUT  DATA  TO  PRINTER 

2974  PRINTER  IS  701 

2977  IF  Iov-0  THEN 

2980  PRINT 

2983  PRINT  USING  "1  OX. ""Data  Set  Number  -  ' 

4A":J,8op.ToidS 

2986  IF  Ihm-0  THEN 

2989  PRINT  USING  *10X,""TC  No:    1      2 

p  ■  m  m 

2992  PRINT  USING  "10X,""Temp  : "• ,8( 1X,MDD.DD) " ; T( 1 ) , T(2) , T(3) ,T(4) , T(5) . T(6)  ,  T( 

7) ,T(8) 

2995  PRINT  USING  "10X.""  Twa  Tliqd   Tliqd2   Tvapr  Psat   Tsump""" 

2998  PRINT  USING  ,10X,2(MDD.DD.  1X) ,  1X.MDD.DD. 1X,2( 1X.MDD.DD) .2X.MDD.D" ;  Tu,  Tld.T 

ld2, Tv.Psat.Tsump 

3001   PRINT  USING  "10X.""  Thetab  Htube      Qdp""" 

3004  PRINT  USING  "1  OX, MDD. 3D,  1X.rtZ.3DE. 1X.MZ.3DE" ; Thetab, Htube, Qdp 

3007  ELSE 

3010  PRINT  USING  "10X,""  Fms    Vu»   Tsat   Tinl  Tdrop  Thetab     q         Uo 

Ho""" 
3013*  PRINT  USING  "1 0X,4(2D.DD. IX) .Z.3D, IX,DD.DD. 1X,3(rtZ.3DE, IX) " ;Fms, Vwdsat, Ti 
nlet, Tdrop, Thetab.Qd 


,DDD.2X,""Bullc  Oil  7. 
3      4      5 


\DD.D,SX,1 

6      7 


3016 

3019 

3022 

3025 

3028 

3031 

3034 

3037 

3040 

3043 

3046 

3049 

3052 

tab 

3055 

3053 

hetab 

3061 

3064 

3067 

3070 

3073 

3076 

3079 

3082 

3085 

3088 

3091 

3094 

3097 

3100 

3103 


END  IF 

END  IF 

IF  Iov-1  THEN 

IF  J-1  THEN 

PRINT 

IF  Ihm-0  THEN 

PRINT  USING  "10X.*"  RUN  No  Oil/.   Tsat    Htube      Qdp       Thetab""" 

ELSE 

PRINT  USING  "10X.""   FMS    OIL/.   TSAT    HTUBE      QDP       THETAB""" 

END  IF 

END  IF 

IF  Ihm-0  THEN 

PRINT  USING  " 1 2X , 3D . 4X , DD . 2X . MDD . DD . 3 ( 1 X , MZ . 3DE ) " ; J . Bop , Tsa t , H  t ube , Qdp ,  The 


ELSE 

PRINT  USING 


'  1 2X , 3D , 4X . DD . 2X , MDD . DD , 3 ( 1 X , MZ . 3DE ) " ; Fms , Bop , Tsa  t , H tube , Qdp . T 


END  IF 

END  IF 

IF  Im-0  THEN 

BEEP 

INPUT  "OK  TO  STORE  THIS  DATA  SET  < 1-Y, 0-N)?" ,0k 

END  IF 

IF  Ok-1  OR  Im-1  THEN  J-J-H 

IF  Ok-1  AND  Im-0  THEN 

IF  Ihm-0  THEN  OUTPUT  $Fi lei :Bop, ToldS.Emt (*> , Vr ,Ir 

IF  Ihm-1  THEN  OUTPUT  9File1  ;Bop, ToidS.Emf (* ) ,Fms 

END  IF 

IF    Iuf-1    THEN  OUTPUT  $Ufile:Vw.Uo 

IF    Ire-1    THEN   OUTPUT   iSRet  ile;Fms.Reu 

IF    dm- 1    OR   Ok-!)    AND    Iplot-I    THEN   OUTPUT   SP  lot*.  Qdp,  Thetab 

IF   Im-0    THEN 
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3106  BEE? 

3109  INPUT  "WILL  THERE  BE  ANOTHER  RUN  (  1  -Y  .  0-N I ?" ,Go_on 

3112  Nrun-J 

31 15  IF  Go  on-0  THEN  3130 

3118  IF  Go_on<>0  THEN  Repeat 

3121  ELSE 

3124  IF  J<Nrun+1  THEN  Repeat 

3127  END  IF 

3130  IF  Im-0  THEN 

3133  BEEP 

3136  PRINT  USING  "1  OX. ""NOTE:  "".ZZ.""  data  runs  were  stored  in  file  ",10A";J- 
1  ,D2_fiieS 

3139  ASSIGN  9File1  TO  * 

3142  OUTPUT  *File2:Nrun-1 

3145  ASSIGN  SFilel  TO  D1_f ileS 

3148  ENTER  9Fi lei :DateS.Ldtcl ,Ldtc2. Ett 

3151  OUTPUT  9File2; DateS, Ldtcl ,Ldtc2.Itt 

3154  FOR  1-1  TO  Nrun-1 

3157  IF  Ihoi-0  THEN 

3160  ENTER  3Fi lei ;3op, ToldS.Emf (*) , Vr . Ir 

3163  OUTPUT  *File2:Bop,Told$,Emf(*),Vr,Ir 

3166  ELSE 

3169  ENTER  JFi lei :Bop. ToldS.Emf (*> .Fms 

3172  OUTPUT  $File2;Bop, ToldS.Emf (»),Fms 

3175  END  IF 

3178  NEXT  I 

3181  ASSIGN  SFilel  TO  * 
3184  PURGE  •DUMMY- 
SI  87  END  IF 
3130  BEEP 
3193  PRINT 
3196*  IF  Iplot-1  THEN  PRINT  USING  "  1  OX.  ""NOTE:  ",11."'    X-Y  pairs  were  stored  in 

plot  data  file  "",1 

3199  ASSIGN  *File2  TO  * 

3202  ASSIGN  *Plot  TO  * 

3205  IF  Iuf-I  THEN  ASSIGN  ?Ufile  TO  + 

3208  IF  Ire-1  THEN  ASSIGN  SRefile  TO  * 

3211  CALL  Stats 

3214  BEEP 

3217  INPUT  'LIKE  TO  PLOT  DATA  ( I-Y.0-N)?" ,0k 

3220  IF  0k-1  THEN  CALL  Plot 

3223  SUBEND 
3226! 

3229!  CURVE  FITS  OF  PROPERTY  FUNCTIONS 

3232  DEF  FNKcu(T) 

3235!  0FHC  COPPER  250  TO  300  K 

3238  Tk-T+273.15     !C  TO  K 

3241  K-434-.!12*Tic 

3244  RETURN  K 

3247  FNEND 

3250  DEF  FNMu(T) 

3253!  170  TO  360  <  CURVE  FIT  OF  VISCOUSITY 

3256  Tk-T+273.1S     !C  TO  K 

3259  t1u-EXP<-4.4636+(1011  .47/Tk)  )*1  .OE-3 

3262  RETURN  Mu 

3265  FNEND 

3268  DEF  FNCp(T) 

3271!  130  TO  400  K  CURVE  FIT  OF  Cp 

3274  Tlc-T+273.15     !C  TO  K 

3277  Cp-. 40188+1  .S5007E-3*Tk-H  .SI494E-6+Tk  * 2-6 . S7853E- 1 0*Tlc  "3 
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3230 

3233 

3236 

3289 

3292 

3295 

3298 

3301 

3304 

3307 

3310 

3313 

3316 

3319 

3322 

3325! 

3328 

3331 

3334 

3337 

3340 

3343 

3346 

3349 

3352 

3355 

3358 

3361 

3364 

3367 

3370 

3373 

3376 

3379 

3382 

3385 

3388 

3391 

3394 

3397 

3400 

3403 

3406 

3409 

3412 

3415 

3413 

3421 

3424 

3427! 

3430 

3433 

3436 

3439 

3442 

3445 

3448 

3451 

3454 

3457 


,5+1 .  119828*X"2 


Icai 


Cp»Cp*1000 

RETURN  Cp 

FNEND 

DEF  FNRho(T) 

Tk-T+273. 1S     !C  TO  < 

X-1-(1.8*Tlc/753.95)   !K  TO  R 

Ro-36. 32+61 .  14641 4*X*( 1/3) +  16.41 8015*X  +  1 7. 476838*X* 

Ro-Ro/. 062423 

RETURN  Ro 

FNEND 

DEF  FNPr(T) 

Pr-FNCp <T) *FNMu ( T ) /FNK ( T) 

RETURN  Pr 

FNEND 

DEF  FNK(T) 

T<360  K  WITH  T  IN  C 

K-.071-.Q00261+T 

RETURN  K 

FNEND 

DEF  FNTanh(X) 

P-EXP(X) 

Q-1/P 

Tanh-(P-Q)/(P+Q) 

RETURN  Tanh 

FNEND 

DEF  FNTvsv(V) 

COM  /Cc/  C(7) 

T-C(0) 

FOR  1-1  TO  7 

T-T+C(I)*V*I 

NEXT  I 

IF  Ical-1  THEN 

T-T-6.7422934E-2+T+(  9.  0277043E-3-T+( -9.325991 7E-5) ) 

ELSE 

T-T+8.S26897E-2+T*(3.76199E-3-T+5.0689259E-5) 

END  IF 

RETURN  T 

FNEND 

DEF  FNBeta(T) 

Rop-FNRho(T+.1) 

Rom-FNRho(T-. I ) 

Beta— 2/  (  RoD+Rom )  *  (  Rop-Rom )  /  .  2 

RETURN  Beta 

FNEND 

DEF  FNHfg(T) 

Hfg-1.3741344E+5-T*(3.3094361E+2+T*1. 2165143) 

RETURN  Hfg 

FNEND 

DEF  FNPsat(Tc) 

0  TO  30  deg  F  CURVE  FIT  OF  Psat 

Tf-1.8*Tc+32 

Pa-5.945S2S+Tf*( . 15352082+Tf*( 1 .4840963E-3+Tf*9.6150671E-6) ) 

Pg-Pa-14.7 

IF  Pg>0  THEN 

Psat-Pg 

ELSE 

Psat»Pg*29. 92/14. 

END  IF 

RETURN  Psat 

FNEND 


!  +-PSI6,  — in  Ha 
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3460  DEF   FNHsmoothU.3op.Isat) 

3463  DIM  A(5) ,B(5) ,C(5) .0(5) 

3466  DATA  .20526. .25322. .319048. .55322. . 79909. 1 . 00258 

3469  DATA  .7451 5. . 72992. . 731 99, . 71 225 . .68472. .54197 

3472  DATA  .41 092. . 1 7726. .25142. .54806. .31916 . 1 .0845 

3475  DATA  .71403. . 72913. . 72565. .696691 , .665867 . .61 889 

3478  READ  A(*),B<*) ,C(*) ,D<*) 

3481  IF  Bop<6  THEN  I-Bop 

3484  IF  Bop-6  THEN  1-4 

3487  IF  Bop- 10  THEN  1-5 

3490  IF  Isat-1  THEN 

3493  Hs-EXP(A(I)+B(I)*LOGU)) 

3496  ELSE 

3499  Hs»EXP(C(I)+D(I)*LOGU)) 

3502  END  IF 

3505  RETURN  Hs 

3508  FNEND 

3511  DEF  FNPoly(X) 

3514  COM  /Cply/  A(  10, 1 0) ,C( 1 0) .B(4) .Nop, Iprnt.Opo, Ilog.If n, I join.Njoin 

3517  Xl-X 

3520  Poly-B(O) 

3523  FOR  1-1  TO  Nop 

3526  IF  Ilog-1  THEN  Xl-LQG(X) 

3529  Poly-Poly+B(I)*Xl*I 
3532  NEXT  I 

3535  IF  Ilog-1  THEN  Poly-EXP(Poly ) 

3538  RETURN  Poly 
3541  FNEND 

3544  SUB  Poly 

3547  DIM  R(  TO)  .St  TO)  ,Sy(  12)  ,Sx(  12)  ,Xx(  100)  .YydOO) 

3550  COM  /Cply/  A(  1  0  . 1  0  )  ,C(  1  0)  ,B(  4)  ,N.  Iprnt  ,0po.  Hog, If  n,  I  join.Njoin 

3553  COM  /Xxyy/  Xp< 25) , Yp( 25 ) 

3556  FOR  1-0  TO  4 

3559  8(I)-0 

3562  NEXT  I 

3565  BEEP 

3568  INPUT  "SELECT  ( 0-FILE. 1 -KEYBOARD. 2-PROGRAM)  \  Im 

3571  Im-Im+1 

3574  BEEP 

3577  INPUT  "ENTER  NUMBER  OF  X-Y  PAIRS". Mo 

3530  IF  Im-1  THEN 
3583  BEEP 

3536  INPUT  "ENTER  DATA  FILE  NAME" ,D_f ileS 

3539  BEEP 

3592  INPUT  "LIKE  TO  EXCLUDE  DATA  PAIRS  ( 1 -Y.O-N)?' . led 

3595  IF  Ied-1  THEN 

3598  BEEP 

3601  INPUT  "ENTER  NUMBER  OF  PAIRS  TO  BE  EXCLUDED" , Ipex 

3604  END  IF 

3607  ASSIGN  SFile  TO  D  fileS 

3610  ELSE 

3613  BEEP 

3616  INPUT  "WANT  TO  CREATE  A  DATA  FILE  ( 1 -Y ,0-N) ?" , Yes 

3619  IF  Yes-1  THEN 

3622  BEEP 

3625  INPUT  "GIVE  A  NAME  FOR  DATA  FILE" ,D_f ileS 

3628  CREATE  BOAT  D_fileS,S 

3631  ASSIGN  SFiie  TO  D_fileS 

3634  END  IF 

3637  END  IF 
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3640  BEEP 

3643  INPUT  "ENTER  THE  ORDER  OF  POLYNOMIAL" ,N 

3646  FOR  1-0  TO  N*2 

3649  Sy(I)-0 

3652  Sx(I)-0 
3655  NEXT  I 

3653  IF  led- I  AND  Im-1  THEN 
3661  FOR  1-1  TO  Ipex 

3664  ENTER  3File;X,Y 

3667  NEXT  I 

3670  END  IF 

3673  FOR  1-1  TO  Np 

3676  IF  Im-1  THEN 

3679  IF  Opo-2  THEN  ENTER  9File;X,Y 

3632  IF  0po<2  THEN  ENTER  «File;Y,X 

3635  IF  Opo-1  THEN  Y-Y/X 

3688  IF  Ilog-1  THEN 

3691  IF  Opo-2  THEN  Xt-X/Y 

3694  X-LOG(X) 

3697  IF  Opo-2  THEN  Y-LOG(Xt) 

3700  IF  0po<2  THEN  Y-LOG(Y) 

3703  END  IF 

3706  END  IF 

3709  IF  Im-2  THEN 

3712  BEEP 

3715  INPUT  "ENTER  NEXT  X-Y  PAIR",X.Y 

3718  IF  Yes-1  THEN  OUTPUT  «File;X.Y 

3721  END  IF 

3724  IF  Im<3  THEN 

3727  Xx(I)-X 

3730  Yy(I)-Y 

3733  ELSE 

3736  X-Xp(I-1)  •  . 

3739  Y-Yp(I-l) 

3742  END  IF 

3745  R(0)-Y 

3748  Sy(0)-Sy(0)+Y 

3751  S(l)-X 

3754  Sx(1)-Sx(1)+X 

3757  FOR  J-1  TO  N 

37S0  R(J)-R(J-1 )*X 

3763  Sy(J)-Sy(J)+R(J) 

3766  NEXT  J 

3769  FOR  J-2  TO  N*2 

3772  SU)-SU-1)*X 

3775  Sx(J)-Sx(J)*SU) 

3778  NEXT  J 

3781  NEXT  I 

3784  IF  Yes-1  AND  Im-2  THEN 

3787  BEEP 

3790  PRINT  USING  "12X,DD.""  X-Y  pairs  were  stored  in  file  "" . 1 OA" ;Np.D_f ileS 

3793  END  IF 

3796  Sx(0)-No 

3799  FOR  1-0  TO  N 

3802  C(I)-Sy(I) 

3805  FOR  J-0  TO  N 

3808  A(I,J)-Sx(I+J) 

3811  NEXT  J 

3814  NEXT  I 

3817  FOR  1-0  TO  N-1 

3820  CALL  Divide(I) 
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3823  CALL  Subtract(  J>1 ) 

3826  NEXT  I 

3829  B(N)-C(N)/A(N,N> 

3832  FOR  1-0  TO  N-1 

3835  B(N-1-I)-C(N-1-I) 

3838  FOR  J-0  TO  I 

3341  B(N-l-I)-B(N-1-I)-A(N-t-I.N-J)*B(N-J) 

3844  NEXT  J 

3847  B(N-I-I)-B(N-1-I)/A(N-1-I,N-1-I) 

3850  NEXT  I 

3853  1PRINTER  IS  701 

3856  ! PRINT  B(*) 

3859  [PRINTER  IS  705 

3862  IF  Iprnt-0  THEN 

3865  PRINT  USING  "12X, ""EXPONENT    COEFFICIENT" 

3868  FOR  1-0  TO  N 

3871  PRINT  USING  * 1SX,DD.SX,MD. 7DE" ;I,B( I) 

3874  NEXT  I 

3877  PRINT  "  " 

3880  PRINT  USING  "12X.""DATA  POINT      X  Y      Y(CALCULATED)  DISCR 

EPANCY""" 

3883  FOR  1-1  TO  Np 

3886  Yc-B(O) 

3889  FOR  J-!  TO  N 

3892  Yc-Yc+8(J)*Xx(I)'J 

3895  NEXT  J 

3898  D-Yy(I)-Yc 

3901  PRINT  USING  "1SX,3D.4X.4<MD.SDE, IX) " ; I,Xx( I) , Yy ( I) , Yc.D 

3904  NEXT  I 

3907  END  IF 

3910  ASSIGN  SFile  TO  * 

3913  SUBEND 

3916  SUB  Divide(M) 

3919  COM  /Cpiy/  A( 10, 1 0 ) ,C< 1 0 ) ,B( 4 ) ,N, Iprnt.Opo. Hog, If n, I join.N join 

3922  FOR  I-M  TO  N 

3925  Ao-A(I.M) 

3928  FOR  J-M  TO  N 

3931  A(I.J)-A(I,J)/Ao 

3934  NEXT  J 

3937  C(I)-C(I)/Ao 

3940  NEXT  I 

3943  SUBEND 

3946  SUB  Subtract(K) 

3949  COM  /Cply/  AU0.10)  ,C(10) ,B(4) ,N, Iprnt ,Opo,Ilog, If n. I/oin.N join 

3952  FOR  I-K  TO  N 
3955  FOR  J-K-1  TO  N 

3953  A(I.J)-A(K-1 ,J)-A(I.J) 
3961  NEXT  J 

3964  C(I)-C(K-1)-C(I) 

3967  NEXT  I 

3970  SUBEND 

3973  SUB  Plin 

3976  COM  /Cply/  AdO.10)  ,C(  TO)  ,B(4)  ,N.  Iprnt  ,0po.  Hog,  If  n,  I  join.N  join 

3979  COM  /Xxyy/  Xx (25) , Yy ( 25 ) 

3982  PRINTER  IS  705 

3985  BEEP 

3988  INPUT  "WANT  TO  PLOT  Uo  vs  Vtu?  (  1  -Y,  0-N) " ,  Iuo 

3991  IF  Iuo-0  THEN 

3994  BEEP 

3997  INPUT  '"SELECT  (0-h/hOX  same  tube.  1  -h(HF)  /h(sm)  "  ,  Irt 
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4000  BEEP  . 

4003  INPUT  "SELECT  h/h  RATIO  ( 1-FILE.O-CQMPUTED) " . Ihrat 

4006  IF  lhrat=0  THEN 

4009  BEEP 

4012  INPUT  "WHICH  Tsat  ( 1-6.7,0  — 2.2)  "  ,Isat 
40 15  END  IF 

4013  Xmin-0 
4021  Xmax-10 
4024  Xstep-2 

4027  IF  Irt-0  THEN 

4030  Ynun-0 

4033  Ymax-1.4 

4036  Ystep-.2 

4039  ELSE 

4042  Ymin-0 

4045  Ymax-15 

4048  Ystep-5 

4051  END   IF 

4054  ELSE 

4057  Qpo-2 

4060  Ynun-0 

4063  Ymax-12 

4066  Ystep-3 

4069  Xmn-fl 

40  72  Xmax-4 

4075  Xstep-1 

4079  END   IF 

4081  IF    Ihrat-1    THEN 

4084  Ynun-0 

4087  Ymax-18 

4090  Ystep-3 

4093  Xmin-0 

4096  Xmax-9 

4099  Xstep-2 

4102  END    IF 

4105  BEEP 

4108  PRINT   "IN:SP1:IP  2300.2200,3300,6300:" 

4111  PRINT    *SC   0, 100.0, 100: TL   2.0:" 

4114  Sfx-100/(Xmax-Xnun) 

4117  Sfy-100/(Ymax-YminJ 

4120  PRINT    "PU   0,0   PD" 

4123  FOR   Xa-Xmin    TO   Xmax   STEP    Xstep 

4126  X-(Xa-Xnun)*Sfx 

4129  PRINT   "PA":X,\0:    XT;" 

4132  NEXT   Xa 

4135  PRINT    "PA    100,0;PU:" 

4138  PRINT    "PU   PA   0,0   PD" 

4141  FOR    Ya-Ynun    TO    Ymax    STEP    Ystep 

4144  Y-(Ya-Ynun)*Sfy 

4147  PRINT    "PA   0,";Y,"YT" 

4150  NEXT   Ya 

4153  PRINT   "PA   0,100    TL   0   2" 

4156  FOR   Xa-Xnun    TO   Xmax    STEP    Xstep 

4159  X-(Xa-Xmin)*Sfx 

4162  PRINT    "PA":X,\100:    XT" 
4165  NEXT  Xa 

4163  PRINT    "PA    100,100   PU   PA    100,0   PD" 
4171  FOR    Ya-Ynun    TO    Ymax   STEP    Ystep 
4174  Y-(Ya-Ymin)*Sfy 

4177  PRINT    "PD   PA    100. " ,Y, "YT" 

4180  NEXT   Ya 
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4183 

4136 

4183 

4132 

4135 

4138 

4201 

4204 

4207 

4210 

4213 

4216 

4213 

4222 

4225 

4228 

4231 

4234 

4237 

4240 

4243 

4246 

4243 

4252 

4255 

4258 

4261 

4264 

4267 

a 

4270 

4273 

4276 

4273 

4232 

4235 

4288 

4231 

4234 

4237 

4300 

4303! 

4306 

4303 

4312! 

4315 

4313 

4321 

4324 

4327 

4330 

4333 

4336 

4333 

4342 

4345 

4348 

4351 

4354 

4357 


PRINT  "PA  100. 100  PIT 
PRINT  "PA  0,-2  SR  1 .5,2" 
FOR  Xa-Xmin  TO  Xmax  STEP  Xstep 
X-(Xa-Xmin)»Sfx 
PRINT  "PA":x,"  .0:" 

IF  Iuo-0  THEN  PRINT  "CP  -2. -1 ;IB" ;Xa: " " 
IF  Iuo-1  THEN  PRINT  "CP  -1 .5. -1 ;LB" ;Xa; 
NEXT  Xa 

PRINT  "PU  PA  0,0" 
FOR  Ya-Ynun  TO  Ymax  STEP  Ystep 
IF  ABS(YaK1.E-5  THEN  Ya-0 
Y-(Ya-Ymin)*Sfy 
PRINT  "PA  0,";Y."" 
IF  Iuo-0  THEN  PRINT  'CP  -4, 
THEN  PRINT  "CP  -3, 


25;LB":Ya:"' 
25;LB";Ya:*' 


'Oil  Percent" 

THEN 


IF  Iuo-1 

NEXT  Ya 

X label S- 

IF  Iuo-0 

IF  Irt-0  THEN 

YlabelS«"h/h07." 

ELSE 

Ylabel$-"h/hsmooth" 

END  IF 

PRINT  "SR  t.S.2:PU  PA  50,-10  CP" : -LEN(XlabelS) /2: "0:LB" ;XlabelS: "" 

PRINT  "PA  -11,50  CP  0,";-LEN(Ylabel$)/2*5/6;"DI  0 , 1 :LB" ; YlaoeiS; "" 

PRINT  "CP  0,0" 

ELSE 

PRINT 

PRINT 


■spo:sp2" 

'SR  1 .2.2. 4;PU  PA 


•8,35:DI  0.1:LBU:PR  1,0.5:LBo:PR  -1.0.5;LB  CkU/m 


PRINT  "PR  -1,0.5:SR  1 . ! .5:LB2:SR  1.5.2:PR  .5..5:LB.:PR  .S.0:LBK)" 

PRINT  "PA  42.-10;DI  1,0:LBV:PR  .4,-1;LBu>;PR  1  ,  .5;LB(m/sJ " 

PRINT  "SP0:SP1" 

END  IF 

Ipn-0 

BEEP 

INPUT  "UANT  TO  PLOT  DATA  FROM  A  FILE  ( 1-Y.O-N)?" ,0kp 

Icn-0 

IF  Okp-1  THEN 

BEEP 

INPUT  "ENTER  THE  NAME  OF  THE  DATA  FILE" ,0  f i ieS 

IF  Iuo-0  THEN 

BEEP 

INPUT  "SELECT  (0-LINEAR.  1-L0G(X,  Y) "  .Hog 

END  IF 

ASSIGN  *File  TO  D_fileS 

BEEP 

INPUT  "ENTER  THE  BEGINNING  RUN  NUMBER", fid 

BEEP 

INPUT  "ENTER  THE  NUMBER  OF  X-Y  PAIRS  STORED" .Npairs 

IF  Iuo-0  AND  Ihrat-0  THEN 

BEEP 

INPUT  "ENTER  DESIRED  HEAT  FLUX'\Q 

END  IF 

BEEP 

PRINTER  IS  1 

PRINT  USING  *4X. ""Select  a  symbol:""" 

PRINT  USING  "4X,""1   Star    2   Plus  sign"*" 

PRINT  USING  "4X,""3   Circle  4   Square""" 

PRINT  USING  "4X,""5   RomDus""" 
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4360  PRINT  USING  "4X.""6   Right-side-up  tr langie"  ■", 

4363  PRINT  USING  "4X,""7   Up-3ide-down  triangle 

4366  INPUT  Sym 

4363  PRINTER  IS  705 

4372  PRINT  "PU  DI" 

4375  IF  Sym-1  THEN  PRINT  "SM*" 

4378  IF  Sym-2  THEN  PRINT  "SM+" 

4331  IF  Sym-3  THEN  PRINT  *SMo" 

4384  Nn-4 

4387  IF  Ilog-1  THEN  Nn-! 

4390  IF  Md>1  THEN 

4393  FOR  1-1  TO  (ttd-1) 

4396  ENTER  »FiU;Xa,Ya 

4399  NEXT  I 

4402  END  IF 

4405  IF  Ihrat-0  THEN 

4408  Ql-Q 

4411  IF  Ilog-1  THEN  Q-LQG(Q) 

4414  END  IF 

4417  FOR  I-!  TO  Npairs 

4420  IF  Iuo-0  AND  Ihrat-0  THEN 

4423  ENTER  3Fi le;Xa,B(*) 

4426  Ya-B(O) 

4429  FOR  K-!  TO  Nn 

4432  Ya-Ya+B(K)*Q*K 

4435  NEXT  K 

4438  END  IF 

4441  IF  Iuo-1  OR  Ihrat-1  THEN 

4444  ENTER  SFiietXa.Ya 

4447  IF  Iuo-t  THEN  Ya-Ya/1000 

4450  END  IF 

4453  IF  Iuo-0  AND  Ihrat-0  THEN 

4456  IF  Ilog-1  THEN  Ya-EXP(Ya) 

4459  IF  Ilog-0  THEN  Ya-Ql/Ya 

4462  IF  Irt-0  THEN 
4465  IF  Xa-0  THEN 

4463  Yo-Ya 
4471  Ya-I 
4474  ELSE 
4477  Ya-Ya/Yo 
4480  END  IF 
4483  ELSE 

4486  Hsm-FNHsmooth(Q.Xa.Isat) 

4489  Ya-Ya/Hsm 

4492  END  IF 

4495  END  IF 

4498  Xx(I-1)-Xa 

4501  Yy(I-1)-Ya 

4504  X-(Xa-Xmin)*Sfx 

4507  Y-(Ya-Ymin)*Sfy 

4510  IF  Sym>3  THEN  PRINT  USM" 

4513  IF  Sym<4  THEN  PRINT  "SR  1.4,2.4" 

4516  PRINT  "PA-.X.Y."* 

4519  IF   Sym>3    THEN   PRINT    *SR    1.2.1.6" 

4522  IF   Sym-4   THEN  PRINT    *UC2. 4. 99. 0 , -8.-4.0.0,3.4,0, ; * 

4525  IF   Sym-5    THEN   PRINT    "UC3.0 .99. -3. -6,-3.6.3,6.3, -6: 

4528  IF   Sym-6    THEN  PRINT    "UC0. 5. 3, 99. 3, -8. -6.0,3,3: " 

4531  IF  Sym-7   THEN   PRINT    "UC0 , -5.3,99, -3,3,6,0 , -3. -8; " 

4534  NEXT    I 

4537  BEEP 

4540  ASSIGN  3File   TO   «■ 
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4543  END  IF 

4546  PRINT  "PU  SM" 

4549  BEEP 

4552  INPUT  "WANT  TO  PLOT  A  POLYNOMIAL  < l-Y,0-N)?",Okp 
4555  IF  Qlcp-1  THEN 

4553  BEEP 

4561  PRINTER  IS  1 

4564  PRINT  USING  "4X. ""Select  line  type:""" 

4567  PRINT  USING  "SX.""0     Solid  line 

4570  PRINT  USING  "6X.""1     Dashed""" 

4573  PRINT  USING  "6X,""2.,,S  Longer  line  -  dash*"" 

4576  INPUT  Ipn 

4579  PRINTER  IS  705 

4582  BEEP 

4585  INPUT  "SELECT  ( 0-LINEAR, 1 -LQG(X, Y) ) " , Ilog 

4588  Iprnt-1 

4591  CALL  Poly 

4594  IF  Iuo-1  THEN 

4597  BEEP 

4600  INPUT  "DESIRE  TO  SET  X  Lower  and  Upper  Limit  ( 1-Y ,0-N)?" , Ixlim 

4603  IF  Ixlim-0  THEN 

4606  Xll-0 

4609  Xul-7 

4612  END  IF 

4615  IF  Ixlim-1  THEM 

4618  BEEP 

462!  INPUT  "ENTER  X  Lower  Limit", Xll 

4624  BEEP 

4627  INPUT  "ENTER  X  Upper  Limit", Xul 

4630  END  IF 

4633  END  IF 

4636  FOR  Xa-Xll  TO  Xul  STEP  Xstep/25 
4639  Icn-Icn+1 

4642  Ya-FNPoly(Xa) 

4645  IF  Iuo-1  THEN  Ya-Ya/1000 

4648  Y»(Ya-Ymin)*Sfy 

4651  X-(Xa-Xnun)*Sfx 

4654  IF  Y<0  THEN  Y-0 

4657  IF  Y>100  THEN  GOTO  4637 

4660  Pu-0 

4663  IF  Ipn-!  THEN  Idf-Icn  MOD  2 

4666  IF  Ipn-2  THEN  Idf-Icn  MOD  4 

4669  IF  Ipn-3  THEN  Idf-Icn  MOD  8 

4672  IF  Ipn-4  THEN  Idf-Icn  MOD  16 

4675  IF  Ipn-5  THEN  Idf-Icn  MOD  32 

4678  IF  Idf-1  THEN  Pu-1 

4631  IF  Pu-0  THEN  PRINT  "PA"  .X, Y, "PD" 

4634  IF  Pu-1  THEN  PRINT  "PA"  .X, Y. "PU" 

4637  NEXT  Xa 
4690  PRINT  "PU" 
4693  GOTO  4285 
4696  END  IF 
4699  BEEP 

4702  INPUT  "WANT  TO  QUIT  ( 1-Y,Q-N)?" ,Iquit 

4705  IF  Iquit-1  THEN  4711 

4708  GOTO  4285 

4711  PRINT  "PU  SPO" 

4714  SUBEND 

4717  SUB  Stats 

4720  PRINTER  IS  701 
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4723 

J-0 

4726 

K«0 

4729 

BEEP 

4732 

IF  Iplot-1  THEN  ASSIGN  SFile  TO  P_file$ 

4735 

BEEP 

4738 

INPUT  "LAST  RUN  No?( 0-QUIT) " ,Nn 

4741 

IF  Nn-0  THEN  4849 

4744 

Nn-Nn-J 

4747 

Sx-0 

4750 

Sy-0 

4753 

Sz-0 

4756 

Sxs-Q 

4759 

Sys-0 

4762 

Szs-0 

4765 

FOR  1-1  TO  Nn 

4768 

J-J+1 

4771 

ENTER  SFiletQ.T 

4774 

H-Q/T 

4777 

Sx-Sx+Q 

4780 

Sxs-Sxs+Q"2 

4783 

Sy-Sy+T 

4786 

Sys-Sys+T*2 

4789 

Sz-Sz+H 

4792 

Szs-Szs+H-2 

4795 

NEXT  I 

4798 

Qave-Sx/Nn 

4801 

Tave»Sy/Nn 

4804 

Have"Sz/Nn 

4807 

Sdevq-SQR(ABS( (Nn*Sxs-Sx "2) /(Nn+(Nn-1 ) ) ) ) 

4810 

Sdevt-SOR(ABS( (Nn*Sys-Sy *2) /(Nn*( Nn-1 )))) 

4813 

Sdevh-SQR(ABS( (Nn*Szs-Sz*2) /(Nn+(Nn-1 ) ) ) ) 

4816 

Sh-100*Sdevh/Have 

4819 

Sq-!00*Sdevq/Oave 

4822 

St-100*Sdevt/Tave 

4825 

IF  K-1  THEN  4843 

4828 

PRINT 

4831 

PRINT  USING  "MIX,  ""DATA  FILE:  ""  ,  I4A"  ;Fi  leS 

4834 

PRINT 

4837 

PRINT  USING  "1 IX, ""RUN  Htube     SdevH   Qdp 

SdevQ   Thetab  SdevT' 


4840  K-1 

4843  PRINT  USING  "  1 1X,DD,2(2X,D.3DE. IX.3D.2D) .2X.DD.3D, 1X.3D.2D" ; J,Have,Sh,Qave 

,Sq, Tave.St 

4846  GOTO  4735 

4849  ASSIGN  SFilel  TO  * 

4852  PRINTER  IS  1 

4855  SUBEND 

4858  SUB  Coet 

4861  COM  /Cply/  A( 10, 1 0) ,C( 10 ) ,B(4) .N, Iprnt .Opo, Hog, If n, Ijoin,N join 

4864  BEEP 

4867  INPUT  "GIVE  A  NAME  FOR  CROSS-PLOT  FILE'.CpfS 

4870  CREATE  BOAT  Cp+"S.2 

4873  ASSIGN  SFile  TO  CpfS 

4876  BEEP 

4879  INPUT  "SELECT  (  0-LINEAR.  1  -L0G(X.  Y)  )  "  ,  Hog 

4882  BEEP 

4885  INPUT  "ENTER  OIL  PERCENT  (-1 -STOP) ", Bop 

4888  IF  BoP<0  THEN  4900 

4891  CALL  Poly 

4894  OUTPUT  wFile;Bop.B<*> 

4897  GOTO  4882 
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4900 

4903 

4906 

4909 

4912 

4915! 

4913 

4921 

4924 

4927 

4930 

4933 

4936 

4939 

4942 

4945 

4948 

4951 

4954 

4957 

4960 

4963 

4966 

4969 

4972 

4975 

4978 

4981 

4984 

4987 

4990 

4993 

4996 

4999 

5002 

5005 

5008 

5011 

5014 

5017 

5020 

5023 

5026 

5029 

5032! 

5035 

5038 

5041 

5044 

5047 

5050 

5053 

5056 

5059 

5062 

5065 

5063 

5071  ! 

5074! 

5077 


ASSIGN  'sFile  TO  * 

SUBEND 

SUB  WilsontCf ,Ci) 

COM  /Mil/  D2.Di,Do,L,Lu.:<cu 

DIM  Emt(l2) 

WLISON  PLOT  SUBROUTINE  DETERMINE  CF  AND  CI 

BEEP 

INPUT  "ENTER  DATA  FILE  NAME\File5 

BEEP 

PRINTER  IS  1 

PRINT  USING  "4X. ""Select  option: 


PRINT 
PRINT 
PRINT 
INPUT 

PRINTER  IS 
IF  Icfix-0 
IF  lcfix>0 
IF  Icfix-1 
IF  Icfix-2 
PRINTER  IS 
INPUT  "Want 
IF  Iccoef -1 


USING  "4X.""  0  Vary  Cf  and  Ci""" 
USING  "4X.""  1  Fix  Cf  Vary  Ci""" 
USING  "4X.""  2  Vary  Cf  Fix  Ci"" 
"ENTER  OPTION", let ix 

701 

THEN  4960 

THEN  BEEP 

THEN  INPUT  "ENTER  Cf.Csf 

THEN  INPUT  "ENTER  CI",Ci 

1 
To  Vary  Coef f ?( 1 -Y.O-N) " , Iccoef 
THEN  INPUT  "ENTER  COEFF",R 


PRINTER  IS  701 

IF  Icfix-0  OR  Icfix-2  THEN  Cfa-.Q04 

IF  Icfix-1  THEN  Cfa-Csf 

Cia-Ci 

Xn-.3 

Fr-.3 

Jj-0 

Rr-3. 

IF  Iccoef-!  THEN  Rr-R 

PRINTER  IS  1 

PRINT  Do.Di,Kcu 

ASSIGN  *File  TO  FileS 

ENTER  *File:Mrun. Dates. Ldtcl ,Ldtc2.Itt 

Ru-Do*L0G ( Do/D i ) / ( 2*Kcu ) 

Sx-0 

Sy-0 

Sxy-0 

Sx2-0 

Sy2»0 

FOR  1-1  TO  Nrun 

ENTER    2>File:Bop.ToldS.Emf  (*)  .Fms 

CONVERT  EMF'S  TO  TEMPERATURE 

FOR  J-1  TO  5 

T(J)-FNTvsv(Emf (J)  ) 

NEXT  J 

Tsat-(T(1)+T(2) )*.5 

Tava-T(5) 

Grad-37.3853+. 104388*Tava 

Tdrop-Emf ( 7)*1 .E+6/( 1 0.*Grad) 

Tavgc-T(5) -Tdrop*.5 

IF   ABS(Tavg-TavgcJ>.01    THEN 

Tavg-( Tavg+Tavgc) *.5 

GOTO   5050 

END    IF 

Compute   properties   of    water 
Kw-FNKuj(Tavg) 
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5080  rtuwa-FNMuu(TavgJ 

5083  Cpw-FNCpw(TavgJ 

5086  Prtii-FNPru*  Tava) 

5089  Rhou-FNRhow(Tavg) 

5092! 

5095!  Compute  properties  of  Freon-114 

5098  Lmtd-Tdrop/L0G((T(5)-Tsat)/(T(5)-Tdrop-TsatJ ) 

5101  IF  Jj-0  THEN 

5104  Tu-Tsat+Fr*Lmtd 

5107  Thetab-Tw-Tsat 

5110  Jj-1 

5113  END  IF 

5116  Tf-(Tui+Tsat)*.S 

5119  Rho-FNRho<Tf ) 

5122  Mu-FNMu(Tf) 

5125  K»FNK(Tf) 

5128  Cp-FNCp(Tf) 

5131  Beta»FNBeta( Tf) 

5134  Hfg-FNHfg(Tsat) 

5137  Ni-Mu/Rho 

5140  Alpha-K/(Rho*Cp) 

5143  Pr-Ni/Alpha 

5146! 

5149!  Analysis 

5152!  COMPUTE  MOOT 

5155  A-PI*(Do'2-Di*2)/4 

5158  P-PI*0o 

5161  ltdot-3.9657E-3-t-Fms*<3.61955E-3-Fms*(8.32006E-6-Fms*(  1  .23688E-7-Fms*4. 31897 

E-10))) 

5164  Q-Mdot*Cpu»*Tdrop 

5167!  COMPUTE  NATURAL -CONNECTIVE  HEAT-TRANSFER  COEFFICIENT 

5170!  FOR  UNENHANCED  END(S) 

5173  Hbar-190 

5176  Fe-(Hbar+P/(Kcu*A) ) '.5*Lu 

5179  Tanh-FNTanh(Fe) 

5182  Theta-Thetab*Tanh/Fe 

5185  Xx-(9.81*8eta*Thetab+0o*3*Tanh/(Fe*Ni+Alpha) ) ".166667 

5188  Yy-( 1 +(.559/Pr)' (9/1 S))' (8/27) 

5191  Hbarc-K/Do*< .6+.387*Xx/Yy ) '2 

5194  IF  ABS((Hbar-Hbarc)/Hbar)>.001  THEN 

5197  Hbar-(Hbar+Hbarc)*.S 

5200  GOTO  5176 

5203  END  IF 

5206! 

5209!  COMPUTE  HEAT  LOSS  RATE  THROUGH  UNENHANCED  ENDS 

5212  Q1-(Hbar+P*Kcu+A)' .5*Thetab*Tanh 

5215  Qc-Q-2*Q1 

5218  As=PI*D2*L 

5221 !  COMPUTE  ACTUAL  HEAT  FLUX 

5224  Qdp-Qc/As 

5227  IF  Icfix-0  OR  Icfix>1  THEN  Csf-!/Cf * ( 1 ./Rr) 

5230  Thetab-Csf/Cp*Hfg*(Odp/(Mu*HfgJ*( .01 4/(9.81 *Rho) ) * .5) * ( 1 /Rr ) *Pr* 1  .  7 

5233  Ho-Qdp/Thetab 

5236  Omega-Ho/Cf 

5239  Uo-Q/(PI*Do*L*Lmtd) 

5242  V.u-Mdot/(Rhoui*PI*D  i*2/4) 

5245  Reui-Rhour*Vw*Di/Muua 

5248  Twi-Tw+Q*Ru/(PI*Do+L) 

5251  Gama-Kui/Di+Rew"  .3*Pru"  (1/3.  )*(Muua/FNMuwl  Twi  )  )  "  .14 

5254!    PRINTER    IS    1 

5257  Yu-  (  1  .  /Uo-Ru )  Omega 
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5260  Xw-Omeaa*Do/(Gama*Di) 

52S3  Sx-Sx+Xw 

52S6  Sy-Sy+Yu 

5263  Sxy-Sxy+Yu»*Xu» 

5272  Sx2-Sx2+Xu,*Xw  "  ■   ♦ 
5275  Sy2-Sy2+Yu>*Yu 

5273  NEXT  I 

5281  ASSIGN  *File  TO  * 

5284  r1-(Sx*Sy-Nrun*Sxy)/(Sx*Sx-Nrun+Sx2) 

5287  C-(Sy-Sx*M)/Nrun 

5230  IF  Icfix-0  OR  Icfix-3  OR  Icfix-4  THEN 

5233  Cic-t/M 

5236  Cfc-1/C 

5233  END  IF 

5302  IF  Icfix-1  THEN 

5305  Cic-l/N 

5308  Cfc-Cf 

5311  END  IF 

5314  IF   Icfix-2   THEN 

5317  Cic-Ci 

5320  Cfc-1/C 

5323  END  IF 

5326  IF  ABSUCi-Cic)/Cic)>.001    OR  ABS( (Cf-Cfc) /Cfc) > .001    THEN 

5323  Ci-(Ci+Cic)*.S 

5332  Cf-(Cf+CfcJ*.S 

5335  PRINTER  IS  1 

5338  PRINT  USING  "2X,""  Csf  ■  •■.f1Z.3DE.2X.""  Ci  -  ""  ,f1Z.3DE":Csf  ,Ci 

5341  PRINTER  IS  701 

5344  GOTO  5002 

5347  END  IF 

5350  PRINT 

5353  PRINTER  IS  701 

5356  PRINT  USING  "23X.""   Cf  Ci""" 

5353  PRINT  USING  "3X.  "ASSUMED      "  ,MZ:3DE .3X .MZ.3DE" tCfa.Cia 

5362  PRINT  USING  *3X. -"CALCULATED    "" ,MZ.3DE.3X,J1Z.3DE" ;Csf .Ci 
5365  PRINT 

5363  Sum2-Sy2-2*M*Sxy-2*C*Sy+W*2*Sx2+-2J*r1*C*Sx+Nrun*C*2 
5371  PRINT  USING  "10X.""Sum  of  Squares  -  " " . Z.3DE" :Sum2 
5374  PRINT  USING  *10X. ""Coefficient  -  "".D.DDD":Rr 
5377  SUBEND 

5380  DEF  FNMuu(T) 

5383  A-247.3/(T>133.15) 

5386  Mu-2.4E-5+10*A 

5383  RETURN  Mu 

5332  FNEND 

5335  DEF  FNCpw(T) 

5338  Cpu-4 .211 20858-T* ( 2 . 26826E-3- T* ( 4 . 4236 1 E-S+2 . 7 1 428E-7*T ) ) 

5401  RETURN  Cpw*1000 

5404  FNEND 

5407  DEF  FNRhow(T) 

5410  Ro-933.52346+T*( .01263-T*<5. 4825 13E-3-T*1 .234147E-5) ) 

5413  RETURN  Ro 

5416  FNEND 

5419  DEF  FNPru(T) 

5422  Pru»-FNCpu(T)*FNMuu(T)/FNKui(T) 
5425  RETURN  Prui 

5423  FNEMD 

5431  DEF  FNKu(T) 

5434  X-(T+273.IS)/273.1S 

5437  Kw--.92247*X*(2.3395-X*( 1 .3007-X*( .52577-. 07344*X) ) ) 
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5440  RETURN  Ku 

5443  FMEMD 

5446  SUB  Plot 

5449  COM  /Cply/  A( 10,10) ,C( 10) ,B<4) .Nop, Iprnt.Opo.Ilog.Ifn.Ijoin.Njoin 

5452  DIM  Bs(3) 
5455  INTEGER  H 

5453  PRINTER  IS  1 
5461  BEEP 

5464  INPUT  "SELECT  HEATING  MODE  (0-ELECTRIC,  1 -UATER) " , Ihm 

5467  Idv-0 

5470  BEEP 

5473  INPUT  "LIKE  DEFAULT  VALUES  FOR  PLOT  ( 1 -Y.0-N)?" , Idv 

5476  Opo-0 

5479  BEEP 

5482  PRINT  USING  "4X. ""Select  Option:""" 

5485  PRINT  USING  "6X,""0   q  versus  delta-T""" 

5488  PRINT  USING  "6X.""1   h  versus  delta-T""" 

5491  PRINT  USING  "6X,""2  h  versus  q""" 

5494  INPUT  Opo 

5497  BEEP 

5500  INPUT  "SELECT  UNITS  (0-SI, 1 -ENGLISH) ", Iun 

5503  PRINTER  IS  705 

5506  IF  Idv<>1  THEN 

5509  BEEP 

5512  INPUT  "ENTER  NUMBER  OF  CYCLES  FOR  X-AXIS", Cx 

5515  BEEP 

5518  INPUT  "ENTER  NUMBER  OF  CYCLES  FOR  Y-AXIS", Cy 

5521  BEEP 

5524  INPUT  "ENTER  MIN  X-VALUE  (MULTIPLE  OF  10)".Xmin 

5527  BEEP 

5530  INPUT  "ENTER  MIN  Y-VALUE  (MULTIPLE  OF  10)".Ymin 

5533  ELSE 

5536  IF  Opo-0    THEN 

5539  IF    Ihm-1    THEN 

5542  Cy-I 

5545  Cx-2 

5548  Xnun-1 

5551  Ymin-10000 

5554  ELSE 

5557  Cy-3 

5560  Cx-3 

5563  Xmin-.1 

5566  Ymin-100 

5569  END    IF 

5572  END   IF 

5575  IF  0po-1    THEN 

5578  IF    Ihm-1    THEN 

5581  Cy-2 

5534  Cx-2 
5587  Xmin-.1 
5590  Ymin-100 
5593  ELSE 
5596  Cy-3 
5599  Cx-2 
5602  Xnun-1000 
5605  Ymin-100 
5608  END   IF 
5611  END    IF 

5614  IF  Opo-2   THEN 

5617  IF    Ihm-1    THEN 

5620  Cy-2 
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5623 

Cx-1 

5626 

Xnun-10000 

5623 

Ynun-100 

5632 

ELSE 

5635 

Cy-3 

5638 

Cx-3 

5641 

Xnun- .  1 

5644 

Ymin-100 

5647 

END  IF 

5650 

END  IF 

5653 

END  IF 

5656 

BEEP 

5659 

PRINT  *IN;SPt;IP  2300,2200.8300.6800 

.  * 

5662 

PRINT  "SC  0, 100.0, 1 00: TL  2,0:" 

5665 

Sfx-100/Cx 

5663 

Sfy-100/Cy 
BEEP 

5671 

5674 

INPUT  'UANT  TO  BY-PASS  CAGE?  (1-Y.0- 

N) 

5677 

IF  Ibyp-1  THEN  6049 

5630 

PRINT  "PU  0.0  PD" 

5683 

Nn-9 

5636 

FOR  1-1  TO  Cx-1 

5689 

Xat-Xmin*1Q*(I-1) 

5692 

IF  I-Cx+1  THEN  Nn-1 

5695 

FOR  J-1  TO  Nn 

5698 

IF  J-1  THEN  PRINT  "TL  2  0" 

5701 

IF  J-2  THEN  PRINT  "TL  1  0' 

5704 

Xa-Xat*J 

5707 

X-LGT(Xa/Xmin)*Sfx 

5710 

PRINT  "PA*;X,\0:  XT;" 

5713 

NEXT  J 

5716 

NEXT  I 

5719 

PRINT  'PA  100,0;PU:" 

5722 

PRINT  "PU  PA  0,0  PD' 

5725 

Nn-9   " 

5728 

FOR  1-1  TO  Cy+1 

5731 

Yat-Ymin*10*(I-1 ) 

5734 

IF  I-Cy-H  THEN  Nn-1 

5737 

FOR  J-1  TO  Nn 

5740 

IF  J-1  THEN  PRINT  "TL  2  0* 

5743 

IF  J-2  THEN  PRINT  "TL  1  0' 

5746 

Ya-Yat*J 

5749 

Y-LGT(Ya/Ynun)*Sfy 

5752 

PRINT  "PA  0,":Y,"YT" 

5755 

NEXT  J 

5758 

NEXT  I 

5761 

PRINT  "PA  0,100  TL  0  2" 

5764 

Nn-9 

5767 

FOR  I-I  TO  Cx+1 

5770 

Xat-Xmin*10"(I-1 ) 

5773 

IF  I-Cx-I  THEN  Nn-1 

5776 

FOR  J-1  TO  Nn 

5779 

IF  J-1  THEN  PRINT  'TL  0  2" 

5782 

IF  J>1  THEN  PRINT  "TL  0  1" 

5785 

Xa-Xat*J 

5788 

X-LGT(Xa/Xnun)*Sfx 

5791 

PRINT  "PA";X,".!00:  XT- 

5794 

NEXT  J 

5797 

NEXT  I 

5800 

PRINT  "PA  100,100  PU  PA  100,0  PD" 

Ibyp 
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5303  Nn-9 

5306  FOR  1-1  TO  Cy  +  1 

5809  Yat- Ymin*10*(I-1) 

5812  IF  E-Cy+1  THEN  Nn- I 

5815  FOR  J-1  TO  Nn 

5818  IF  J- I  THEN  PRINT  "TL  0  2" 

5821  IF  J>1  THEN  PRINT  'TL  0  1" 

5824  Ya-Yat*J 

5827  Y-LGT(Ya/Ymm)*Sfy 

5830  PRINT  "PD  PA  100,", Y, "YT" 

5833  NEXT  J 

5836  NEXT  I 

5839  PRINT  "PA  100,100  PU" 

5842  PRINT  "PA  0,-2  SR  1.5,2" 

5845  Ii-LGT(Xnun) 

5848  FOR  1-1  TO  Cx+1 

5851  Xa-Xmin*10*(I-n 

5854  X-LGT(Xa/Xnun)*Sfx 

5357  PRINT  "PA";X.",0:" 

5860  IF  Ii>-0  THEN  PRINT  "CP  -2, -2:LB10 :PR  -2.2:LB" : Ii :  "" 

5863  IF  li<0  THEN  PRINT  "CP  -2. -2;LB1 0;PR  0 ,2;LB" ; Ii : *" 

5866  Ii-Ii+1 

5869  NEXT  I 

5872  PRINT  "PU  PA  0.0" 

5875  Ii-LGT(Yirun) 

5878  Y10-10 

5881  FOR  1-1  TO  Cy  +  1 

5884  Ya-Ymin*10"(I-1) 

5887  Y-LGT(Ya/Ymin)*Sfy 

5890  PRINT  "PA  0,";Y."" 

5893  PRINT  "CP  -4,-.25:LB1 0 ;PR  -2,2:LB" : Ii ; '" 

5896  li-.Ii+l 

5899  NEXT  I 

5902  BEEP 

5905  INPUT  "WANT  USE  DEFAULT  LABELS  ( 1 -Y.O-N)?- . Idl 

5908  IF  Idloi  THEN 

5911  BEEP 

5914  INPUT  "ENTER  X-LABEL" .XlabelS 

5917  BEEP 

5920  INPUT  "ENTER  Y-LABEL" , YlabelS 

5923  END  IF 

5926  IF  0po<2  THEN 

5929  PRINT  "SR  1 ,2:PU  PA  40.-14:" 

5932  PRINT  "LB(T:PR  -1.6,3  PD  PR  1.2,0  PU:PR  .5, -4;LBu*o;PR  .5,1;" 

5935  PRINT  "LB-T:PR  .5, -1 ;LBsat :PR  .5,1;" 

5938  IF  Iun-0  THEN 

5941  PRINT  "LB)  CIO" 

5944  ELSE 

5947  PRINT  "LB)  (F)" 

5950  END  IF 

5953  END  IF 

5956  IF  Opo-2  THEN 

5959  IF  Iun-0  THEN 

5962  PRINT  "SR  !.S,2;PU  PA  40,-14;LBq  (U/m;SR  1,I.5;PR  0.5,l:LB2:SR  t.5.2;PR 
0.5.-1 ;LB)" 

5965  ELSE 

5963  PRINT  "SR  1.5.2:PU  PA  34,-14:LBq  (Btu/hr:PR  .5,.5;LB.:PR  .5, -.5:" 
5971  PRINT  "LBft;PR  .5,1 ;SR  1,1.5;LB2;SR  1.5.2:PR  .5,-l:LB);" 

5974  END  IF 

5977  END  IF 

5980  IF  Opo-0  THEM 
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5383  IF  Iun-0  THEN 

5386  PRINT  "SR  1.5.2:PU"PA  -12.40:DI  0.l:LBq  (W/m;PR  - 1 ,0.5  :SR  !.!.5:LB2:SR  I 

.5.2:PR  1,.S:L8)" 

5383  ELSE 

5332  PRINT  "SR  f.5.2;PU  PA  -12.32:01  0.1:LBq  (Btu/hr:PR  -.5, .5;LB. :PR  .5.-5;" 

5335  PRINT  "LBft:SR  1,1.S;PR  -1,.5:LB2:PR  1,.5;SR  1.5.2:LB)" 

5338  END  IF 

6001  END  IF 

6004  IF  0po>0  THEN 

6007  IF  Iun-0  THEN 

6010!  PRINT  "SR  1.5,2;PU  PA  -12,38:DI  O.ULBh  (U/m;PR  -1,.5;SR  1,1.5;LB2:SR  1. 

5.2:PR  .5. .5;" 

6013  PRINT  "SR  1.2,2.4:PU  PA  -12.37:01  0,l;LBh;PR  !.0.5:L3o:PR  -1.0.5;LB  C W/m 


6016 

5013 

6022 

6025 

6023 

6031 

6034 

6037 

6040 

6043 

6046 

6043 

6052 

6055 

6058 

6061 

6064 

5067 

6070 

6073 

6076 

6073 

6082 

6085 

6088 

6031 

6034 

6037 

6100! 

6103! 

6106 

6103 

6112 

6115 

6118 

6121 

6124 

6127 

6130 

6133 

6136 

6133 

6142 

6145 

6148 

6151 


PRINT  "PR  -1,.5:SR  1,1.S:LB2:SR  1.5,2:PR  .S..S:LB.:PR  .S,0:LBK)" 

ELSE 

PRINT  "SR  1.5.2:PU  PA  -12.28:01  O.I;LBh  (Btu/hr:PR  -.5. .S:LB. :PR  .5,-5; 

PRINT  "LBftjPR  -1,.5:SR  1,1.5:LB2:SR  1.5.2:PR  .S,.S:LB.:PR  .5,.5:LBF)" 

END  IF 

END  IF 

IF  Idl-0  THEN 

PRINT  "SR  1.S.2:PU  PA  50.-16  CP" : -LEM(XlabelS) /2: "0:LB" :XlabelS: "" 

PRINT  "PA  -14.50  CP  0, " : -LEN( YlabelS) /2*5/6: "01  0, 1 JLB" : YlaoelS; "" 

PRINT  "CP  0,0  01" 

END  IF 

Ipn-0 

XI 1-1 .E+6 

Xul— I.E+6 

Icn-fl 

Ifn-0 

I join-1 

BEEP 

INPUT  "WANT  TO  PLOT  DATA  FROM  A  FILE  ( W,0-N)?",0lr 

IF  Ok-1  THEN 

BEEP 

INPUT  "ENTER  THE  NAME  OF  THE  DATA  FILE" ,D_f i leS 

ASSIGN  SFile  TO  D  fileS 

BEEP 

BEEP 

INPUT 

BEEP 

INPUT 

BEEP 

INPUT  "CONNECT  DATA  UITH  LINE  ( 1 -Y.O-N) ?" , Icl 

BEEP 

PRINTER  IS  1 

PRINT  USING  "4X. ""Select  a  symool:""" 

PRINT  USING  "SX 

PRINT  USING  "6X 

PRINT  USING  "6X 

PRINT  USING  "6X 

PRINT  USING  "6X 

INPUT  Sym 

PRINTER  IS  705 

PRINT  "PU  DI" 

IF  Sym-1  THEN  PRINT  "SM*" 

IF  Sym-2  THEM  PRINT  "SM+" 

IF  Sym-3  THEN  PRINT  "SMo" 

IF  Md>!  THEN 

FOR  1-1  TO  (Md-1 ) 


"ENTER  THE  BEGINNING  RUN  NUMBER", Md 

"ENTER  THE  NUMBER  OF  X-Y  PAIRS  STORED" .Npairs 


"1 

Star     2   Plus  sign""" 

■3 

Circle   4   Square""* 

"S 

Romous""" 

*6 

Right-side-up  triangle"" 

"7 

Up-side-down  triangle""" 
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6154  ENTER  ?File:Ya.Xa 

S157  NEXT  I 

SI 60  END  IF 

6163  FOR  1-1  TO  Noairs 

6166  ENTER  9Fil«:Ya.Xa 

6169  IF  1-1  THEN  Q1-Ya 

6172  IF  I-Npairs  THEN  Q2=Ya 

6175  IF  Qpo-1  THEN  Ya-Ya/Xa 

6178  IF  Opo-2  THEN 

6181  Q-Ya 

6184  Ya-Ya/Xa 

6187  Xa-Q 

6190  END  IF 

6193  IF  Xa<Xll  THEN  Xll-Xa 

6196  IF  Xa>Xul  THEN  Xul-Xa 

6199  IF  Iun-1  THEN 

6202  IF  0po<2  THEN  Xa-Xa*1 .8 

6205  IF  0po>0  THEN  Ya-Ya*.I76t 

6208  IF  Opo-O  THEN  Ya-Ya*.317 

6211  IF  Opo-2  THEN  Xa-Xa*.317 

6214  END  IF 

6217  X-LGT(Xa/Xnun)*Sfx 

6220  Y-LGT(Ya/Ynun)*Sfy 

6223  Kj-0 

6226  CALL  Symb(X. Y.Sym.Icl ,Kj) 

6229  GOTO  6258 

6232  IF  Sym;3  THEN  PRINT  "SM" 

6235  IF  Sym<4  THEN  PRINT  *SR  1.4,2.4" 

6238  IF  Icl-0  THEN 

6241  PRINT  "PA",X,Y,"" 

6244  ELSE 

6247  PRINT  "PA" ,X, Y, "PD" 

6250  END  IF 

6253  IF  Sym>3  THEN  PRINT  "SR  1.2.1.6" 

6256  IF  Sym-4  THEN  PRINT  "UC2 .4,99. 0 . -8. -4, 0  .0 ,8  .4  .0  J " 

6259  IF  Sym-5  THEN  PRINT  "UC3.0 .99. -3. -6. -3,6.3.5,3.-6: " 

6262  IF  Sym-6  THEN  PRINT  "UCO .5.3.99.3. -8. -6. 0 .3,3; " 
6265  IF  Sym-7  THEN  PRINT  "UCO, -5.3,99, -3,3,6,0,-3, -8; " 

6263  NEXT  I 
627!  PRINT  "PU" 
6274  BEEP 

6277  INPUT  "WANT  TO  LABEL?  ( 1 -Y,0-N) " , Ilab 

6280  IF  Ilab-1  THEN 

6283  PRINT  •SP0:SP2B 

6286  BEEP 

6289  IF  Klab-0  THEN 

6292  Xlab-5 

6295  Ylab-85 

6298  INPUT  "ENTER  INITIAL  X,Y  LOCATIONS" ,Xlab. Ylab 

6301  Xtt-Xlab-5 

6304  Ytt-Ylab+8 

6307  PRINT  "SR  1 ,1.5" 

6310  BEEP 

6313  INPUT  "SELECT  LABEL  TYPE  (0-FILE  NAME,  1-TUBE  TYPE) " ,If ntt 

6316  IF  Ifntt-0  THEN 

6319  PRINT  "Sf1:PA",Xtt,Ytt,"LB       7.   Heat  File" 

6322  Ytt-Ytt-3 

6325  PRINT  "SMcPA"  ,Xtt,  Ytt.  "LB       7.        Heat  File" 

6323  ELSE 

6331  PRINT  "SM:PA".Xtt.Ytt."LB       7.        Heat   Tube" 
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boo  / 

6340 

6343 

6346 

5349 

63S2 

6355 

6358 

6361 

6364 

6367 

6370 

6373 

6376 

6379 

6382 

6385 

6388 

6391 

6394 

6397 

6400 

6403 

6406 

6409 

6412 

6415 

6418 

6421 

6424 

6427 

6430 

6433 

6436 

6439 

6442 ! 

6445 

6448 

6451 

6454 

6457 

6460 

6463 

5466 

6469 

6472 

6475 

6478 

6481 

6484 

6487 

6490 

6493 

6496 

6499 

6502 

6505 

6508 

65!  1 

6514 


Ytt-Ytt-3 

PRINT  "PA'.Xtt.Ytt, 

END  IF 


"LB 


Oil   Flux   Type' 


THEN  PRINT  "SM*" 
THEN  PRINT  "SM+" 
THEN  PRINT 


"SMo" 


IF  Sym-I 

IF  Sym-2 

IF  Sym-3 

Klab-1 

END  IF 

KJ-1 

CALL  Symb(Xlab.Ylab.Sym,Icl,Kj) 

PRINT  "SR  1  ,1.S:SM" 

IF  Syn><4  THEN  PRINT  "PR  2,0" 

BEEP 

INPUT  "ENTER  BOP", Bop 

IF  Ifntt-1  THEN 

BEEP 

INPUT    "ENTER   TUBE   TYPE'.TtypeS 

END    IF 

IF  Bop<10  THEN  PRINT  "PR  2.0:L3" :8opt "" 

IF  BoP>9  THEN  PRINT  "PR  ,5,0;LB" ;8oP; "" 

Ihf-0 

IF  Q1>02  THEN  Ihf-1 

IF  Ihf-0  THEN  PRINT  "PR  4,0:LBInc" 

IF  Ihf-1  THEN  PRINT  "PR  4,0;LBDec" 

IF  Ifntt-0  THEN 

PRINT  "PR  2,0;LB";D_file$:"" 

pi  CC 

PRINT  "PR  2,Q;LB";TtypeS:"* 

END  IF 

PRINT  "SP0;SP1:SR  1.5.2" 

Ylab-Ylab-5 

END  IF 

BEEP 

ASSIGN  8Fila  TO  ♦ 

Xll-Xll/1.2 

Xul-Xui*1.2 

GOTO  3040 

END  IF 

PRINT  *PU  SM" 

BEEP 

INPUT  "WANT  TO  PLOT  A  POLYNOMIAL  ( 1 -Y , 0-N) 7" ,Go_on 

IF  Go_on-1  THEN 

BEEP 

PRINTER  IS  1 

PRINT  USING  "4X. 

PRINT  USING  "SX. 

PRINT  USING  "SX, 

PRINT  USING  "SX, 

INPUT  Ipn 

PRINTER  IS  705 

BEEP 

INPUT  "SELECT  (  O-LIN.  1  -L0G(X.  YD  *  ,  Hog 

Iprnt-1 

CALL  Poly 

IF  Ifn-0  THEN 

N join- 1 

BEEP 

INPUT  "ENTER  NUMBER  OF  FILES  TO  JOIN?"  ,Njom 

END  IF 

I join-0 

IF    Ifn<N;oin    THEN    Ijoin-I 


"Select    line    type:""" 

"0  Solid    line 

'"1  m    Dashed""" 

"2,,, 5   Longer    line    -    dash*"" 


128 


G517  IF  Ifn-0  OR  Ijoin-1  THEN 

6520  FOR  Ij-0  TO  3 

6523  8s(Ij)-Bs(Ij)+8(Ij) 

6526  NEXT  Ij 

6523  Ifn-Ifn+1 

6532  END  IF 

6535  IF  Njom-Ifn  THEN 

6538  FOR  Ij-0  TO  3 

6541  8(Ij)-Bs(Ij)/Njoin 

6544  Bs(Ij)-0 

6547  NEXT  Ij 

5550  Ifn-0 

6553  ELSE 

6556  GOTO  6067 

6559  END  IF 

6562  BEEP 

6565  INPUT  "ENTER  Y  LOWER  AND  UPPER  LIMITS" ,Yll,Yul 

6568  FOR  Xx-0  TO  Cx  STEP  Cx/200 

6571  Xa-Xnun+10'Xx 

6574  IF  Xa<Xll  OR  Xa>Xui  THEN  6655 

6577  Icn-Icn+1 

6580  Pu-0 

6583  IF  Ipn-1  THEN  Idf-Icn  MOD  2 

6586  IF  Ipn-2  THEN  Idf-Icn  HOD  4 

6589  IF  Ipn-3  THEN  Idf-Icn  MOD  8 

6592  IF  Ipn-4  THEN  Idf-Icn  MOD  16 

6595  IF  Ipn-5  THEN  Idf-Icn  MOD  28 

6598  IF  Idf-1  THEN  Pu-1 

6601  IF  Opo-0  THEN  Ya-FNPolyCXa") 

6604  IF  Opo-2  AND  Ilog-0  THEN  Ya-Xa/FNPoly ( Xa) 

6607  IF  Opo-2  AND  Ilog-1  THEN  Ya-FNPoly ( Xa) 

6610  IF  Opo-1  THEN  Ya-FNPoiy (Xa) 

6613  IF  Ya<Ymin  THEN  6655 

6616  IF  Ya<Yll  OR  Ya>Yul  THEN  6655 

6619  IF  Iun-1  THEN 

6622  IF  0po<2  THEN  Xa-Xa*1 .3 

6625  IF  0po>0  THEN  Ya-Ya*.1761 

6628  IF  Opo-0  THEN  Ya-Ya+.317 

6631  IF  Opo-2  THEN  Xa«Xa*.317 

6634  END  IF 

6637  Y-LGT(Ya/Ymin)*Sfy 

6640  X-LGT(Xa/Xmin)*Sfx 

6643  IF  Y<0  THEN  Y-0 

6646  IF  Y>100  THEN  GOTO  6655 

6649  IF  Pu-0  THEN  PRINT  "PA" ,X. Y, "PD" 

6652  IF  Pu-1  THEN  PRINT  'PA* ,X, Y, "PU" 

6655  NEXT  Xx 

6658  PRINT  "PU" 

6661  GOTO  6067 

6664  END  IF 

6667  BEEP 

6670  INPUT  "WANT  TO  PLOT  REILLY'S  SMOOTH-TUBE  DATA?  ( 1-Y ,0-N) " , Irly 

6673  IF  Opo-0  OR  Opo-1  THEN 

6676  Xll-3 

6679  Xui-20 

6682  END  IF 

6685  IF  Opo-2  THEN 

6688  Xll-10000 

6691  Xul-100000 

6694  END  IF 
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6637  IF  Irly-1  THEM 

6700  Yl  1-20 

6703  YuWO 

6706  BEEP 

6709  INPUT  "ENTER  LQUER  AND  UPPER  Y-LIMITS  FOR  PLOTTING" . Yll ,Yul 

6712  FOR  Xx-0  TO  Cx  STEP  Cx/200 
6715  Xa-Xnun*10'Xx 

6713  IF  Xa<Xll  OR  Xa>Xul  THEN  6757 
6721  Xl-LOG(Xa) 

6724  IF  Opo-0  THEN  Yl— 2.5402837E-1+X1*(4.9720151-X1*2.S134787E-1 ) 

6727  IF  Opo-1  THEN  Yl— 2. S402337E-1 +X1*(3. 9720151 -Xi*2.S134787E-1 ) 

6730  IF  Opo-2  THEN  Yl— 3. 7073801E- 1 +X1*<3. 7259190E-1 -Xl*6.3826842E-3> 

6733  Ya-EXP(Yl) 

6736  Y-LGT(Ya/Ynun)*Sfy 

6739  X-LGT(Xa/Xmin)*Sfx 

6742  Ipu-0 

6745  IF  Y<Y11  THEN  Ipu-1 

6748  IF  Y>Yul  THEN  GOTO  6757 

6751  IF  Ipu-0  THEN  PRINT  "PA" ,X, Y, "PD" 

6754  IF  Ipu-1  THEN  PRINT  "PA" ,X. Y, "PU" 

6757  NEXT  Xx 

6760  PRINT  "PU" 

6763  END  IF 

6766  BEEP 

6769  INPUT  "WANT  TO  PLOT  ROHSENOW  CORRELATION?  ( 1  -Y,0-N) ■ , Irohs 

6772  IF  Irohs-1  THEN 

6775  Yl  1- 15 

6778  Yui-80 

6781  BEEP 

6734  INPUT  "ENTER  Tsat  (Deg  O'.Tsat 
6787  Csf-.0040 

6790  BEEP 

6793  INPUT  "ENTER  Csf  (DEF-0.004) * .Csf 

6796  Tf-Tsat+2 

6799  FOR  Xx-0  TO  Cx  STEP  Cx/200 

6802  Xa-Xmin+10'Xx 

6805  IF  Xa<Xll  OR  Xa>Xul  THEN  6889 

6308  Xl-LOG(Xa) 

6811  IF  0po<2  THEN  Tf-Tsat+Xa/2 

6814  Rho-FNRho(Tf ) 

6817  K-FNK(Tf) 

6820  Mu-FNMu('Tf) 

6823  Cp-FNCp(Tf) 

6826  Hfg-FNHfg(Tsat) 

6329  Ni-Mu/Rho 

6832  Pr-Cp*Mu/K 

6335  Omega-Cst+Hfg/Cp*( C. 014/(9.81 *Rho) )* .5/ (Wu*Hfg) ) * ( 1 ./3)*Pr  '  1 .7 

6833  IF  Opo-0  THEN  Ya-< Xa/OmegaJ "3 
6841  IF  Opo-I  THEN  Ya-(Xa/Omega J *3/Xa 
6844  IF  Qpo-2  THEN  Ya-Xa' (2./3) /Omega 
6847  IF  Opo-2  THEN 

6850  Tfc-Tsat+Xa/Ya*.5 

6353  IF  ABS(Tf-Tfc)>.01  THEN 

6856  Tf-(Tf+TfcJ*.5 

6859  GOTO  6314 

6362  END  IF 

6365  END  IF 

6863  Y-LGT(Ya/Ynun)+Sfy 

6871  X-LGT(Xa/Xmin)*Sfx 

6874  Ipu-0 
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7060  PRINT 

7063  PRINT  USING  "10X.""  Fms      Tin      Tev      vw"2     Tdrop" "" 

7066  IF  K-fl  THEN 

7069  PRINTER  IS  701 

7072  PRINT 

7075  PRINT  "  Month,  date  and  time:  ":Date£ 

7073  IF  Itt-0  THEN  PRINT  USING  "1  OX. ""Tube  Type:  Wieland  Smooth 

7081  IF  Itt-1  THEM  PRINT  USING  "1  OX . ""Tube  Type:  High  Flux""" 

7084  IF  Itt-2  THEN  PRINT  USING  "10X,""Tube  Type:  Turbo-B""" 

7087  PRINT 

7090  PRINT  USING  "10X,""  Fms      Tin      Tev      Vu*2      Tdrop""" 

7093  PRINTER  IS  1 

7096  K-! 

7099  END  IF 

7102  BEEP 

71  OS  INPUT  'ENTER  FLOWMETER  READING". Fms 

7108  OUTPUT  709; "AR  AF0  AL4  VR1 " 

7111  FOR  L-0  TO  4 

7114  OUTPUT  709; "AS  SA" 

7117  IF  L>0  AND  L<4  THEN  7141 

7120  S-0 

7123  FOR  1-0  TO  9 

7126  ENTER  709:E 

7129  S-S+E 

7132  NEXT  I 

713S  IF  L-0  THEN  Emt (0) -ABS(S/10) 

7138  IF  L-4  THEN  Emt ( 1 ) -ABS(S/1 0) 
7141  NEXT  L 

7144  OUTPUT  709: "AR  AF20  AL20  VRt" 

7147  OUTPUT  709:"AS  SA" 

7150  Etp-0 

7153  FOR  1-0  TO  9 

7156  ENTER  709:Et 

7159  Etp-Etp+Et 

7162  NEXT  I 

7165  Etp-Etp/10 

7168  Tin-FNTvsv(Emf (1  )) 

7171  Tev-FNTvsv(Emf (0) ) 

7174  Grad-37.38S3+.!04388*Tin 

7177  Mdot-3.9657E-3+Fms*<3.61955E-3-Fms*<8.S2006E-6-Fms*( 1 .23638E-7-Fms*4.31897 

E-10)  )) 

7180  Vw«Mdot/(1000*PI*Di"2)*4 

7183  Tdrop-Etp*1 .E+6/< 1 0*Grad) 

7186  PRINT  USING    "10X.3(DD.DD.4X) .  IX. Z.DD.4X.MZ.4D" ;Fms. Tin, Tev. vuT2. Tdrop 

7139  BEEP 

7192  INPUT  "WANT  TO  ACCEPT  THIS  DATA  SET?  ( 1 -Y,0-N) *  ,0k 

7195  J-J>1 

7198  IF  Ok-O  THEN 

7201  J-J-1 

7204  GOTO  7102 

7207  ELSE 

7210  OUTPUT  *File:Fms.Emt(*),Etp 

7213  PRINTER  IS  701 

7216  PRINT  USING  "10X,3(DD.DD,4X) , 1X,Z.DD,4X.MZ.4D" ;Fms, Tin. Tev, Vw*2, TdroD 

7219  PRINTER  IS  1 

7222  BEEP 

7225  INPUT  "WILL  THERE  BE  ANOTHER  DATA  SET7  ( 1-Y.0-N) " ,Go  on 

7223  IF  Go  on-1  THEN  7102 
7231  END  IF 

7234  ASSIGN  SFile  TO  * 
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7.227  PRINTER    IS    701 

7240  PRINT 

7243  PRINT  USING  "1  OX. ""NOTE:  ",11,"    data  sets  are  stored  in  file  "". 15A" ; J,F 

iieS 

7246  PRINTER  IS  1 

7249  SUBEND 

72S2  SUB  Uoprt 

72S5  PRINTER  IS  1 

72S8  BEEP 

7261  INPUT  'Enter  Uo  File  Name" .FileS 

7264  BEEP 

7267  INPUT  "Number  of  Data  Runs'.Nrun 

7270  INPUT  "Do  You  Uant  a  Plot  Fiie?( 1 -Y.0-N) ■ .Iplot 

7273  BEEP 

7276  IF  Iplot-1  THEN 

7279  INPUT  "Give  Plot  File  Name'.P  fileS 

7282  CREATE  BOAT  P_file$,4 

728S  ASSIGN  «Plot  TO  P_fileS 

7288  END  IF 

7291  PRINTER  IS  701 

7294  PRINT 

7297  PRINT 

7300  PRINT  USING  "10X.""   Water  Vel        Uo""" 

7303  ASSIGN  SFiie  TO  FileS 

7306  IF  Iplot-1  THEN  ASSIGN  »Filel  TO  P_fileS 

7309  FOR  1-1  TO  Nrun 

7312  ENTER  »File:Vu.Uo 

731S  IF  Iplot-1  THEN  OUTPUT  SFilel tVu.Uo 

7318  PRINT  USING  M5X,D.DD.SX.MZ.3DE";Vw,Uo 

7321  NEXT  I 

7324  ASSIGN  »File  TO  * 

7327  ASSIGN  *File1  TO  * 

7330  PRINT  USING  "10X,""NOTE:  ".11,"    data  sets  are  stored  in  file  '", ISA" ;Nru 

n. FileS 

7333  IF  Iplot-1  THEN 

7336  PRINT  USING  "10X.*"NOTE:  "  ,11,"    X-Y  Pairs  are  stored  in  file  "",15A";Nru 

n.P_f iieS 

7339  END  IF 

7342  PRINTER  IS  1 

734S  SUBEND 

7348  SUB  Select 

7351  COM  /Idp/  Idp 

73S4  BEEP 

7357-  PRINTER  IS  1 

7360  PRINT  USING  "4X. ""Select  option:""" 

7363  PRINT  USING  "SX.""  0  Talcina  data  or  re-processing  old  data""" 

7366  PRINT  USING  "SX.""  1  Plottina  data  on  Log-Log  """ 

7369  PRINT  USING  "6X.""  2  Plotting  data  on  Linear""" 

7372  PRINT  USING  "SX.""  3  Make  cross-plot  coefft  file""* 

737S  PRINT  USING  "SX,""  4  Re-circulate  water""" 

7378  PRINT  USING  "6X.""  5  Purge""" 

7381   PRINT  USING  "SX.""  6  T-Drop  correction""" 

7384  PRINT  USING  "6X.""  7  Print  Uo  File""" 

7387  PRINT  USING  "SX.""  8  Modify  X-Y  file""" 

7390   PRINT  USING  "6X,""  9  Move""" 

7393  PRINT  USING  "6X.""10  Comb""" 

7396  PRINT  USING  "6X.""11  Fixup""" 

7399   INPUT  Idp 

7402   IF  Idp-0  THEN  CALL  Main 

740S   IF  Idp-1  THEN  CALL  Plot 
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7408  IF  Idp-2  THEN  CALL  Plin 

741 1  IF  Idp-3  THEN  CALL  Coef 

7414  IF  Idp-4  THEN  CALL  Main 

7417  IF  Idp-5  THEN  CALL  Purg 

7420  IF  Idp-6  THEN  CALL  Tdcn 

7423  IF  Idp-7  THEN  CALL  Uoprt 

7426  IF  Idp-8  THEN  CALL  Xymod 

7429  IF  Idp-9  THEN  CALL  Move 

7432  IF  Idp-10  THEN  CALL  Comb 

7435  IF  Idp-11  THEN  CALL  Fixup 

7438  SUBEND 

7441  SUB  Xymod 

7444  PRINTER  IS  1 

7447  BEEP 

7450  INPUT  'ENTER  FILE  NAME'.FileS 

7453  ASSIGN  SFilel  TO  File* 

7456  BEEP 

7459  INPUT  "ENTER  NUMBER  OF  X-Y  PAIRS", Np 

7462  BEEP 

7465  INPUT  "ENTER  NEW  FILE  NAME",File2S 

7468  CREATE  BOAT  File2S,S 

7471  ASSIGN  $Fiie2  TO  File2S 

7474  BEEP 

7477  INPUT  "ENTER  NUMBER  OF  X-Y  PAIRS  TO  BE  DELETED" .Ndel 

7480  IF  Ndei-0  THEN  7492 

7483  FOR  1-1  TO  Ndel 

7486  BEEP 

7489  INPUT  "ENTER  DATA  SET  NUMBER  TO  BE  DELETED" ,Nd( I) 

7492  NEXT  I 

7495  FOR  J-1  TO  Np 

7498  ENTER  9File1 :X.Y 

7501  FOR  1-1  TO  Ndel 

7504  IF  Nd(I)-J  THEN  7516 

7507  NEXT  I 

7510  OUTPUT  8File2;X,Y 

7513  PRINT  J.X.Y 

7516  NEXT  J 

7519  PRINTER  IS  701 

7522  ASSIGN  8File1  TO  * 

7525  ASSIGN  9Fiie2  TO  * 

7528  SUBEND 

7531  SUB  Move 

7534!  FILE  NAME:  MOVE 

7537! 

7540  DIM  Bop(66) ,A(66) .B(66) ,CC66),D(66) ,E(66) ,F(66),G(66).H(66) ,J(66) .K(66) ,L( 

66) ,M(66) 

7543  DIM  Told$(66)CI4],N(66),vr(66),Ir(66) 

7546  BEEP 

7549  INPUT  "OLD  FILE  TO  MOVE" .D2_f i leS 

7552  ASSIGN  *File2  TO  D2_fileS 

7555  ENTER  3File2:Nrun,DateS.Ldtc1 ,Ldtc2. Itt 

7558  FOR  1-1  TO  Nrun 

7561  ENTER  SFi le2:8op( I) , ToldS< I) 

7564  ENTER  *File2:rt(  I)  ,BU)  ,C(I)  ,D(  I)  ,E(  I)  ,F(  I)  ,G(  I)  ,H(  I) ,  J(  I)  ,K(  I)  ,L(  I)  ,M<  I)  ,N 

(I) 

7567  ENTER  »Fi le2:Vr( I) ,Ir( I) 

7570  NEXT  I 

7573  ASSIGN  ?File2  TO  * 

7576  BEEP 

7579  INPUT  "SHIFT  DISK  AND  HIT  CONTINUE". Ok 
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7532  BEE? 

7535  INPUT  'INPUT  BDAT  SIZE", Size 

7588  CREATE  BDAT  D2_f  ileS.Sue 

7591  ASSIGN  ?FileI  TO  D2_fUeS 

7594  OUTPUT  «Fi  lei :Nrun.DateS.Ldtc1 ,Ldtc2. Itt 

7597  FOR  1-1  TO  Nrun 

7500  OUTPUT  9File1:Bop<I).ToldS(I) 

7603  OUTPUT  aFilet:A(I),3(I),C(I).D(I).E(I),F(I),G(I),H(I),J(I),!<(I).L(I).M(I) 

N(I) 

7806  OUTPUT  SFilel ;Vr( I) , Ir( I) 

7609  NEXT  I 

7612  ASSIGN  iFilel  TO  * 

7615!  RENAME  "TEST"  TO  D2_fileS 

7613  BEEP  2000.. 2 
7621  BEEP  4000.. 2 
7624  BEEP  4000,. 2 

7627  PRINT  'DATA  FILE  MOVED" 

7630  SUBEND 
7633  SUB  Comb 

7636!  FILE  NAME:  COMB 

7639! 

7642  DIM  ErafU2) 

7645  BEEP 

7648  INPUT  'OLD  FILE  TO  FIXUP" ,D2_f i leS 

7651  ASSIGN  8File2  TO  D2_fileS 

7654  D1  fileS-"TEST" 

7657  CREATE  BDAT  D1_file$.30 

7660  ASSIGN  ?Filel  TO  D1  fileS 

7663  ENTER  AFile2:Nrun,DateS,Ldtc1 .Ldtc2. I tt 

7666  IF  K-0  THEN  OUTPUT  »Fi lei ;Nrun,DateS.Ldtcl ,Ldtc2, Itt 

7669  FOR  1-1  TO  Nrun 

7672  ENTER  3File2:8op. ToldS.Emt ( *) ,Vr ,Ir 

7675  OUTPUT  9Fi lei :Bop, ToldS.Emt (*) ,Vr . Ir 

7678  NEXT  I 

7631  ASSIGN  §FLle2  TO  * 

7634!  RENAME  'TEST"  TO  D2_f:le$ 

7687  BEEP  2000.. 2 

7690  BEEP  4000.. 2 

7693  BEEP  4000,. 2 

7696  BEEP 

7699  INPUT  "UANT  TO  ADD  ANOTHER  FILE  ( 1-Y,0-N)?",Oka 

7702  IF  Qka-1  THEN 

7705  K-l 

7708  BEEP 

7711  INPUT  "GIVE  NEW  FILE  NAME",Nfile$ 

7714  ASSIGN  »File2  TO  NfileS 

7717  GOTO  7663 

7720  END  IF 

7723  ASSIGN  9File2  TO  * 

7726  SUBEND 

7729  SUB  Fixup 

7732!  FILE:  FIXUP 

7735!  DATE:  February  18,  1986 

7738! 

7741  DIM  Emt(12) 

7744  BEEP 

7747   INPUT  "OLD  FILE  TO  FIXUP" ,D2  fi ieS 

7750  ASSIGN  8File2  TO  D2_file$ 

7753  D1    file$«"TEST" 

7756  CREATE   BDAT   D1_fileS,20 

7759  ASSIGN  SFilel    TO   Dl    file* 
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7762  ENTER  *Fi  le2:Nrun.DateS.l_dtc1  .Ldtc2.  Itt 

7765  Nruno-27 

77S3  OUTPUT  ©File! :Nruno.DateS,Ldtc1 ,Ldtc2.Itt 

7771  FOR  1-1  TO  Nrun 

7774  ENTER  3FUe2:Bop,  ToldS,Emf  (*)  .  Vr.Ir 

7777  IF  1-1  THEN  7783 

7780  OUTPUT  9Fi lei :8op,ToldS,Emt (*) ,Vr , Ir 

7783  NEXT  I 

7786  ASSIGN  3File2  TO  * 

7739  ASSIGN  9File1  TO  * 

7792!  RENAME  "TEST"  TO  D2_fileS 

7795  BEEP  2000.. 2 

7798  BEEP  4000.. 2 

7801  BEEP  4000,. 2 

7804  SUBEND 
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APPENDIX  C 
EXAMPLES  OF  REPRESENTATIVE  DATA  RUNS 

The       two    following      data    runs       (TXE222    and      GTB185)    are 

representative       samples    of    the    data    taken    over    the    course    of 

the    investigation.    The    first    data    run    for    the    Thermoexce 1 -E 

surface      was    for    a    decreasing    heat-flux    condition    and    an    oil 

concentration      of    three    percent.    The    second    data    run    was    for 

the      GEWA-T    surface      was    also      for       a      decreasing      heat-flux 

condition    but    with    an    oil    concentration    of    one    percent. 

Month,  date  ami    tinm  : 05 :0S: !  t : 4 1  : 20 

NOTE:  Program  name  :  DRP4 
Disk  numoer   ■  13 
Old  file  nans:  TXE222 

This  data  set  taken  on  :  02: 18: 13: 42:S8 
Tube  Numoer:    10 

Data  Set  Number  -    I   Bulk  Oil  7.  -  3.0     02:18:13:48:44 
TC  No:    1      2      3      4      5      S      7      8 

Temo  :  8.71  9.45   9.00   8.48  9.70   8.12   10.21    8.S9 

Tua  Tliqd  TUqd2   Tvapr   Psat    Tsump 

8.74  2.23    2.17    2. OS   -1.77   -18.9 

Thetab  Htube      Qdp 

S.S41  1.4S5E--04  9.5I9£*04 

Data  Set  Number  -   2   Bulk  Oil  7.   -      3.0  02: 18: 13:49:S7 

TC  No:     1      2      3  4  5  5      7      8 

Temo  :  3. SI  9.01    8.42  8.00    9. S3    8.29   10.33    8. SI 

Tua  Tliqa  Tltqd2   Tvapr  Psat    Tsump 

8.S5  2.24    2.22    2.09  -1.73   -18.3 

Thetao  Htuoe      Qdp 

S.319  1.S08E-04  9.S31E+04 

Data  Set  Numoer  -    3   Bulk  Oil  7.   -   3.0  02: 18: 13:S0:S  I 

TC  No:  1      2      3  4      5      S      7      8 

Temo  :  9.S8   8.38   8.33  7.39   9.34   8.07   10.22   8.S4 

Tua    Tliqd    TUqd2   Tvapr  Psat    Tsump 

8. S3  2.26    2.19    2.07  -1.74   -18.8 

Thetao  Htube      Qdp 

S.30S  I.5I1E+04  9.S26E-I-04 

Data  Set  Number  -    4   Bulk  Oil  7.  -  3.0     02:18:14:02:37 

TC  No:     1  2      3      4      S      6      7      8 

Temp  :  S.93  7.05   S.71    S.48  7.23   S.24    7. S3   S.SS 

Tua    Tliqd  Tliqd2   Tvapr   Psat    Tsump 

S.SS  2.2S  2.21    2.14   -1.73   -18. S 

Thetab  Htube      Qdp 

4.433  1.414E+04   S.2S9E+04 

Data  Set  Number  -   S   Bulk  Oil  7.   -   3.0  02:18:14:03:13 

TC  No:    1  2      3  4  5      S      7      3 

Temo  :  6.32  7.18   S.30  S.SS   7.35   6.2S   7.65   8.52 

Tua     Tliqd  Tliod2   Tvapr  Psat    Tsump 

6.71  2.28  2.25    2.18  -1.69   -18.4 

Thetab  Htuoe      Qdp 

4.444  1.408E+04   6.2S8E-04 

Data  Set  Number  -   6   Bulk  Oil  7.    -  3.0     02: 18:  14:  1 3:S6 

TC  No:    1      2      3      4  5      6      7      8 

Temo  :  5. SO   5.65   5.41    S.26  5.58   4.90   5.84   5.06 

Tua     Tliqd    Tliqd2   Tvapr   Psat  Tsump 

S.28  2.32    2.28    2.35   -1.66  -16.3 

Thetab  Htuoe      Qdp 

2.981  I.304E+04   3.388E*04 
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02:18:14: 14:29 
S  7  8 

4.88        5.32        5.02 


02:13:14:20:33 

S  7  8 

4.07        4. 58        4.06 


02:18:14:21:05 
6  7  8 

4.10        4.65        4. It 


02:18:14:26:17 
6  7  8 

3.51        3.85        3.47 


02:18:14:26:48 
6  7  8 

3.44       3.80        3.39 


02:18:14:33:11 
6  7  3 

2.98       3.14       2.30 


02:18:14:33:43 
6  7  8 

2.90        3.08        2.31 


02:18:14:40:43 
6  7  3 

2.71        2.79        2.61 


Data   Set   Number    -      IS      Bulk    Oil    7.   -  3.0  02:18:14:41:19 

TC  No:  1  2  3  4  5  5  7  8 

Temp    :        2.66        2.70        2.65        2.67  2.74        2.66        2.71        2.54 
Tua           TUqd        TUqd2      Tvapr      Psat         Tsump 
2.65        2.39          2.00          3.78      -1.77      -18.0 
The  tab      Htube               Qdp 

.458      1.145E-04     5.244E*03 


Data  Set 

Number 

7   Bulk  Oil  Z> 

3.0 

TC  No: 

i 

2      3      4' 

5 

Temp  : 

5^4S 

5.60   5.37   5.31 

5.61 

Tua 

Tligd 

TUqd2   Tvapr   Psat 

Tsumo 

S.26 

2.29 

2.26    2.35   -1.68 

-16.2 

Thetab 

Htube 

Qdp 

2.387 

1.299E+04  3.381E+04 

Data  Set 

Numoer 

8  Built  Oil  2  - 

3.0 

TC  No: 

1 

2      3      4 

5 

Temp  : 

4.25 

4.33   4.22   4.16 

4.49 

Tua 

TUqd 

TUqd2   Tvapr   Psat 

Tsump 

4.20 

2.24 

2.16    2.64   -1.77 

-15.3 

Thetab 

Htube 

Qdp 

1.398 

1.176E- 

•■04  2.349E+04 

Data  Set 

Numoer 

9  Bulk  Oil  %   - 

3.0 

TC  No: 

1 

2      3      4 

5 

Temp  : 

4.33 

4.41   4.26   4.23 

4.52 

Tua 

TUqd 

TUqd2   Tvapr   Psat 

Tsump 

4.25 

2.28 

2.23    2.71   -1.71 

-15.8 

Thetab 

Htube 

Qdp 

2.004 

1.169E-04  2.344£*04 

Data  Set 

Numoer 

-   10   Bulk  Oil  ■/.   ■ 

3.0 

TC  No: 

1 

2      3-4 

5 

Temp  : 

3.68 

3.68   3.59   3.54 

3.80 

Tua 

Tliqd 

TUqd2   Tvapr   Psat 

Tsump 

3.59 

2.33 

2.27    2.95   -1.66 

-16.4 

ThetaD 

Htube 

Qdp 

1.293 

1.137E- 

►04   1.469E*04 

Data  Set 

Number 

-   11   Bulk  OH  X   - 

3.0 

TC  No: 

•1 

2      3      4 

5 

Temp  : 

3.56 

3.61   3.55   3.51 

3.75 

Tua 

TUqd 

TUqd2   Tvapr   Psat 

Tsuqp 

3.53 

2.28 

2.21    2.34  -1.72 

-16.4 

Thetab 

Htube 

Qdp 

1.285 

1 . 1 42E 

♦04   1.467E*04 

Data  Set 

Number 

-   12   Bulk  Oil  7.  - 

3.0 

TC  No: 

1 

2      3      4 

5 

Temp  : 

3.06 

3.07   3.02   3.00 

3.12 

Tua 

riiqd 

TUqd2   Tvapr   Psat 

Tsump 

3.01 

2.43 

2.11    3.43   -1.69 

-17.1 

Thetab 

Htuoe 

Qdp 

.739 

1.1 20E 

♦04  3.283E*03 

• 

Data  Set 

;  Numoer 

-   13  Bulk  Oil  I   - 

3.0 

TC  No: 

1 

2      3      4 

5 

Temp  : 

3.00 

3.0!    2.94   2.92 

3.05 

Tua 

TUqd 

TUqd2   Tvapr   Psat 

Tsump 

2.34 

2.36 

2.05    3.41   -1.77 

-17.2 

Thetab 

Htuoe 

Qdp 

••  .738 

1 . 1 36E 

♦04  8.386E*03 

Data  Se< 

:  Number 

-   14   Bulk  Oil  7.  - 

3.0 

TC  No: 

1 

2      3      4 

5 

Temp  : 

2.70 

2.74   2.73   2.71 

2.30 

Tua 

TUqd 

TUqd2   Tvapr   Psat 

Tsump 

2.71 

2.41 

2.03    3.79   -1.75 

-17.3 

Thetab 

Htube 

Qdp 

.488 

I.078E-H 

14  5.257E+03 
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Month,  date  and  time  :OS:06: 1 1 : 47 : 44 

NOTE:  Program  name  :  DRP4 
Disk  nunoer   -  10 
Old  file  name:  GTB18S 

This  data  set  taken  on  :  I  1 :2S  :  I  1  :3S  :33 
Tube  Numoer:    9 

Data  Set  Numoer  -    I   Bulk  Oil  7.  -  1.0     11 : :2S:  IS : 46: 12 

TC  No:    1      2      3      4  S      S      7      8 

Temp  :  9. 27  3. SI    7.33   7.37  7.93   7.32   9.35   7.81 

Tua  Tliqd  Tligd2   Tvapr   Psat  Tsumo 

7.97  2.23    2.23    2.24   -1.73  -16.6 

Thetao  Htuoe      Qdp 

S.739  I.S81E+Q4  9.073E-04 

Data  Set  NuraDer  -    2   Sulk  Oil  7.  -  1.0      11  :25:  IS:  46: 48 

TC  No:  1      2      3      4  S      6      7      3 

Temo  :  3.30  8.SS   7.36   7.90  7.38   7.35   9.88   7.34 

Tua     Tligd  Tligd2   Tvapr   Psat  Tsumo 

3,00  2.28    2.27    2.28   -1.68  -16.6 

Thetao  Htuoe      Qdp 

S.728  1.S82E*04  9.062E+04 

Data  Set  Numoer  -   3  Bulk  Oil  7.  -  1.0     1 1 :2S:  1S:S4:0S 

TC  No:  1      2      3      4  S      6      7      8 

Temp  :  6.93  7.13   6.62   6.52  6.64   6. SO    7.78   6.S4 

Tua     Tligd  Jligd2   Tvapr   Psat  Tsump 

S.S6  2.21    2.22    2.77   -1.7S  -17.0 

Thetao  Htuoe      Odo 

4.448  1.348E*04  S.995E-04 

Data  Set  Numoer  -   4  Bulk  Oil  Z   -  1.0     1 1 :25: 1S:S4:S5 
TC  No:     1      2      3      4      S      S      7      8 

Temp  :  6.95   7. IS   6.51    6.62  5.52   5.49   7.77   6. SO 
Twa     Tligd    Tligd2   Tvapr   Psat    Tsump 

6.66  2.22    2.23    2.31   -1.74   -16.3 

Thetao  Htuoe       Qdp 

4.430  1.3S3E-04  5.993E-04 

Data  Set  Numoer  -   5   Bulk  Oil  Z   -  1.0      11:25:16:03:48 
TC  No:    1      2      3      4      5      6      7      8 

Temp  :  5.09   6.19   S.32   S.31  5.77   5.74   6.37   5.56 

Twa  Tligd    Tliad2   Tvapr   Psat    Tsump 

S.32  2.24    2.23    3.18   -1.73   -17.4 

Thetao  HtuOe      Qdp 

3.S83  1.036E-04   3.712E-04 

Data  Set  Numoer  -   6  Bulk  OH  Z  -  4.0     1 1 :2S:  16:04:22 

TCNo:1      2      3      4  S      6      7      8 

Temp  :    5.11  6.23   S.34   S.34  5.30   5.77   6.42   5.73 

Tua     Tligd  Tligd2   Tvapr   Psat    Tsump 

S.85   2.27    2.26    3.22  -1.59  -17.4 
ThetaD   Htube      Qdp 
3.S88   1.034E<-04  3.7T2E*04 

Oata  Set  Numoer  -    7   3ulk  Oil  7.  -  1.0      11:25:16:10:15 
TC  No:     1      2      3      4      S      5      7      3 

Temp  :    S.41  S.S2   S.19   S.32  S.2S   5. IS   S.S0   5.12 

Twa     Tligd  Tligd2   Tvapr   Psat    Tsump 

S.24   2.24    2.20    3.61   -1.74   -17.7 
Thetao   Htuoe      Qdp 
3.020   7.320E-03  2.211E-04 

Data  Set  Numoer  -   8   Bulk  Oil  7.  -  1.0      11:25:16:12:11 

TC  No:  1      2      3      4  S      6      7      8 

Temp  :  S.2S  5.32   4.98   S.13  S.06   4.95   5.29   4.90 

Tua     Tligd  Tligd2   Tvapr   Psat    Tsump 

S.J4  2.14    2.05    3.S6   -1.88  -17.7 

Thetao  Htube      Qdp 

2.949  7.480E+03   2.206E*04 

Data  Set  Numoer  -   9   Bulk  Oil  7.  -  1.0      11:25:16:17:23 

TC  No:  1      2      3      4      S      5      7      3 

Temp  :  S.05    5.14   4.31    S.05  4.34    4.81    4.95    4.78 
Tua     Tligd    Tligd2   Tvapr   Psat    Tsump 

4.30  2.25    2.22    3.4S   -1.73   -17.7 

Thetao  Htube      Qdp 

2.663  S.114EH-03   1 .362E+04 

Data  Set  Number  -   10   Bulk  Oil  Z   -  1.0      1  I :2S: 16: 18:04 

TC  No:  1      2      3      4      S      5      7      8 

Temp  :  5.10   5.20   4.33   S.09  4.99   4.85   5.00   4.85 
Tua     Tligd    Tligd2   Tvapr   P3at    Tsump 

4.95  2.30    2.27    3.44   -1.57   -17.6 

Thetab  HtuDe      Qdp 

2.657  S.I03E-03   1.3S6E+04 
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Data  Set  Number  -   11   Sulk  OH  X  -   t.O     H :2S: 16:24:S4g 

lemp°-  «!t5  4238   4^2  4?34    4. 75    4. SI    4.S2    4.S8 

Twa  Tligd  Tliqd2   Tvapr   Psat   Tsumo 

4.55  2.30    2.29    3. IS   -1  .56   -IS. 3 
Thetab   Htube      Qdp 

2. 360   3.273E-03   7.72SE-03 

Data  Set  Number  -   12  Bulk  Oil  Z  -   1.0     11:25:16:25:38 
ic  II..     i       ?      3      4      5      o      '       o 

Temp  i    4.76  4.38    4.42   4.34    4.74    4. SI    4. SO    4. S3 

Tya    Tligd  Tliqd2   Tvapr  Psat    Tsump 

4.56  2.30    2.29    3.13   -1.66   -16.0 
Thetab   Htuoe      Qdp 

2.360   3.264E+03   7.702E--03 

Data  Set  Number  -   13  talk  OH  Z  -   t.O  11:25:16:31:21 

rr  «*»•  1       7               3  4      5      5       <                « 

Temp  I  4.45    «7S4   4.09  4.60    4.39    4.34    4.36    4.28 

lua  Tliqd   Tliqd2   Tvapr  Psat    Tsump 

4.37  2.27    2.27    2.95  -1.69   -14.9 

Thetab  Htube      Qdp 

2.094  2.36IE*03  4.944E+03 

Data  Set  Numoer  -   14  Bulk  Oil  Z  -  1.0     11 :2S: 16:31 :S9 

rr  u..  ?      3      4  5      6      7  8 

Temp  i  4.46  47S5    4.12   4.61  4.43    4.38    4.38    4.32 

"T«a  Tliqd  Tliqd2   Tvapr  Psat  Tsump 

4.39  2.30    2.29    2.94  -1.66  -14.3 

Thetab  Htuoe      Qdp 

2.096  2.352E*03   4.931E*03 

Data  Set  Numoer  -   IS  Bulk  Oil  Z  -  1.0     1 1 :25: 16:37 :S3 

TC  No:  1      2      3      4  S      6      7      8 

Temp  :  4.19  4.29   3.34    4.38  4.07    4.06   4.00    3.99 

Tua    Tliqd  Tliqd2   Tvapr   Psat  Tsump 

4.09  2.23    2.23    2.33   -1.73  -14.8 

Thetab  Htuoe      Qdp 

1.363  1.582E*03  2.946E-03 

Data  Set  Number  -   16  Bulk  Oil  7.   -   1.0  11:25:16:38:28 

TCNo:1      2      3  4      S      6      7      8 

Temp  :    4.20    4.30   3.83  4.36    4.07    4.05    3.98    3.99 

Tim    TUqd    Tliqd2   Tvapr  Psat    Tsump 

4.09   2.23    2.22    2.82  -1.74   -14.3 
Thetab   Htube      Qdp 
1.362   1.S32E*03  2.946E»03 

Data  Set  Numoer  -   17   Bulk  Oil  Z  -   1.0  11:25:16:44:56 

TC  No:  -1  2      3  4      5      6      7      8 

Temp  :  3.97  4.04   3.57  4.17   3.77   3.31    3.67   3.74 

Twa  Tliqd  Tliqd2   Tvapr  Psat    Tsump 

3.35  2.19  2.19    2.71  -1.78   -IS. 2 

Thetao  Htube      Qdp 

1.S60  1.I46E*03   1.901E-03 

Data  Set  Number  -   18   Bulk  Oil  7.   -  1.0      11:25:16:45:33 

TC  No:    1  2      3      4  S      6      7      8 

Temp  :  3.97  4.04   3.68   4.19  3.78   3.82   3.66   3.72 

tua  Tliqd  Tliqd2   Tvapr   Psat    Tsump 

3.35  2.19  2.19    2.71   -1.77   -15.2 

Thetab  Htube      Qdp 

1.660  1.147E*03   1.904E+03 

Data  Set  Number  -   19   Bulk  Oil  7.  -  1.0      1  1  :25: 16:S1  :  19 

TC  No:    1      2      3      4  S      6      7      8 

Temp  :  3.69  3.73   3. SO   3.77  3.40   3.46   3.25   3.32 

Twa     Tliqd  Tliqd2   Tvapr   Psat  Tsump 

3. SI  2.19    2.18    2. 58   -1.78  -15.4 

Thetab  Htube      Qdp 

1.326  3.S73E*02   1 . I 37E+03 

Data  Set  Number  -   20   Bulk  Oil  7.   -   1.0  1  1  :25: 16:S3:  16 

TCNo:1      2      3  4      5      S      7      8 

Temp  :    3.77   3.81    3.57  3.84   3.46   3.S4   3.32   3.38 

Tua    Tliqd    Tliqd2   Tvapr  Psat    Tsump 

3. 58   2.25    2.25    2.62  -T .71   -15.7 
Thetab   Htube      Qdp 
1.335  8.630E+02   1.152E+03 
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APPENDIX  D 

UNCERTAINTY  ANALYSIS 

The  uncertainty  of  the  heat- transfer  coefficient  at  60 
kW/m2  and  8  kW/m2  for  data  runs  of  the  GEWA-T  tube  (GTB185 
with  one  percent  oil  concentration)  and  of  the  Thermoexcel- 
E  ( TXE222  with  three  percent  oil  concentration)  are 
analyzed  below.  The  method  of  analysis  is  based  on  Kline- 
McClintock  CRef.  20]  method  of  uncertainty  analysis. 
The  heat- transfer  coefficient  is: 

h  =  (qc  )/  (Tw  o  -  Ts  a  t  )  (D.  1  ) 

and 

Two  -  T9  at  =  Tu  ,  -CQc   In  (D2  /D,  )  ]/(2  7T  k  L)   -  TsaT     (D.2) 
where 

h  =  heat- transfer  coefficient 
qc  =  heat  flux  corrected  for  end  loss 
TM o  =  average  outer  tube  wall  temperature 
T9  a t  =  saturation  temperature 
TM i  =  average  inner  tube  wall  temperature 
Qc  =  heat  input  corrected  for  end  losses 
Di  =  tube  internal  diameter 
D2  =  tube  external  diameter 
k  =  tube  wall  thermal  conductivity 
L  =  length  of  tube's  heated  surface 
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To   begin  the   analysis,  a  dummy  variable  is  assigned  to  the 
conduction  term  as  follows: 

DELT  =CQC   In  (D*/D, )]/(  2  J  k  L)  (D.3) 

In    accordance   with      Kline    and   McClintock,   the 
uncertainty  of  the  heat- transfer  coefficient  is: 


5h  = 

h 


qc 
qc 


6tu  r 


(D.  4) 


Two  T  3    A    T 


5 DELT        1   2     + 


Tu    o         -       Ts    A   T 


_ _<5Tsat_ 

Tw   0         —        Ts   A    T 


The  uncertainty  for  the  conduction  term  is  estimated  since 
the  logarithmic  error  is  negligible  when  compared  to  the 
other  error  of  the  analysis.  The  uncertainty  of  the 
conduction   term  is  as  follows: 


Sdelt 

DELT 


Sqc' 

Qc 


5k 

k 


ft' 


where   Qc  =  qc   7T  D:   L 

The  uncertainty  of  this  term  is  as  follows 


(D.  5) 


(D.6) 


5qc 

Qc 


5qc 

qc 


5d2 

D2 


L 


(D.  7) 


Table  4  contains  a  listing  of  the  computed  terms  of 
equations  D.l  through  D.7  for  selected  data  runs  for  the 
GEWA-T  and  the  Thermoexce 1 -E  tubes  at  specified  heat-flux 
setting  of  60  kW/m2  and  8  kW/m2 .  The  algorithm  in  the  data- 
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TABLE  4 
UNCERTAINTY  ANALYSIS 


Fi le  Name  GTB185 

Heat  Flux    60  kW/m2      8  kW/m2 

Oi  1  %  1%  1% 


TXE222 
60  kW/m2       8  kW/m2 
3%  3% 


0.015 


0.015 


0.015 


0.015 


i.k 
k 


0.03 


0.03 


0.03 


0.03 


6L 


L 

0.0025 

0.0025 

0.0025 

0.0025 

DELT 

0.486 

0.064 

0.486 

0.  064 

5DELT 

DELT 

0.069 

0.520 

0.069 

0.526 

T„  ,   (  C) 

6.845 

4.678 

6.904 

2.964 

T«,  (  C) 

0.49 

0.065 

0.49 

0.065 

TSftT  (  C) 

2.215 

2.295 

2.265 

2.205 

TsAT   (K> 

0.  1 

0.  1 

0.  1 

0.  1 

5  Tm  i  _.. 

T«  o  ~Ts  A  T 

0.111 

0.028 

0.  Ill 

0.  088 

6  DELT 

TM  o  _Ts  A  T 

0.008 

0.014 

0.008 

0.046 

5T§  A  i 

T*  0  ~ Ts  A  T 

0.023 

0.042 

0.023 

0.  135 

*<lc_ 

qc 

0.02 

0.02 

0.02 

0.02 

5h 

h 

0.  115 

0.056 

0.  115 

0.  169 

h(W/m2  K) 

13480 

3273 

14140 

11200 
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reduction  program  used  this  approach  and  assumed  that  all 
the  eight  wall  thermocouples  were  used  to  calculate  the 
heat- transfer  coefficient. 
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